Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,208,550 T→A 45.4% K95* (AAA→TAA)  mug ← G/U mismatch‑specific DNA glycosylase; xanthine DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,208,5500TA45.4% 87.8 / inf 238K95* (AAA→TAA) mugG/U mismatch‑specific DNA glycosylase; xanthine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base T (69/61);  new base A (58/50);  total (127/111)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 243 total aligned reads displayed.

TAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAA  >  CP009273/3208401‑3208696
                                                                                                                                                     |                                                                                                                                                  
 aaaCCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGGTAATGTTCATttgtt                                                                                                                                                   >  2:4081131/1‑150 (MQ=255)
  aaCCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAAtttatt                                                                                                                                                  <  1:1051929/150‑1 (MQ=255)
   aCCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTc                                                                                                                                                 <  1:1969120/150‑1 (MQ=255)
   aCCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTc                                                                                                                                                 <  1:1969046/150‑1 (MQ=255)
   aCCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTc                                                                                                                                                 <  1:2374464/150‑1 (MQ=255)
     cGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCaa                                                                                                                                               >  1:713489/1‑150 (MQ=255)
      gCTGGGATTTGGCAGCAACCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAAt                                                                                                                                              >  2:530321/1‑150 (MQ=255)
       gtGGGTTGTGGCAGCTCCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCTCCGCGCTGGCTGAATGCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCTATc                                                                                                                                             <  2:1073968/149‑1 (MQ=255)
       cTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATc                                                                                                                                             <  2:2237790/150‑1 (MQ=255)
       cTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATc                                                                                                                                             <  2:3350446/150‑1 (MQ=255)
         gggATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAg                                                                                                                                           >  1:1880856/1‑150 (MQ=255)
              ttGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTAc                                                                                                                                      <  1:1542376/150‑1 (MQ=255)
              ttGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTAc                                                                                                                                      >  2:2947548/1‑150 (MQ=255)
              ttGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTAc                                                                                                                                      <  1:3173123/150‑1 (MQ=255)
                  cagcaCCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTAcgccc                                                                                                                                  <  1:2207447/150‑1 (MQ=255)
                      aCCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCt                                                                                                                              >  2:2545868/1‑150 (MQ=255)
                      aCCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCt                                                                                                                              >  2:246973/1‑150 (MQ=255)
                           aaTCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCgtg                                                                                                                         <  2:414553/150‑1 (MQ=255)
                              cTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAg                                                                                                                      >  2:1650456/1‑150 (MQ=255)
                               tGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGc                                                                                                                     >  1:1982622/1‑150 (MQ=255)
                                  gTCGAACCACTGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTcc                                                                                                                  <  2:1919091/150‑1 (MQ=255)
                                  gTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTcc                                                                                                                  >  1:4054046/1‑150 (MQ=255)
                                  gTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTcc                                                                                                                  <  2:1359302/150‑1 (MQ=255)
                                   tCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCt                                                                                                                 <  1:2615471/150‑1 (MQ=255)
                                           aTGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGaa                                                                                                         >  1:2888229/1‑150 (MQ=255)
                                               tGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACtt                                                                                                     <  2:3463929/150‑1 (MQ=255)
                                               tGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACtt                                                                                                     <  2:3463985/150‑1 (MQ=255)
                                                 aGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCa                                                                                                   >  1:3855730/1‑150 (MQ=255)
                                                 aGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCa                                                                                                   <  2:2793258/150‑1 (MQ=255)
                                                       tgtttCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCt                                                                                             <  2:1490287/150‑1 (MQ=255)
                                                           tCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCa                                                                                         >  1:301012/1‑150 (MQ=255)
                                                              cccACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCg                                                                                      >  2:60864/1‑150 (MQ=255)
                                                                cACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGtt                                                                                    <  2:2831476/150‑1 (MQ=255)
                                                                cACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGtt                                                                                    <  2:544656/150‑1 (MQ=255)
                                                                 aCTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTg                                                                                   >  1:2892591/1‑150 (MQ=255)
                                                                       cACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGACACTTCATTGGCTTGCACCGTTGGACGGt                                                                             >  2:4397525/1‑150 (MQ=255)
                                                                       cACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGt                                                                             >  2:4397616/1‑150 (MQ=255)
                                                                       cACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGt                                                                             >  2:1181544/1‑150 (MQ=255)
                                                                         ccGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCt                                                                           >  1:75031/1‑150 (MQ=255)
                                                                           gcgcTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTAc                                                                         <  2:891235/150‑1 (MQ=255)
                                                                           gcgcTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTAc                                                                         >  2:588916/1‑150 (MQ=255)
                                                                             gctggctgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCa                                                                       >  2:211330/1‑150 (MQ=255)
                                                                             gctggctgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCa                                                                       >  2:210412/1‑150 (MQ=255)
                                                                              ctggctgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAg                                                                      <  2:2594808/150‑1 (MQ=255)
                                                                              ctggctgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAg                                                                      <  2:4532329/150‑1 (MQ=255)
                                                                              ctggctgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAg                                                                      <  2:4532276/150‑1 (MQ=255)
                                                                                 gctgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGttt                                                                   <  1:3285816/150‑1 (MQ=255)
                                                                                 gctgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGttt                                                                   <  1:3095067/150‑1 (MQ=255)
                                                                                 gctgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGttt                                                                   <  1:4199908/150‑1 (MQ=255)
                                                                                  ctgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTg                                                                  <  2:433611/150‑1 (MQ=255)
                                                                                   tgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTgg                                                                 >  2:1839244/1‑150 (MQ=255)
                                                                                   tgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTgg                                                                 >  2:1812608/1‑150 (MQ=255)
                                                                                   tgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTgg                                                                 >  1:904843/1‑150 (MQ=255)
                                                                                    gAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGt                                                                >  2:4483954/1‑150 (MQ=255)
                                                                                    gAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGt                                                                >  1:4280005/1‑150 (MQ=255)
                                                                                     aaTCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTg                                                               >  1:2693260/1‑150 (MQ=255)
                                                                                       tCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGAc                                                             >  1:655665/1‑150 (MQ=255)
                                                                                       tCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGAc                                                             <  1:1177932/150‑1 (MQ=255)
                                                                                       tCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGAc                                                             >  1:2980154/1‑150 (MQ=255)
                                                                                       tCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGAc                                                             <  2:2784710/150‑1 (MQ=255)
                                                                                             ttCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCcaca                                                       >  1:107804/1‑150 (MQ=255)
                                                                                             ttCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCcaca                                                       >  1:4055210/1‑150 (MQ=255)
                                                                                                atatGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACg                                                    >  2:1316005/1‑150 (MQ=255)
                                                                                                atatGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACg                                                    <  1:802938/150‑1 (MQ=255)
                                                                                                   tGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATa                                                 >  2:340893/1‑150 (MQ=255)
                                                                                                       tGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATcc                                             <  2:3715247/150‑1 (MQ=255)
                                                                                                        gTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCca                                            >  1:271406/1‑150 (MQ=255)
                                                                                                            gCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGATTGGTGACGCCACAACGATAATCcagcag                                        <  2:3614963/150‑1 (MQ=255)
                                                                                                             cccAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGa                                       >  1:1933979/1‑150 (MQ=255)
                                                                                                             cccAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGa                                       <  2:364259/150‑1 (MQ=255)
                                                                                                               cAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATg                                     <  2:2062060/150‑1 (MQ=255)
                                                                                                                aGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGc                                    >  2:4475917/1‑150 (MQ=255)
                                                                                                                  aaTCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATAAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCtg                                  <  2:570686/150‑1 (MQ=255)
                                                                                                                     cGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGc                               >  2:4300209/1‑150 (MQ=255)
                                                                                                                       ccAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcct                             >  2:3895814/1‑150 (MQ=255)
                                                                                                                        cAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcctc                            <  2:3878429/150‑1 (MQ=255)
                                                                                                                         aaCGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcctcc                           >  1:2630142/1‑150 (MQ=255)
                                                                                                                         aaCGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCca                                                         >  1:1624769/1‑120 (MQ=255)
                                                                                                                          aCGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcctcct                          >  2:1597064/1‑150 (MQ=255)
                                                                                                                          aCGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcctcct                          >  2:850259/1‑150 (MQ=255)
                                                                                                                           cGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTg                         <  1:1451295/150‑1 (MQ=255)
                                                                                                                             ccTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCg                       <  1:2337805/150‑1 (MQ=255)
                                                                                                                             ccTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCg                       <  1:36019/150‑1 (MQ=255)
                                                                                                                                 cGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCtt                   >  2:3010520/1‑150 (MQ=255)
                                                                                                                                 cGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCAACGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCtt                   >  2:729571/1‑150 (MQ=255)
                                                                                                                                    cTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCaa                <  2:2631023/150‑1 (MQ=255)
                                                                                                                                    cTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCaa                >  2:2606553/1‑150 (MQ=255)
                                                                                                                                    cTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCaa                <  2:82342/150‑1 (MQ=255)
                                                                                                                                     tGATAATCTACAATTTTTACAATCAGCTTACGCCCGCCAGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAAc               <  1:533217/150‑1 (MQ=255)
                                                                                                                                      gATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACt              <  2:3177923/150‑1 (MQ=255)
                                                                                                                                      gATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACt              >  2:3916745/1‑150 (MQ=255)
                                                                                                                                       aTAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTg             >  1:2667878/1‑150 (MQ=255)
                                                                                                                                        tAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGa            >  2:1813237/1‑150 (MQ=255)
                                                                                                                                         aaTCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAc           >  1:2017564/1‑150 (MQ=255)
                                                                                                                                         aaTCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGATTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAc           >  1:1804050/1‑150 (MQ=255)
                                                                                                                                         aaTCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCCTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAc           >  1:490425/1‑150 (MQ=255)
                                                                                                                                          aTCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAcg          <  1:1085415/150‑1 (MQ=255)
                                                                                                                                          aTCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAcg          >  1:806240/1‑150 (MQ=255)
                                                                                                                                              tCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAcggtcg      <  2:1735208/150‑1 (MQ=255)
                                                                                                                                                aaTTTTTTCAACCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAcggtcggt    >  2:1783894/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAA  >  CP009273/3208401‑3208696

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