Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,173 C→A 5.4% intergenic (‑107/+488) 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,1730CA 5.4% 727.6 / 2.8 256intergenic (‑107/+488)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (73/157);  new base A (16/10);  total (89/167)
Fisher's exact test for biased strand distribution p-value = 4.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 258 total aligned reads displayed.

CAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGA  >  CP009273/360025‑360321
                                                                                                                                                    |                                                                                                                                                    
cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                                     <  2:516205/150‑1 (MQ=255)
   accgacaaccaacagagaaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                                                  >  2:1505435‑M1/77‑150 (MQ=255)
   accgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                                                  >  2:3915599‑M1/77‑150 (MQ=255)
    ggCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACCGTCAAATCCGccc                                                                                                                                                 >  1:3390425/1‑148 (MQ=255)
     cgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTcc                                                                                                                                                <  2:3132888‑M1/76‑1 (MQ=255)
     cgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCGACGTTCAAGTCCGCTcc                                                                                                                                                <  2:3520162‑M1/76‑1 (MQ=255)
       ttACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAACCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                                                                                              <  1:795287/150‑1 (MQ=255)
        tACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcgg                                                                                                                                             <  2:2316712/150‑1 (MQ=255)
        tACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcgg                                                                                                                                             <  2:4769268/150‑1 (MQ=255)
        caaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcgg                                                                                                                                             <  1:3243439‑M1/79‑1 (MQ=255)
         aaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttattggatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc                                                                                                                                            <  2:587548‑M1/80‑1 (MQ=255)
         aCCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc                                                                                                                                            >  2:2690011/1‑150 (MQ=255)
          ccATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGGGCTTCGCAACGTTCAAATCCCCCCCcggcg                                                                                                                                           >  2:4552647/1‑150 (MQ=255)
          aacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggacttAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCACCGGTCCAGTCCGCGCCcggcg                                                                                                                                           <  2:4337781‑M1/81‑1 (MQ=255)
               cAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                                                                                      <  2:3195400/150‑1 (MQ=255)
               cAGCGCCACCATCCAGAGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                                                                                      <  2:4726415/150‑1 (MQ=255)
                 gcgcCACCATCCAGTCCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGt                                                                                                                                    <  2:3365072/150‑1 (MQ=255)
                  ataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTGCCAATCCGCACCCGGCGGATTTGTc                                                                                                                                   <  2:1867763‑M1/89‑1 (MQ=255)
                  ataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                                   <  1:1945431‑M1/89‑1 (MQ=255)
                   gcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                                                                                  >  2:3546023/1‑150 (MQ=255)
                     caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                                <  1:1992283/150‑1 (MQ=255)
                     caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGGCCAAGTTAAAATCCGCGCCCGGCGGGTTTGTccaa                                                                                                                                >  1:1424823/1‑148 (MQ=25)
                     aaacgaaaggcccagtctttcgactgagcctttcgatttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                                >  2:1569438‑M1/59‑150 (MQ=255)
                      accaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTAAAATACGCTACCGGCGGATTTGTCCTaa                                                                                                                               >  2:3025782/1‑149 (MQ=255)
                      aacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                                               >  1:2294506‑M1/58‑150 (MQ=255)
                       ccaTCCTGTGCTGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCTGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                                                                                                              <  1:4751483/150‑1 (MQ=255)
                         aTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCTATCCCCAATTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                                            <  1:4413347/150‑1 (MQ=255)
                          tCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATGCGTGGCTGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAg                                                                                                                           <  2:16174/150‑1 (MQ=255)
                            aggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGga                                                                                                                         >  1:4436708‑M1/52‑150 (MQ=255)
                            aggcccagtctttcgactgagcctttcgttgtatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGGTGAAATCCGCTCCCGGCGGATTTGTCCTAATCAGga                                                                                                                         <  2:1277415‑M1/99‑1 (MQ=255)
                             ggcgccgtctttcgactgagcctttggttttattgtctgcgcggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                                        >  2:3721728‑M1/51‑150 (MQ=255)
                             aGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCATGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                                        >  2:2656127/1‑150 (MQ=255)
                               tGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTTCGTAGCAACGGTCAAATACCCTACCCGCGGATTTTTGCTACACAGGagag                                                                                                                      >  1:1062941/1‑150 (MQ=21)
                                gCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                                     <  2:4851640/150‑1 (MQ=255)
                                ccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGGTCAACACCCCACCCCGCGGCTTTGTGCTACACAGGAGAGc                                                                                                                     <  1:1672364‑M1/103‑1 (MQ=255)
                                ccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                                     <  2:1150188‑M1/103‑1 (MQ=255)
                                ccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCAGCCCGCCACCCTACGGGCCGTTGCTTAGCAACCTTCAAATACGCTCCCGGCGGCTTTGTCCTACTCAGGAGAGc                                                                                                                     <  1:2606558‑M1/103‑1 (MQ=255)
                                 cagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTGGCAACCGTCAAATCCGCTCCCCGCGGATTTGTACtactcaggtgtgcg                                                                                                                    <  1:3540460‑M1/104‑7 (MQ=255)
                                   gtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAAGTCAAATACGCTACCCGCCGCTTTGTCCAAAACACGAGAGCGtt                                                                                                                  <  1:2743732‑M1/106‑1 (MQ=255)
                                     ctaactacttaccctatctttttatttgacgcctgggtaataAATCCCTGGCAGTTTATGGCGGCCGTCCTCCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                                                >  1:3025782‑M1/43‑150 (MQ=255)
                                            ctgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGGGCGTTCACCGacaa                                                                                                         <  2:605682‑M1/115‑1 (MQ=255)
                                              gagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                                       >  2:4011684‑M1/34‑150 (MQ=255)
                                               agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaaca                                                                                                      >  1:289581‑M1/33‑150 (MQ=255)
                                               agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaaca                                                                                                      <  2:289581‑M1/118‑1 (MQ=255)
                                                gcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGaccaacaa                                                                                                     <  2:2593597‑M1/119‑1 (MQ=255)
                                                  agcttaaaaaCATTCGATTTCCTATCCCCATTTCCTGGCTGTTTAGGGGGGGCGTCCTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACGAacaaca                                                                                                   <  1:915556/143‑1 (MQ=11)
                                                   tttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCTGCAAATACGCTACCGGCGGCTTTTTCCAACTCAGGAGAGCGGTGACCGACAAACAACAg                                                                                                  <  1:2404236‑M1/122‑1 (MQ=255)
                                                   gCTGTAAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                                  <  2:2074127/150‑1 (MQ=255)
                                                    cTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGa                                                                                                 <  2:4277290/150‑1 (MQ=255)
                                                     tgacAAACATGAGAATTCCGATCACCAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGAGTTCACCGACAAACAACAGAt                                                                                                <  2:4778547/147‑1 (MQ=255)
                                                     tGTCAAACTGGAGAATTCCGATCCCCAATTCCTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                                                <  1:1466432/150‑1 (MQ=255)
                                                     tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTTCTTCGCAACGGGAAAATCCGCACCCGGCGGATTTGTCCGACAAAGGCGAGCGTTGACCGACAAAAAAAAGAt                                                                                                >  1:753565/1‑150 (MQ=17)
                                                       gttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCACTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                                              >  2:2728159‑M1/25‑150 (MQ=255)
                                                         aaaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                                            <  2:4241469/150‑1 (MQ=255)
                                                          ttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTAAAATACGCTCCCCGCGGAGTTGTGCTACAAAGGAGTGCGTTCACCCACAAACAACAcaaaaaaa                                                                                           <  1:650219‑M1/129‑2 (MQ=255)
                                                          ttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAAAGTTAAAATACCCACCCCGCGGATTTGTGCTACACAGGAGAGCGGTCACCGACAAACAAACGAtaaaaa                                                                                           <  1:3984247‑M1/129‑2 (MQ=255)
                                                          aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                                           <  2:3061135/150‑1 (MQ=255)
                                                           ccATGTGAATCCCGACCCCCATTTCGTGGCTTTTTTGGGGGGGGGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                                          <  1:56437/149‑1 (MQ=11)
                                                            cATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTGCCAAAGGTCAAATACGCACCCCGCGGATTTTTGCTAATGAGGAGATCGTTGACCGAAAAAAAACAGAAAAAAAGa                                                                                         >  1:1024547/1‑150 (MQ=17)
                                                               tgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTCAAATCCGCTCCCGGCGGATTTTTCCTACTCATGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAg                                                                                      <  1:1577983‑M1/134‑1 (MQ=255)
                                                               gagaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAg                                                                                      <  1:1542832/150‑1 (MQ=255)
                                                                 atgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACCAACAACAGATAAAACGAAAGGc                                                                                    <  2:1998998‑M1/136‑1 (MQ=255)
                                                                  tgcgtgggtaattAATTCGTGGCAGTTTATGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                                                   >  2:3251507‑M1/14‑150 (MQ=255)
                                                                   aTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                                                                                  <  2:3424895/150‑1 (MQ=255)
                                                                    ttCCGATCCCCAATTCCGGGCAGTTTTTGGGGGGCGTCAGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                                                 <  1:2779369/150‑1 (MQ=14)
                                                                      ccGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGt                                                                               <  2:902001/150‑1 (MQ=40)
                                                                        gcgccccAAGTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                                             <  2:3657406/147‑1 (MQ=16)
                                                                         aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAAAGTTAAAAAACGCTCCCGGCGGCTTTGTGCTAAACAGGAGAGCGGTACCCCACAAACAACAGATAAAAAGAAAGGCCCAGTCtt                                                                            >  2:1231585/1‑150 (MQ=11)
                                                                         aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                                            <  1:2164209/150‑1 (MQ=38)
                                                                         aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGGCAACGGTCAAATACGCTCCCGGCGGATTTGTACAACTCAGGTGAGCGTTCACCGACAAACCACAGATAAAACGAAAGGCCCAGTcgt                                                                            >  2:1128408/1‑148 (MQ=11)
                                                                          tCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTTCATCGCAAACAGAAAAACCGCTACCCGCGGATTTTTCCTAAACAGGGGAGCGTTCCCCGACAAACAAAACATAAAACGAAAGGCCCACTCatt                                                                           >  1:1540843/1‑150 (MQ=11)
                                                                           aattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                                          <  2:2805273‑M1/146‑1 (MQ=255)
                                                                              tacttCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCGACCCTCCGGGCCGTTGCGTCGCTAGGTGCAAATCCGCTGCCGGCGGTTTTGTCCTACTCTGGAGAGCGTTCACGGACAACCAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                                       <  2:253501/147‑1 (MQ=11)
                                                                              cAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                                       <  2:1866234/150‑1 (MQ=33)
                                                                              cAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                                       <  1:1874143/150‑1 (MQ=33)
                                                                               aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                                                                                      >  2:2706406/1‑102 (MQ=33)
                                                                                 ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                                    >  2:4641650/1‑150 (MQ=31)
                                                                                  tCCGGGCAGTTTCTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAg                                                                   <  2:4496041/150‑1 (MQ=11)
                                                                                            ctctgccgggctTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGttt                                                         <  1:2922103/144‑1 (MQ=0)
                                                                                                         ccggCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTg                                            <  1:2436062/147‑1 (MQ=1)
                                                                                                         cccgCCCCCCACCCTCCGGGCCGTCGCCTCGCAACGGACAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTg                                            <  2:1263996/147‑1 (MQ=1)
                                                                                                           ggCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                                          <  2:3983809/149‑1 (MQ=11)
                                                                                                             cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCGACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGggga                                        <  2:2687484/150‑1 (MQ=17)
                                                                                                              ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                                       <  2:1540992/150‑1 (MQ=32)
                                                                                                               ccccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAGCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                      <  1:2866340/148‑1 (MQ=18)
                                                                                                               cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                      <  1:1169433/150‑1 (MQ=33)
                                                                                                               cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                      <  1:3896899/150‑1 (MQ=33)
                                                                                                               cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCGTTCGTTTTATTTGATGCCGggggtct                                      >  2:1738563/1‑147 (MQ=9)
                                                                                                               cGCCACCCTCCGGGCCGTTGCTTCGCAAAGGTCAAATCCGCTCCCGGCGGATTTGTCCTACACAGGAGAGCGTTCACCGACAAACAAAAGATAAAACGAAAGGCCCAGTCTTTAGACTGATCCTTTCGTTTTATTTGATGCCTGgggaag                                      >  2:360929/1‑149 (MQ=11)
                                                                                                                    catcccGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCTGCGGATTTGTCCTACTCAGGAGTGCGTTCACCGACAAAAAACAGATAAAACGAAAGGCCCAGTCTTTCTACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattggga                                 <  1:4175548/146‑1 (MQ=11)
                                                                                                                           ggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaatt                          <  1:2713642/150‑1 (MQ=255)
                                                                                                                           ggCCGTTGCGTCGCAACGTACACCTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaatt                          <  1:1437857/150‑1 (MQ=255)
                                                                                                                               gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTg                      <  1:1127441/150‑1 (MQ=255)
                                                                                                                                ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGTATTAATTCTTGa                     >  2:409949/1‑150 (MQ=255)
                                                                                                                                   cTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                                      >  1:4846406/1‑130 (MQ=34)
                                                                                                                                     acGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACGCAGGAGAGCGTTCACCGACACACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCGGGGGAATTGGGAATTAATTAGTGAAGACg                >  1:996902/2‑150 (MQ=14)
                                                                                                                                      cGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAAGTGGGAATTAATTCTTGAAGACGa               >  1:41138/1‑150 (MQ=255)
                                                                                                                                        cAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTGGTTTTATTTGATGGCTGGGGAAGTGGGAAGTATTGCTTGATGACGata             >  2:3834377/1‑148 (MQ=14)
                                                                                                                                        cAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAAGTGGGAATTAATTCTTGAAGACGaaa             >  2:3772942/1‑150 (MQ=255)
                                                                                                                                          aCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAgg           <  1:4765679/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGACGAAAGGGCCTCGTGA  >  CP009273/360025‑360321

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: