Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,156 C→G 8.8% intergenic (‑90/+505) 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,1560CG 8.8% 486.5 / 7.5 175intergenic (‑90/+505)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (60/97);  new base G (14/4);  total (74/101)
Fisher's exact test for biased strand distribution p-value = 1.89e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 178 total aligned reads displayed.

GCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGA  >  CP009273/360009‑360302
                                                                                                                                                   |                                                                                                                                                  
cctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGc                                                                                                                                                    <  2:250778‑M1/53‑1 (MQ=255)
  tactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTtgcgt                                                                                                                                                  <  1:160661‑M1/55‑3 (MQ=255)
   cTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGcgtg                                                                                                                                                 >  1:811234/1‑145 (MQ=255)
     tcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGc                                                                                                                                               <  2:2837453‑M1/58‑1 (MQ=255)
       aggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCaa                                                                                                                                             <  2:2523501‑M1/60‑1 (MQ=255)
       aCCGCTTGCCAGCGGCTTACCATCCAGCGCCACCAACCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa                                                                                                                                             <  2:434780/148‑1 (MQ=255)
        ccGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGTTCGCCAc                                                                                                                                            >  1:104682/1‑148 (MQ=255)
         cGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATTATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTGGCAACg                                                                                                                                           >  1:56290/1‑148 (MQ=38)
          agagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCACCGt                                                                                                                                          <  2:196716‑M1/63‑1 (MQ=255)
           gagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgtgttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                         <  2:1650922‑M1/64‑1 (MQ=255)
           cTTGCCAGCGGCTTACCATCCAGCGCCACCATCCCGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                         <  2:2151021/148‑1 (MQ=255)
           cTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                         >  2:246188/1‑148 (MQ=255)
             tGCCAGCGGCTTACCATCCCGCGCCACCATCCAGGGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                                                                                                                                       <  1:1185586/148‑1 (MQ=255)
              gCCAGCGGCTTAACATCCCGCCCCACCATCCAGTGCTGGCGCTCCGATTTTACGCTGCCAAACTTGAGATTCCCGAGCCCCAATTCGGGGCTTTTTTTGGGGGGCGTCCTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaa                                                                                                                                      <  1:1702598/148‑1 (MQ=16)
              cgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaa                                                                                                                                      <  2:2352011‑M1/67‑1 (MQ=255)
               ccAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaaa                                                                                                                                     <  1:3555071/148‑1 (MQ=255)
               ccAGCGGCTTACCAACCAGCGCCCCCATCCAGTGCAGGAGCTCCGATGATAAGCTTTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaaa                                                                                                                                     <  1:3555036/148‑1 (MQ=255)
                  gCGGCTTACCATCCAGCGCCACCATCCAGTGCCGGAGCTCCGATGATCAGCTGTCAAACATGAGAATTCCGATCCCCCATTCCTGGCAGGTTATGGCGGGCGTCCTGCCCGCCACCCTCCGTGCCGTGGCGTCGCAACGTTCAAATcc                                                                                                                                  >  1:1764447/1‑148 (MQ=37)
                  gCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                  >  2:2784220/1‑148 (MQ=255)
                      gacaaccaccagatcaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                                              >  2:2664637‑M1/74‑148 (MQ=255)
                      cTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGGAGCAACGTTAAAATCCGcac                                                                                                                              >  1:2579581/1‑146 (MQ=37)
                      cTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                                              >  2:1340227/1‑148 (MQ=255)
                       acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                                                                           >  2:2661277‑M1/73‑134 (MQ=255)
                       acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                                                                           >  2:2661278‑M1/73‑134 (MQ=255)
                       acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                                                                           <  1:2661277‑M1/62‑1 (MQ=255)
                       acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                                                                           <  1:2661281‑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)
                             tCCACCCCCACCATCCAGTGCAGGAGCGCCGATTATAAGCTGAAAAACATGATAATTCCGATCCCCACTTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                                                                                                       <  1:859646/148‑1 (MQ=21)
                              ccAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                      <  2:3213774/148‑1 (MQ=255)
                              acagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                      <  1:69619‑M1/83‑1 (MQ=255)
                               cAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGGGGGGGGCCACGTTGAAATACGCGGCcggcggga                                                                                                                     >  1:1241140/1‑146 (MQ=37)
                                 gataacaccacaggcccagtcgttcgactgagcctgtcgatttatttgaggcctggggaagtAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                                                                   >  2:1840840‑M1/63‑148 (MQ=255)
                                      accaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGTCAACGATCAAAAACCCACCCGGCGGCtttttgcc                                                                                                              >  1:2598268/1‑143 (MQ=25)
                                       ccaTCCCGGGCAGGAGCTCCGATGATAAGCTGACAAACATGAGAATGCCGATCCCCATGTCCGGGCAGTTTAGGGGGGGCGTCGGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                             <  1:3393060/148‑1 (MQ=18)
                                         aTCCAGTGCAGGAGCGCCGATGAAAAGCTGACAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCCGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                                                                                           <  2:100744/148‑1 (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)
                                                  agtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCGCCCGGCGGATGTGTGCTACTCCGGCGAGc                                                                                                  <  1:2024324‑M1/103‑1 (MQ=255)
                                                    gAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                <  2:1342610/148‑1 (MQ=255)
                                                    gAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTCGCCACGTTGAAATGAACGCCACGCGGAGTTGTGCGTCTCAGGAGTGCGt                                                                                                >  1:964206/1‑148 (MQ=21)
                                                      gCTCCGATGAGACGCTGTCAAACATCACAATGCCGACCCCCAATGCCGGGCAGTTTCGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                              <  1:400536/148‑1 (MQ=14)
                                                       cTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                             <  2:172699/148‑1 (MQ=255)
                                                         cgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                           <  2:3380786‑M1/110‑1 (MQ=255)
                                                         ccGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCTTCAAATCCCCTCCCCGCGGATTTGTCCTAATCCGGAGAGCGTTAAcc                                                                                           >  2:910435/1‑148 (MQ=18)
                                                           tatcttAAGCTGTCAAAATTTAGTCTTCCGCTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGa                                                                                         <  2:3037884/143‑1 (MQ=14)
                                                            ctgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACCCTCCCGGCGGCTGTGTCCTACTCAGGAGAGCGGTCACCGac                                                                                        <  2:2685424‑M1/113‑1 (MQ=255)
                                                              gatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCCACGTGGAAATCCGCGCCCCGGGGTTTTGTGATACTAAGGAGACCGTTGACCGacca                                                                                      >  1:2594975/1‑146 (MQ=18)
                                                              gagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACCGTCAAAAACCCTCCCCGCCGATTTGTACAACTCAGGAAAGCGGTAACCCacaa                                                                                      <  2:3273557‑M1/115‑1 (MQ=255)
                                                               atAAGCTGTCAAACGTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaa                                                                                     <  2:888125/148‑1 (MQ=255)
                                                               atAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCTTCAAATACGCACCCGGCGGAAGTGTGCTACTGAGGAGAGCGTTCACCGaccaa                                                                                     >  1:2654307/1‑148 (MQ=18)
                                                               agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTCAAATACGCTCCCCGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaa                                                                                     <  1:2079073‑M1/116‑1 (MQ=255)
                                                                   tttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTCAAATACGCACCCCGCGGAATTTTACAAAGCAGGAGAGCGGTGAACGACAaacaac                                                                                 <  1:3512469‑M1/120‑1 (MQ=255)
                                                                   gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTGGCAACCTTAAAATACGCTCCCGGCGGATTTGTGCTACTAAGGAGAGCGGTCACCGACAaaaaac                                                                                 >  1:543110/1‑148 (MQ=17)
                                                                    cTGTAAAACATGAGATTTCCGAGCCCCATTTCCTGGATGTTTAGGGCGGGCGTCCTCCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                                <  1:3027736/148‑1 (MQ=14)
                                                                     tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                               <  2:3143520/148‑1 (MQ=255)
                                                                         tttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCAAGTTCAAATACGCTCCCCGCGGATTTTTGCTACTCAGGAGAGCGTTCACCGACAAACAACACATaa                                                                           <  1:2938113‑M1/126‑1 (MQ=255)
                                                                         aaaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTGGCAAAGTGCAAATCCGCGCCCGGCGGCTTGGTGCTACTCCGGGGCTGGTTGACCGCCCAAAAACAGAtta                                                                           >  1:2466093/1‑146 (MQ=17)
                                                                          ttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaa                                                                          >  2:69619‑M1/22‑148 (MQ=255)
                                                                          aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTTGGTGGCCACCTTCAAAAAAGCACCCCGCGGCTTTGTCCGACGCAGGCGTGGGGTTACCGCACAACAACACATaaa                                                                          >  1:1272402/1‑148 (MQ=12)
                                                                          aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGGAGCCACCTGAAAAAACCCGCCCCGCGGCTTTGGGCAACACAGGAGAGCGGGCACCCCCCAACAACACATaaa                                                                          >  1:1272338/1‑148 (MQ=12)
                                                                          aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGCGAGCGTTCACCGCCAAACAACAGATaaa                                                                          >  2:2592743/1‑148 (MQ=255)
                                                                            cATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGGCCAACTTGAAATACGCGCCCGCCGCGTTTTTCCCACATCGGCGAGCGGTGAACGACCCACAACCCATAAAAc                                                                        >  1:3347000/1‑148 (MQ=12)
                                                                             tttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAAAGTTCAAATACGCTACCGGCCGAATTTTGCGACAAAGGAGAGCGTTGACCCACAAACAACagaaaaaaag                                                                       <  1:3236504‑M1/130‑3 (MQ=255)
                                                                              tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTGGCGACCGTGAAAAACGCGCCCCGCCGCTGTTTCCGACTCAGGGGAGCGGTGAACGACACACAACACAAAAAAcca                                                                      >  1:3383868/1‑146 (MQ=9)
                                                                              tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACCGTAAAAAACCCGCCCGGCGGCGTTGTGCTACAAAGGGGGGGGGTAACCCACAAAAAACCAATAAAAAGa                                                                      >  1:2827843/1‑148 (MQ=9)
                                                                               ttttgccgggtgaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                                                                     >  2:550270‑M1/17‑148 (MQ=255)
                                                                               gagaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                                                                     <  2:2759330/148‑1 (MQ=255)
                                                                               gagaAGTCCGATCCCCACTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                                                                     <  2:2374299/148‑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)
                                                                                        gATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGt                                                            <  2:2810442/148‑1 (MQ=38)
                                                                                        gATCCCCAATTCCTGGCAGTTTATGGCGGGCCTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCAGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGt                                                            <  2:645144/148‑1 (MQ=255)
                                                                                         ggaattaaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGGCGCAACGGTCCAATCCGGTCCCGGCGGATTTGGGCTACTCCGGCGAGCGTTCCCCGACCCACCCCAGCTAAATCGAACGGCCCCGtg                                                           >  2:2396201/7‑147 (MQ=2)
                                                                                         aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc                                                           <  2:1850625/148‑1 (MQ=38)
                                                                                            cccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCttt                                                        <  1:569562/148‑1 (MQ=35)
                                                                                             ccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCg                                                                                          <  1:2323430/113‑1 (MQ=35)
                                                                                               aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                     <  2:370804/148‑1 (MQ=32)
                                                                                               aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTGCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACCAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                     >  2:1154023/1‑148 (MQ=17)
                                                                                                 ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                   >  2:2717495/1‑148 (MQ=31)
                                                                                                  tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACACAGGAGAGCGTTCACCGACCAACAACAGATAAAACCAAAGGCCCCGTCTTTCGACTg                                                  >  2:2280371/1‑148 (MQ=2)
                                                                                                                             cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                       <  1:196716/148‑1 (MQ=11)
                                                                                                                              ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTTTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGTGCCTTTCGTTTTATTTGTTGCCTgggg                      >  2:583525/1‑148 (MQ=9)
                                                                                                                              ccGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGCTTTGTCCTACTCAGGAGAGCGGTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTGCGTTTTATTTGATGCCTgggg                      >  2:583490/1‑148 (MQ=9)
                                                                                                                               cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGggga                     <  2:104682/148‑1 (MQ=32)
                                                                                                                                gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAACCAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                    <  2:1496404/148‑1 (MQ=17)
                                                                                                                                gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAACCAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                    <  1:1495498/148‑1 (MQ=17)
                                                                                                                                   aCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAAAAGATAAAAAGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGGTTTATTTGTTGCCTGgggcattg                 >  2:3090827/1‑148 (MQ=11)
                                                                                                                                            gCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaat        >  2:270156/1‑148 (MQ=255)
                                                                                                                                            gCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAAGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTTATGCCTGGGGAATTGGGaattaat        >  2:1375124/1‑148 (MQ=255)
                                                                                                                                               gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCt     >  2:1304029/1‑148 (MQ=255)
                                                                                                                                                 tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTg   >  2:2499496/1‑148 (MQ=255)
                                                                                                                                                  gCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgg                <  1:250416/134‑1 (MQ=37)
                                                                                                                                                   |                                                                                                                                                  
GCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGA  >  CP009273/360009‑360302

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: