Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,208,550 T→A 42.5% 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,5500TA42.5% 142.3 / inf 240K95* (AAA→TAA) mugG/U mismatch‑specific DNA glycosylase; xanthine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base T (74/64);  new base A (58/43);  total (133/107)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 245 total aligned reads displayed.

TAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACC  >  CP009273/3208401‑3208699
                                                                                                                                                     |                                                                                                                                                     
       cTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATc                                                                                                                                                <  1:3740696/150‑1 (MQ=255)
       cTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATc                                                                                                                                                >  1:766816/1‑150 (MQ=255)
       cTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGAGAATCTTCAATTTTTTCAATc                                                                                                                                                >  1:468104/1‑150 (MQ=255)
        tGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCa                                                                                                                                               <  2:2722148/150‑1 (MQ=255)
        tGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTCTTCAATCa                                                                                                                                               >  2:16010/1‑150 (MQ=255)
        tGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCa                                                                                                                                               >  2:16141/1‑150 (MQ=255)
         gggATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAg                                                                                                                                              <  1:1835613/150‑1 (MQ=255)
            aTTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCtt                                                                                                                                           <  2:520359/150‑1 (MQ=255)
            aTTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCtt                                                                                                                                           <  2:1996609/150‑1 (MQ=255)
             tttGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTa                                                                                                                                          <  2:1880552/150‑1 (MQ=255)
                   agcaCCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTAcgcccg                                                                                                                                    <  1:92461/150‑1 (MQ=255)
                   agcaCCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTAcgcccg                                                                                                                                    <  1:3789450/150‑1 (MQ=255)
                   agcaCCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTAcgcccg                                                                                                                                    >  1:3623454/1‑150 (MQ=255)
                       cccAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTg                                                                                                                                <  2:802312/150‑1 (MQ=255)
                          aaaTCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCgt                                                                                                                             >  1:647832/1‑150 (MQ=255)
                          aaaTCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCgt                                                                                                                             >  1:3016391/1‑150 (MQ=255)
                           aaTCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCgtg                                                                                                                            >  1:1112602/1‑150 (MQ=255)
                           aaTCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCgtg                                                                                                                            >  2:3444266/1‑150 (MQ=255)
                              cTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAg                                                                                                                         >  2:3602709/1‑150 (MQ=255)
                               tGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGc                                                                                                                        >  2:1360095/1‑150 (MQ=255)
                                   tCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCt                                                                                                                    <  2:1446354/150‑1 (MQ=255)
                                   tCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCt                                                                                                                    >  2:187567/1‑150 (MQ=255)
                                    tgAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTg                                                                                                                   >  1:1576170/2‑150 (MQ=255)
                                     gAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGc                                                                                                                  >  1:533286/1‑150 (MQ=255)
                                      agccAATGGTGAGCGTTTGTTTCCCCCACTGTGCGCCGCGCTGGCTGAATCCCTGTTCATAGGCTTGTTTGCCCAGAATCGCCCACGCCTGCGGCTGAGAAGCTTCAATTTTTTCTATCAGCTTTCGCCCGCCTGCGTGTAGCTCCTGCt                                                                                                                 >  1:4246873/3‑150 (MQ=255)
                                      aaCCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCt                                                                                                                 >  1:947798/1‑150 (MQ=255)
                                       aCCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCtt                                                                                                                <  1:3217592/150‑1 (MQ=255)
                                       aCCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCtt                                                                                                                >  1:3190635/1‑150 (MQ=255)
                                        ccAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCttt                                                                                                               <  2:1806443/150‑1 (MQ=255)
                                          aaTGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGa                                                                                                             >  1:4099117/1‑150 (MQ=255)
                                           aTGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGaa                                                                                                            <  2:3853801/150‑1 (MQ=255)
                                           aTGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGaa                                                                                                            >  1:1211723/1‑150 (MQ=255)
                                           aTGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGaa                                                                                                            <  1:2412280/150‑1 (MQ=255)
                                            tGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGaaa                                                                                                           <  1:3870921/150‑1 (MQ=255)
                                             ggTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAAc                                                                                                          >  1:2261663/1‑150 (MQ=255)
                                              gTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACt                                                                                                         <  1:1723186/150‑1 (MQ=255)
                                              gTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACt                                                                                                         <  2:3306396/150‑1 (MQ=255)
                                              gTGAGCGTTTGTTTCACCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACt                                                                                                         <  1:1960886/150‑1 (MQ=255)
                                                 aGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCa                                                                                                      <  1:2739206/150‑1 (MQ=255)
                                                 aGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCa                                                                                                      <  2:253623/150‑1 (MQ=255)
                                                   cGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCAtt                                                                                                    >  1:3485270/1‑150 (MQ=255)
                                                   cGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCAtt                                                                                                    <  2:3563279/150‑1 (MQ=255)
                                                   cGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCAtt                                                                                                    <  2:3182719/150‑1 (MQ=255)
                                                    gtttgtttCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTg                                                                                                   >  1:3348731/1‑150 (MQ=255)
                                                      ttgtttCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGc                                                                                                 <  2:1253701/150‑1 (MQ=255)
                                                       tgtttCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCt                                                                                                >  1:4370262/1‑150 (MQ=255)
                                                         tttCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTg                                                                                              >  2:376878/1‑150 (MQ=255)
                                                          ttCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGc                                                                                             >  1:2367400/1‑150 (MQ=255)
                                                              cccACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCg                                                                                         <  2:280112/150‑1 (MQ=255)
                                                                cACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGtt                                                                                       <  1:4027545/150‑1 (MQ=255)
                                                                  cTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTgg                                                                                     >  2:486642/1‑150 (MQ=255)
                                                                  cTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTgg                                                                                     >  2:3466950/1‑150 (MQ=255)
                                                                  cTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTgg                                                                                     >  2:1143689/1‑150 (MQ=255)
                                                                  cTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTgg                                                                                     <  2:2962312/150‑1 (MQ=255)
                                                                   tgtgCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGa                                                                                    >  1:1614533/1‑150 (MQ=255)
                                                                   tgtgCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGa                                                                                    >  1:4227514/1‑150 (MQ=255)
                                                                      gCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACgg                                                                                 >  2:247388/1‑150 (MQ=255)
                                                                          cgcgCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTa                                                                             <  2:3080607/150‑1 (MQ=255)
                                                                            cgcTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTAcc                                                                           >  1:3920762/1‑150 (MQ=255)
                                                                              ctggctgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAg                                                                         >  2:965774/1‑150 (MQ=255)
                                                                                   tgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTgg                                                                    >  1:1025501/1‑150 (MQ=255)
                                                                                    gAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGt                                                                   >  2:2993272/1‑150 (MQ=255)
                                                                                    gAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGt                                                                   <  1:1198471/150‑1 (MQ=255)
                                                                                      aTCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGa                                                                 <  2:4062367/150‑1 (MQ=255)
                                                                                      aTCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGa                                                                 <  2:4114962/150‑1 (MQ=255)
                                                                                       tCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGAc                                                                <  1:4382928/150‑1 (MQ=255)
                                                                                          cTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGcc                                                             <  1:2625522/150‑1 (MQ=255)
                                                                                           tGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCca                                                            >  1:2881308/1‑150 (MQ=255)
                                                                                             ttCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCcaca                                                          >  2:1687301/1‑150 (MQ=255)
                                                                                                atatGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACg                                                       <  2:678939/150‑1 (MQ=255)
                                                                                                atatGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACg                                                       >  1:447548/1‑150 (MQ=255)
                                                                                                atatGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACg                                                       >  1:2119427/1‑150 (MQ=255)
                                                                                                    gCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATaa                                                   <  2:1759639/150‑1 (MQ=255)
                                                                                                          ttGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCcagc                                             >  2:3989448/1‑150 (MQ=255)
                                                                                                          ttGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCcagc                                             >  1:2165522/1‑150 (MQ=255)
                                                                                                             cccAGACTCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGa                                          <  1:3531308/150‑1 (MQ=255)
                                                                                                                aGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGc                                       >  1:4012471/1‑150 (MQ=255)
                                                                                                                  aaTCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCtg                                     <  1:1340430/150‑1 (MQ=255)
                                                                                                                  aaTCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCtg                                     <  1:1162502/150‑1 (MQ=255)
                                                                                                                   aTCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTTGCTCCTGATTTGAAACTTTTTTGGCTTGCACAGTTGGACGGTATTCCAGTTTGGTGACGCCTCAACGATAATCCAGCAGATGCtgt                                    <  2:600171/150‑1 (MQ=255)
                                                                                                                    tCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCtgtg                                   <  1:4190954/150‑1 (MQ=255)
                                                                                                                    tCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCtgtg                                   >  1:4071578/1‑150 (MQ=255)
                                                                                                                       ccAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcct                                >  1:3786575/1‑150 (MQ=255)
                                                                                                                        cAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCAGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCAGTGcctc                               >  2:3345532/1‑150 (MQ=255)
                                                                                                                        cAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcctc                               >  2:2420931/1‑150 (MQ=255)
                                                                                                                         aaCGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcctcc                              <  1:1051911/150‑1 (MQ=255)
                                                                                                                         aaCGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcctcc                              <  1:1051459/150‑1 (MQ=255)
                                                                                                                           cGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTg                            <  1:2520579/150‑1 (MQ=255)
                                                                                                                           cGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTg                            <  1:4074137/150‑1 (MQ=255)
                                                                                                                            gCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGc                           >  1:3181174/1‑150 (MQ=255)
                                                                                                                             ccTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCg                          >  1:450587/1‑150 (MQ=255)
                                                                                                                               tGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGTCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGc                        >  2:4255745/1‑150 (MQ=255)
                                                                                                                                 cGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCtt                      >  1:4204631/1‑150 (MQ=255)
                                                                                                                                 cGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCGCGGCTTc                     <  2:1016191/150‑1 (MQ=255)
                                                                                                                                      gATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACt                 <  2:2021391/150‑1 (MQ=255)
                                                                                                                                      gATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACt                 <  2:1211723/150‑1 (MQ=255)
                                                                                                                                          aTCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAcg             <  2:2192476/150‑1 (MQ=255)
                                                                                                                                           tCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAcgg            >  2:3227393/1‑150 (MQ=255)
                                                                                                                                                aaTTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAcggtcggt       <  2:4490933/150‑1 (MQ=255)
                                                                                                                                                     tttCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAAcc  >  1:2918679/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACC  >  CP009273/3208401‑3208699

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: