Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,165 G→C 7.7% intergenic (‑99/+496) lacZ ← / ← lacI pseudogene, truncated/DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273360,1650GC 7.7% 498.0 / 4.5 171intergenic (‑99/+496)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (63/92);  new base C (10/6);  total (73/98)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 178 total aligned reads displayed.

CGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAA  >  CP009273/360018‑360310
                                                                                                                                                   |                                                                                                                                                 
cGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATTATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTGGCAACg                                                                                                                                                   >  1:56290/1‑148 (MQ=38)
  cTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                                 >  2:246188/1‑148 (MQ=255)
  cTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAAAGtt                                                                                                                                                 >  1:1184998/1‑148 (MQ=255)
     gCCAGCGGCTTAACATCCCGCCCCACCATCCAGTGCTGGCGCTCCGATTTTACGCTGCCAAACTTGAGATTCCCGAGCCCCAATTCGGGGCTTTTTTTGGGGGGCGTCCTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaa                                                                                                                                              <  1:1702598/148‑1 (MQ=16)
     cgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaa                                                                                                                                              <  2:2352011‑M1/67‑1 (MQ=255)
      ccAGCGGCTTCCCAACCCGCGCCCCCAGCCAGGGCTGGCGCTCCGATGATAAGCTGTCAAACATGAGAATGCCGATCCCCAATTCCTGGCAGTTTAGGGCGGGGGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaaa                                                                                                                                             <  1:247989/148‑1 (MQ=25)
      ccAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaaa                                                                                                                                             <  1:3555071/148‑1 (MQ=255)
      ccAGCGGCTTACCAACCAGCGCCCCCATCCAGTGCAGGAGCTCCGATGATAAGCTTTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaaa                                                                                                                                             <  1:3555036/148‑1 (MQ=255)
        tcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGTTCAAata                                                                                                                                           <  1:827398‑M1/70‑2 (MQ=255)
        tcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATc                                                                                                                                           <  1:1496404‑M1/70‑1 (MQ=255)
         gCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                          >  2:2784220/1‑148 (MQ=255)
            cgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCCGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                                       <  2:2570966‑M1/74‑1 (MQ=255)
             gacaaccaccagatcaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                                                      >  2:2664637‑M1/74‑148 (MQ=255)
             cTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                                                      >  2:1340283/1‑148 (MQ=255)
             cTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCAAATCCGcgc                                                                                                                                      >  1:2579595/1‑146 (MQ=255)
              acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                                                                                   >  2:2661277‑M1/73‑134 (MQ=255)
              acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                                                                                   >  2:2661281‑M1/73‑134 (MQ=255)
              acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACg                                                                                                                                                   <  1:2661278‑M1/62‑1 (MQ=255)
              acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATccgcgcc                                                                                                                                     <  1:2867577‑M1/76‑4 (MQ=255)
               tACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGGCGTGGGTTCGCACCGGGCAAATCCGCTccc                                                                                                                                    >  1:1056232/1‑148 (MQ=37)
                aaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattTATTCCTGGCAGTTTATGGCGGGCGTGCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGGTGAAATGGGCTCCcg                                                                                                                                   <  2:945614‑M1/78‑1 (MQ=255)
                aaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                                                                                   <  2:1495498‑M1/78‑1 (MQ=255)
                   caacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTGCAAATCCGCTCCcggcg                                                                                                                                <  1:2016663‑M1/81‑1 (MQ=255)
                     ccAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                              <  2:3213774/148‑1 (MQ=255)
                     acagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                              <  1:69619‑M1/83‑1 (MQ=255)
                     acagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCGTGGCAACGTTGAAATGCGcgcccgggggc                                                                                                                              <  1:1922557‑M1/83‑5 (MQ=255)
                              ccaTCCCGGGCAGGAGCTCCGATGATAAGCTGACAAACATGAGAATGCCGATCCCCATGTCCGGGCAGTTTAGGGGGGGCGTCGGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                     <  1:3393060/148‑1 (MQ=18)
                                aTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                                                                                                   <  2:3132497/148‑1 (MQ=255)
                                aTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                                                                                                   <  2:100725/148‑1 (MQ=255)
                                 tCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACtg                                                                                                                  <  2:2804614/148‑2 (MQ=255)
                                 tCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACtg                                                                                                                  <  2:1032750/148‑2 (MQ=255)
                                      tGCAGGAGCTCCGATGATAAGCTGTAAACCATGAGAATTCCGAGCCCCAATTCCGGGCATTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                             <  1:3141219/148‑1 (MQ=18)
                                      cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTcaggtg                                                                                                             <  2:419431‑M1/100‑3 (MQ=255)
                                      cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                             >  1:250778‑M1/49‑148 (MQ=255)
                                        cAGGCGCTCCTATGATACGCTGACAAACATGACAATGCCCAGCCCAAATTCCTGGCTGTTTAGGGCGGGCGTCCGCCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                                                                           <  1:2103819/148‑1 (MQ=17)
                                         agtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCGCCCGGCGGATGTGTGCTACTCCGGCGAGc                                                                                                          <  1:2024324‑M1/103‑1 (MQ=255)
                                            aGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCAGCACCCGGCGGATTTGTCCTACTCAGGCGAGCGtt                                                                                                       >  2:375750/1‑148 (MQ=255)
                                             gCTCCGATGAGACGCTGTCAAACATCACAATGCCGACCCCCAATGCCGGGCAGTTTCGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                                      <  1:400536/148‑1 (MQ=14)
                                                ccGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCTTCAAATCCCCTCCCCGCGGATTTGTCCTAATCCGGAGAGCGTTAAcc                                                                                                   >  2:910435/1‑148 (MQ=18)
                                                   ctgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACCCTCCCGGCGGCTGTGTCCTACTCAGGAGAGCGGTCACCGac                                                                                                <  2:2685424‑M1/113‑1 (MQ=255)
                                                   ctgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACCTTCAAAAACGCTCCCGGCGGATTTGTCCTAATCAGGAGAGCGTTCACCGac                                                                                                <  1:2576470‑M1/113‑1 (MQ=255)
                                                    tgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGaca                                                                                               <  2:1694893‑M1/114‑1 (MQ=255)
                                                      atAAGCTGTCAAACGTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaa                                                                                             <  2:888125/148‑1 (MQ=255)
                                                      atAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAAAGGTCAAAAACGCACCCCCCGGCATAGTGCTACAAAGGTGAGCGTTCACCGacaaa                                                                                             >  1:2654892/1‑148 (MQ=18)
                                                      atAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCTTCAAATACGCACCCGGCGGAAGTGTGCTACTGAGGAGAGCGTTCACCGaccaa                                                                                             >  1:2654307/1‑148 (MQ=18)
                                                          tttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTCAAATACGCACCCCGCGGAATTTTACAAAGCAGGAGAGCGGTGAACGACAaacaac                                                                                         <  1:3512469‑M1/120‑1 (MQ=255)
                                                          gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTGGCAACCTTAAAATACGCTCCCGGCGGATTTGTGCTACTAAGGAGAGCGGTCACCGACAaaaaac                                                                                         >  1:543110/1‑148 (MQ=17)
                                                          gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAaacaac                                                                                         <  2:346117/148‑1 (MQ=255)
                                                           cTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTGAGGAGAGCGTTCACCGACAaaccacc                                                                                        >  2:3207991/1‑147 (MQ=255)
                                                           cTGTAAAACATGAGATTTCCGAGCCCCATTTCCTGGATGTTTAGGGCGGGCGTCCTCCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                                        <  1:3027736/148‑1 (MQ=14)
                                                                tttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCAAGTTCAAATACGCTCCCCGCGGATTTTTGCTACTCAGGAGAGCGTTCACCGACAAACAACACATaa                                                                                   <  1:2938113‑M1/126‑1 (MQ=255)
                                                                tttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                                   >  1:2595035‑M1/23‑148 (MQ=255)
                                                                 ttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                                  >  2:69619‑M1/22‑148 (MQ=255)
                                                                 aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTTGGTGGCCACCTTCAAAAAAGCACCCCGCGGCTTTGTCCGACGCAGGCGTGGGGTTACCGCACAACAACACATaaa                                                                                  >  1:1272402/1‑148 (MQ=12)
                                                                 aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGCGAGCGTTCACCGCCAAACAACAGATaaa                                                                                  >  2:2592743/1‑148 (MQ=255)
                                                                 aaCATGAGAATCCCGACCCCCAATTCCTGGATGTTTTGGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                                  <  1:127176/148‑1 (MQ=14)
                                                                   cATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGGCCAACTTGAAATACGCGCCCGCCGCGTTTTTCCCACATCGGCGAGCGGTGAACGACCCACAACCCATAAAAc                                                                                >  1:3347000/1‑148 (MQ=12)
                                                                    tttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAAAGTTCAAATACGCTACCGGCCGAATTTTGCGACAAAGGAGAGCGTTGACCCACAAACAACagaaaaaaag                                                                               <  1:3236504‑M1/130‑3 (MQ=255)
                                                                    tttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                                <  2:1490111‑M1/129‑1 (MQ=255)
                                                                    tttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                                >  1:1490111‑M1/19‑147 (MQ=255)
                                                                     ttgtaggctggtgaatgAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAAGAGATAAAACGa                                                                              >  2:2967113‑M1/18‑148 (MQ=255)
                                                                     tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTGGCGACCGTGAAAAACGCGCCCCGCCGCTGTTTCCGACTCAGGGGAGCGGTGAACGACACACAACACAAAAAAcca                                                                              >  1:3383868/1‑146 (MQ=9)
                                                                     tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACCGTAAAAAACCCGCCCGGCGGCGTTGTGCTACAAAGGGGGGGGGTAACCCACAAAAAACCAATAAAAAGa                                                                              >  1:2827843/1‑148 (MQ=9)
                                                                     tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACCACAGATACTACGa                                                                              >  2:2731927/1‑148 (MQ=255)
                                                                      ttttgccgggtgaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                                                                             >  2:550270‑M1/17‑148 (MQ=255)
                                                                      gagaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                                                                             <  2:1355560/148‑1 (MQ=255)
                                                                      gagaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                                                                             <  2:2759330/148‑1 (MQ=255)
                                                                      gagaAGTCCGATCCCCACTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                                                                             <  2:2374299/148‑1 (MQ=255)
                                                                       agaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaaa                                                                            <  1:2673883/148‑1 (MQ=255)
                                                                        atgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAg                                                                           <  2:120303‑M1/134‑1 (MQ=255)
                                                                         tgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACCAACAACAGATAAAACGTAAgg                                                                          <  2:2709827‑M1/135‑1 (MQ=255)
                                                                          aTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCGTCAAATACGCTCCCGGCGGATTTGTCCGACTCAGGAGAGCGGTCACCGACAAACAACACATAAAAACAAAgtc                                                                         >  1:550270/1‑146 (MQ=14)
                                                                           ttCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTTCAAATACGCTCCCGGCGGGTTTGTGCGACTCAGGCGTGCGTTGACCCACCAACAACAGATAAAAAGAAATGcc                                                                        >  1:1746952/1‑148 (MQ=14)
                                                                           ttCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                                        >  2:2642501/1‑148 (MQ=255)
                                                                           ttCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGGCAACGTTCAAATGCGCGACCGGCGGGTTTGTCCTACTCAGGAGTGCGTTAACCCACAAACCACAGATAAAAAGAAAGGcc                                                                        >  1:2166723/1‑148 (MQ=14)
                                                                               gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGt                                                                    >  2:160661‑M1/8‑148 (MQ=255)
                                                                               gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGt                                                                    <  2:1728108‑M1/141‑1 (MQ=255)
                                                                               gATCCCCAATTCCTGGCAGTTTATGGCGGGCCTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCAGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGt                                                                    <  2:645144/148‑1 (MQ=255)
                                                                                ggaattaaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGGCGCAACGGTCCAATCCGGTCCCGGCGGATTTGGGCTACTCCGGCGAGCGTTCCCCGACCCACCCCAGCTAAATCGAACGGCCCCGtg                                                                   >  2:2396201/7‑147 (MQ=2)
                                                                                aTCCCCAATTCCTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc                                                                   <  2:1850620/148‑1 (MQ=38)
                                                                                aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc                                                                   <  2:1850625/148‑1 (MQ=38)
                                                                                  aattaaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCCACGGTCAAAAACGCGCCCCGCGGCTTTTTACAACTAAGGAGTGCGTTCACCGACAAACAACAGAAAAAACCAAACGCCCAGTCtt                                                                 >  1:1202293/5‑148 (MQ=11)
                                                                                   cccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCttt                                                                <  1:569562/148‑1 (MQ=35)
                                                                                    ccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCg                                                                                                  <  1:2323430/113‑1 (MQ=35)
                                                                                         tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                          <  2:3366001/148‑1 (MQ=11)
                                                                                         tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACACAGGAGAGCGTTCACCGACCAACAACAGATAAAACCAAAGGCCCCGTCTTTCGACTg                                                          >  2:2280371/1‑148 (MQ=2)
                                                                                                    tagggcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGtt                                               <  2:2139919/145‑1 (MQ=1)
                                                                                                                    cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                               <  1:196716/148‑1 (MQ=11)
                                                                                                                     ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTTTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGTGCCTTTCGTTTTATTTGTTGCCTgggg                              >  2:583525/1‑148 (MQ=9)
                                                                                                                     ccGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGCTTTGTCCTACTCAGGAGAGCGGTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTGCGTTTTATTTGATGCCTgggg                              >  2:583490/1‑148 (MQ=9)
                                                                                                                         cACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                          >  2:3079892/1‑148 (MQ=34)
                                                                                                                             cTCCGGGCCGTTGCGTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACCAACAACAGATAAAACGAAAGGCCCAGTCTTTCGCCTGAGCCGTTCGTTTTATTTGATGGCGGgggaattggga                      >  2:1111089/1‑148 (MQ=11)
                                                                                                                                   gCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaat                >  2:270156/1‑148 (MQ=255)
                                                                                                                                      gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCt             >  2:1304029/1‑148 (MQ=255)
                                                                                                                                        tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTg           >  2:2499496/1‑148 (MQ=255)
                                                                                                                                            tCGCAACCTTCAAATCCGCGCCCGGCGGATTTGTCCTACTCCCGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGa       <  1:876573/148‑1 (MQ=255)
                                                                                                                                            tCGCAAAGTTCAAATCCGCTCCCGGCGGATTTGTCCGACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGa       <  1:875616/148‑1 (MQ=255)
                                                                                                                                                aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCGTGGGGAATTGGGAATTAATTCTTGAAGACGaa   >  2:2840428/1‑148 (MQ=255)
                                                                                                                                                aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTGAAGACGaa   >  2:2840457/1‑148 (MQ=255)
                                                                                                                                                 aCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTGAAGACGaaa  >  1:2372335/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAA  >  CP009273/360018‑360310

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: