Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,208,581 C→A 49.4% E84D (GAG→GAT 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,5810CA49.4% 126.3 / inf 273E84D (GAG→GATmugG/U mismatch‑specific DNA glycosylase; xanthine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base C (81/57);  new base A (76/59);  total (157/116)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 276 total aligned reads displayed.

GCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCTTCCAGAAGCGATT  >  CP009273/3208433‑3208730
                                                                                                                                                    |                                                                                                                                                     
    cGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTg                                                                                                                                                  <  1:282713/150‑1 (MQ=255)
    cGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTg                                                                                                                                                  <  1:2148689/150‑1 (MQ=255)
     gAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGc                                                                                                                                                 >  2:4244922/1‑150 (MQ=255)
       aCCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCtt                                                                                                                                               >  1:1964592/1‑150 (MQ=255)
         cAATGGTGAGCGTTTGTTTCCCCCACTGTGTACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTg                                                                                                                                             >  2:3473500/1‑150 (MQ=255)
         cAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTg                                                                                                                                             >  2:2989677/1‑150 (MQ=255)
           aTGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGaa                                                                                                                                           >  2:4188204/1‑150 (MQ=255)
           aTGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGaa                                                                                                                                           >  2:49376/1‑150 (MQ=255)
             ggTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAAc                                                                                                                                         <  1:2926123/150‑1 (MQ=255)
               tGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACtt                                                                                                                                       >  1:4023161/1‑150 (MQ=255)
               tGAGCGTTTGTTTCCCCCACGGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACtt                                                                                                                                       <  2:229146/150‑1 (MQ=255)
                gAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTc                                                                                                                                      <  2:1766089/150‑1 (MQ=255)
                   cGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCAtt                                                                                                                                   <  2:4008297/150‑1 (MQ=255)
                      ttgtttCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGc                                                                                                                                <  1:1507061/150‑1 (MQ=255)
                       tgtttCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCt                                                                                                                               <  1:957153/150‑1 (MQ=255)
                         tttCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTg                                                                                                                             >  2:1139186/1‑150 (MQ=255)
                             ccccACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCAcc                                                                                                                         <  2:1741629/150‑1 (MQ=255)
                              cccACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCg                                                                                                                        >  1:1042240/1‑150 (MQ=255)
                                cccTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGtt                                                                                                                      <  2:60372/148‑1 (MQ=255)
                                  cTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTgg                                                                                                                    <  1:3070218/150‑1 (MQ=255)
                                  cTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTgg                                                                                                                    <  2:2964153/150‑1 (MQ=255)
                                     tgCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACg                                                                                                                 >  2:3339278/1‑150 (MQ=255)
                                         ccGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCt                                                                                                             >  2:2257402/1‑150 (MQ=255)
                                          cgcgCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTa                                                                                                            >  2:3942936/1‑150 (MQ=255)
                                             gctggctgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCa                                                                                                         <  2:952405/150‑1 (MQ=255)
                                              ctggctgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAg                                                                                                        >  1:1383565/1‑150 (MQ=255)
                                              ctggctgAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAg                                                                                                        >  1:371096/1‑150 (MQ=255)
                                                    gAATCCCTGTTCACATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTggc                                                                                                  <  1:518879/150‑2 (MQ=255)
                                                       tCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGAc                                                                                               <  2:71618/150‑1 (MQ=255)
                                                       tCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGAc                                                                                               <  2:2661294/150‑1 (MQ=255)
                                                          cTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGcc                                                                                            >  1:1100293/1‑150 (MQ=255)
                                                           tGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCca                                                                                           <  1:1823290/150‑1 (MQ=255)
                                                            gTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCcac                                                                                          <  1:3765064/150‑1 (MQ=255)
                                                            gTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCcac                                                                                          >  1:3988543/1‑150 (MQ=255)
                                                             ttCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCcaca                                                                                         <  1:986031/150‑1 (MQ=255)
                                                              tCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACaa                                                                                        >  1:1775570/1‑150 (MQ=255)
                                                                     cTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAAt                                                                                 <  2:2768842/150‑1 (MQ=255)
                                                                      ttGTTTGACCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCATTGAAACTTCATTGGCTTGCACCGTTGGACTGTCTACCAGTTTGGTGACGCCACAACGATAATc                                                                                <  1:2257769/150‑1 (MQ=255)
                                                                       tGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATcc                                                                               <  1:630780/150‑1 (MQ=255)
                                                                       tGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATcc                                                                               <  1:3445856/150‑1 (MQ=255)
                                                                         tttGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGACCCGTTGGACGGTCTACCAGTGTGGTGACGCCACAACGATAATCcag                                                                             >  2:1550007/1‑150 (MQ=255)
                                                                             cccAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGa                                                                         >  1:1245866/1‑150 (MQ=255)
                                                                                aGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGc                                                                      <  2:2671395/150‑1 (MQ=255)
                                                                                 gAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCt                                                                     >  2:132298/1‑150 (MQ=255)
                                                                                       ccAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcct                                                               >  1:1638382/1‑150 (MQ=255)
                                                                                        cAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcctc                                                              <  1:2538197/150‑1 (MQ=255)
                                                                                         aaCGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGcctcc                                                             >  1:2984303/1‑150 (MQ=255)
                                                                                           cGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTg                                                           <  2:1103563/150‑1 (MQ=255)
                                                                                           cGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTg                                                           <  2:3465885/150‑1 (MQ=255)
                                                                                           cGCCTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTg                                                           >  1:3808724/1‑150 (MQ=255)
                                                                                             ccTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCg                                                         >  2:4500415/1‑150 (MQ=255)
                                                                                             ccTGCGGCTGATAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCg                                                         <  2:244406/150‑1 (MQ=255)
                                                                                                gCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCt                                                      <  1:312737/150‑1 (MQ=255)
                                                                                                 cGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCtt                                                     >  1:1691322/1‑150 (MQ=255)
                                                                                                  ggCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTc                                                    <  2:3650345/150‑1 (MQ=255)
                                                                                                   gCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCa                                                   <  1:4188204/150‑1 (MQ=255)
                                                                                                    cTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCaa                                                  <  1:3300444/150‑1 (MQ=255)
                                                                                                     tGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAAc                                                 >  1:3697442/1‑150 (MQ=255)
                                                                                                      gATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACt                                                <  1:865037/150‑1 (MQ=255)
                                                                                                       aTAATCTTCAATTTATTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTg                                               >  1:2154133/1‑150 (MQ=255)
                                                                                                        tAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGa                                              >  2:2375060/1‑150 (MQ=255)
                                                                                                        tAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGa                                              >  2:1663636/1‑150 (MQ=255)
                                                                                                        tAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGa                                              >  2:3789992/1‑150 (MQ=255)
                                                                                                         aaTCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAc                                             <  1:2458099/150‑1 (MQ=255)
                                                                                                         aaTCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAc                                             <  1:148332/150‑1 (MQ=255)
                                                                                                         aaTCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCAGTGGATTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCATCTGAc                                             >  2:4042506/1‑150 (MQ=255)
                                                                                                              tCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGAcggtcg                                        <  1:4560205/150‑1 (MQ=255)
                                                                                                                 aTTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTa                                     >  1:775279/1‑150 (MQ=255)
                                                                                                                  ttttttCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTaa                                    <  2:4557833/150‑1 (MQ=255)
                                                                                                                  ttttttCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTaa                                    >  1:2397223/1‑150 (MQ=255)
                                                                                                                    ttttCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAAc                                  >  1:3143905/1‑150 (MQ=255)
                                                                                                                     tttCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAAcc                                 >  2:356979/1‑150 (MQ=255)
                                                                                                                     tttCAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAAcc                                 <  2:3877617/150‑1 (MQ=255)
                                                                                                                     attCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAAcc                                 >  1:1908382/2‑150 (MQ=255)
                                                                                                                        cAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCATCTGACGGTCGGTAAACCCgg                              <  2:1995426/150‑1 (MQ=255)
                                                                                                                        cAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCgg                              >  1:3869528/1‑150 (MQ=255)
                                                                                                                        cAATCAGCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCgg                              >  1:3869547/1‑150 (MQ=255)
                                                                                                                         aaTCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGc                             >  2:3230667/1‑150 (MQ=255)
                                                                                                                           tCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCt                           >  1:3607657/1‑150 (MQ=255)
                                                                                                                            cAGCTTACTCCCGCCTGCGTGTAGCTCCTTCTTTGAAACTTCATTGGCTTGCAACGTTGTACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTg                          >  2:345488/1‑150 (MQ=255)
                                                                                                                            cAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTg                          >  2:4424450/1‑150 (MQ=255)
                                                                                                                            cAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTg                          >  2:4571085/1‑150 (MQ=255)
                                                                                                                            cAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTg                          >  2:1300746/1‑150 (MQ=255)
                                                                                                                              gCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTCAACCCGGCCTGat                        >  2:4273580/1‑150 (MQ=255)
                                                                                                                              gCTTACGCCCGCCTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGat                        <  1:2788438/150‑1 (MQ=255)
                                                                                                                                ttACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGatat                      >  1:3910489/1‑150 (MQ=255)
                                                                                                                                ttACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGatat                      <  2:3907399/150‑1 (MQ=255)
                                                                                                                                  aCGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGatatat                    >  2:1058299/1‑150 (MQ=255)
                                                                                                                                      ccgccTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCAcc                >  2:533419/1‑150 (MQ=255)
                                                                                                                                      ccgccTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCAcc                >  1:2358133/1‑150 (MQ=255)
                                                                                                                                      ccgccATCGTGAAGCACCTGCTTTGAAACTTCATTGGCTTGCACCGATGGACGGTATACCATTTTGGTGACGCCACAACGATAATCAAGCAGATGCTGTTCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCAcc                >  2:3422212/1‑150 (MQ=255)
                                                                                                                                       cgccTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCt               >  1:2103371/1‑150 (MQ=255)
                                                                                                                                       cgccTGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCt               >  1:264067/1‑150 (MQ=255)
                                                                                                                                         ccTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCTTc             >  2:127985/1‑150 (MQ=255)
                                                                                                                                          cTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCTTcc            >  1:3507742/1‑150 (MQ=255)
                                                                                                                                          cTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCTTcc            >  1:3564401/1‑150 (MQ=255)
                                                                                                                                           tGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCTTCCa           >  2:728206/1‑150 (MQ=255)
                                                                                                                                           tGCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCTTCCa           >  2:2408659/1‑150 (MQ=255)
                                                                                                                                            gCGTGTAGATCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCTCCTTCCAg          <  2:3808373/150‑1 (MQ=255)
                                                                                                                                                    cTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCTTCCAGAAGCGAtt  >  2:2824840/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCACCGCGCTGGCTGAATCCCTGTTCATATGCTTGTTTGCCCAGAATCGCCAACGCCTGCGGCTGATAATCTTCAATTTTTTCAATCAGCTTACGCCCGCCTGCGTGTAGCTCCTGCTTTGAAACTTCATTGGCTTGCACCGTTGGACGGTCTACCAGTTTGGTGACGCCACAACGATAATCCAGCAGATGCTGTGCCTCCTGCGGCTTCAACTGACGGTCGGTAAACCCGGCCTGATATATCACCTTCCAGAAGCGATT  >  CP009273/3208433‑3208730

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: