Predicted mutation
evidence seq id position mutation annotation gene description
RA CP002416 3,132,643 C→T L152L (TTG→TTA)  Clo1313_2666 ← Hedgehog/intein hint domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0024163,132,6430CT98.5% 145.5 / NA 68L152L (TTG→TTA) Clo1313_2666Hedgehog/intein hint domain protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (36/31);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAAATCCTTCCAGAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCCAACCCGGTAATAACTACTGCGGCAGCCACTATCTTGGCAATGGATTTCCAATTCTCTTTACACCACTCCGCAACCGACTTTAAGCCGTCCTTG  >  CP002416/3132497‑3132736
                                                                                                                                                  |                                                                                             
gCAAATCCTTCCAGAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAAcc                                                                                           <  1:1402019/151‑1 (MQ=17)
 cAAATCCTTCCAGAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAAccc                                                                                          >  1:3284820/1‑151 (MQ=17)
 cAAATCCTTCCAGAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAAccc                                                                                          <  1:1651469/151‑1 (MQ=17)
 cAAATCCTGCCAGAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAACGCTCCTGCCAGTATGACTCCCAATTCCCCTCCTGTCAATGCCGCCGCTATCCCTAAccc                                                                                          <  1:3344775/151‑1 (MQ=11)
  aaaTCCTTCCAGAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCg                                                                                         >  1:891230/1‑151 (MQ=17)
  aaaTCCTTCCAGAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCg                                                                                         >  1:2963647/1‑151 (MQ=17)
  aaaTCCTTCCAGAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCg                                                                                         >  1:1401234/1‑151 (MQ=17)
  aaaTCCTTCCAGAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCg                                                                                         >  1:1218925/1‑151 (MQ=17)
   aaTCCTTCCAGAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCgg                                                                                        >  1:2784359/1‑151 (MQ=17)
    aTCCTTCCAGAAACGATCCTCCATTTAGCGCAGCGGCTATTCCTCCAACTGCCCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCCCCTGCCAGTATGACTCCCAATGCCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGt                                                                                       <  1:336493/151‑1 (MQ=11)
     tCCTTCCAGAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGta                                                                                      <  1:508835/151‑1 (MQ=17)
     tCCTTCCAGAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGta                                                                                      >  1:3480640/1‑151 (MQ=17)
     tCCTTCCAGAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGta                                                                                      <  1:2024456/151‑1 (MQ=17)
      ccTTCCAGAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGtaa                                                                                     >  1:418015/1‑151 (MQ=17)
      ccTTCCAGAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGtaa                                                                                     <  1:224046/151‑1 (MQ=17)
      ccTTCCAGAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGtaa                                                                                     >  1:2880147/1‑151 (MQ=17)
           cAGAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAact                                                                                <  1:6143/151‑1 (MQ=17)
           cAGAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAact                                                                                >  1:301666/1‑151 (MQ=17)
            aGAAACGTTCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAacta                                                                               <  1:2971580/151‑1 (MQ=12)
            aGAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAacta                                                                               >  1:2400728/1‑151 (MQ=17)
             gaacctaTCCTCCATTTGTCGCACCGGCTATTCCTCCACCTCCTCCCCCTATCAATCCACCGGCCAATGCTCCCCAGAAGGCTCCTGCCAGTATGCCTCCCAATTCCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAactac                                                                              <  1:217792/145‑1 (MQ=11)
             gAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAactac                                                                              >  1:2950089/1‑151 (MQ=17)
             gAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAactac                                                                              <  1:2063620/151‑1 (MQ=17)
              aaaCGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAactact                                                                             <  1:2211990/151‑1 (MQ=17)
              aaaCGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAactact                                                                             >  1:358983/1‑151 (MQ=17)
              aaaCGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAactact                                                                             <  1:509786/151‑1 (MQ=17)
               aaCGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTg                                                                            <  1:228613/151‑1 (MQ=17)
               aaCGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTg                                                                            >  1:2830876/1‑151 (MQ=17)
                 cttacctccATTTACCGCAGCGCCTATTCCACCAACTGCTCCCCCTATCAGTCCTCCCCCCAATGCTCCCCAGAATGCTCCTGCCAGTATGCCTCCCAATACCCCTCTTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCg                                                                          <  1:856691/147‑1 (MQ=11)
                 cGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCg                                                                          <  1:2924613/151‑1 (MQ=17)
                  gATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCgg                                                                         >  1:2706288/1‑151 (MQ=17)
                   aTCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGc                                                                        >  1:2382809/1‑151 (MQ=17)
                   aTCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGc                                                                        >  1:2394600/1‑151 (MQ=17)
                   aTCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATACTCCGGCCAATGCTCCCCAGAATGCCCCTGCCAGTATGACTCCCCATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGc                                                                        >  1:108296/1‑151 (MQ=11)
                    tcctccATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAACCCTCCGGCCAATGCTCCCCAAAATGCTCCTGCCAGTATGACTCCCAATGCCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCa                                                                       <  1:2874248/151‑1 (MQ=11)
                     cctccATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCCCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAg                                                                      <  1:2767828/151‑1 (MQ=12)
                       tccATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGcc                                                                    <  1:276453/151‑1 (MQ=17)
                        ccATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCa                                                                   >  1:574143/1‑151 (MQ=17)
                        ccATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCa                                                                   <  1:1757250/151‑1 (MQ=17)
                        ccATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCa                                                                   >  1:2012535/1‑151 (MQ=17)
                         cATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCAc                                                                  >  1:2480027/1‑151 (MQ=17)
                           tttATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTa                                                                >  1:1724716/1‑151 (MQ=17)
                            ttATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTAt                                                               >  1:2811645/1‑151 (MQ=17)
                            ttATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTAt                                                               >  1:2001098/1‑151 (MQ=17)
                             tATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAACGCTCCCCAGAATGCCCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTATc                                                              >  1:2011732/1‑151 (MQ=255)
                              aTCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTATCt                                                             <  1:768679/151‑1 (MQ=17)
                              aTCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAACGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTATCt                                                             >  1:1065698/1‑151 (MQ=12)
                               tCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAATTACTGCGGCAGCCCCTATCtt                                                            >  1:2023191/1‑151 (MQ=255)
                                cGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTATCTTg                                                           >  1:1170034/1‑151 (MQ=17)
                                   aGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTATCTTGGCa                                                        <  1:427450/151‑1 (MQ=17)
                                   aGCGGCTATTCCTCCAACTGCTCCCCCTATCAACCCCCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTATCTTGGCa                                                        <  1:1642590/151‑1 (MQ=255)
                                     cGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTATCTTGGCAAt                                                      <  1:3373346/151‑1 (MQ=17)
                                      ggCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTATCTTGGCAATg                                                     <  1:850350/151‑1 (MQ=17)
                                      ggCTATTCCCCCAACTGCTCCCCCTATCACCCCTCCGGCCAATGCTCCCCAGAATGCCCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAGCCCGGTAATAACTACTGCGGCAGCCACTATCTTGGCAATg                                                     <  1:2151551/151‑1 (MQ=11)
                                         tATTCCTCCACCTGCCCCCCCTATCAACCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTATCTTGGCAATGGAt                                                  <  1:3184096/151‑1 (MQ=11)
                                          aTTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTAGGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTATCTTGGCAATGGAtt                                                 >  1:3294120/1‑151 (MQ=12)
                                           ttCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTATCTTGGCAATGGAttt                                                >  1:2300670/1‑151 (MQ=17)
                                             cctccAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTATCTTGGCAATGGATTTcc                                              >  1:1016291/1‑151 (MQ=17)
                                             cctccAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTATCTTGGCAATGGATTTcc                                              >  1:1264741/1‑151 (MQ=17)
                                             cctccAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTATCTTGGCAATGGATTTcc                                              >  1:1287629/1‑151 (MQ=17)
                                              ctccAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTATCTTGGCAATGGATTTCCa                                             >  1:47206/1‑151 (MQ=17)
                                              ctccAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTATCTTGGCAATGGATTTCCa                                             >  1:16885/1‑151 (MQ=17)
                                                   ccTGCTCCCCCTATCAACCCTCCGGCCAGCGCTCCCCAGAATGCCCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAGCCCGGTAATAACTACTGCGGCAGCCACTATCTTGGCAATGGATTTCCAATtct                                        <  1:2029297/150‑1 (MQ=11)
                                                    cTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTATCTTGGCAATGGATTTCCAATtctc                                       >  1:2090872/1‑151 (MQ=14)
                                                     tGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTATCTTGGCAATGGATTTCCAATTctct                                      <  1:578707/151‑1 (MQ=11)
                                                     tGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTATCTTGGCAATGGATTTCCAATTctct                                      <  1:1972058/151‑1 (MQ=11)
                                                     tGCTCCCCCTAACAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATCCCCCTCCTGTCAATGCCGCCGCTATCCCTAACCCGGTAATAACTACTGCGGCAGCCACTATCTTGGCAATGGATTTCCAATTctct                                      <  1:47205/151‑1 (MQ=2)
                                                                                         gAATGCTCCTGCCAGTATGACTCCCAATATCCCTCCTGTCAATGTCTTTGCTATCCCCTACCCGGTAATAACTACTGCGCTAGTTACTACCTTGGCTATGGATTTCCAATTCTCTTTACACCACTCAGCAACCGACTTCAAGCCGTccctg  >  1:2324597/1‑148 (MQ=2)
                                                                                                                                                  |                                                                                             
GCAAATCCTTCCAGAAACGATCCTCCATTTATCGCAGCGGCTATTCCTCCAACTGCTCCCCCTATCAATCCTCCGGCCAATGCTCCCCAGAATGCTCCTGCCAGTATGACTCCCAATACCCCTCCTGTCAATGCCGCCGCTATCCCCAACCCGGTAATAACTACTGCGGCAGCCACTATCTTGGCAATGGATTTCCAATTCTCTTTACACCACTCCGCAACCGACTTTAAGCCGTCCTTG  >  CP002416/3132497‑3132736

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑