Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,618,536 C→A 90.9% R267R (CGG→AGG)  yhhI → putative transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,618,5360CA90.9% 401.7 / 40.9 184R267R (CGG→AGG) yhhIputative transposase
Reads supporting (aligned to +/- strand):  ref base C (10/7);  new base A (82/85);  total (92/92)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 184 total aligned reads displayed.

AGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTG  >  CP009273/3618391‑3618682
                                                                                                                                                 |                                                                                                                                                  
  agCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTc                                                                                                                                                >  2:2792499/1‑148 (MQ=17)
   gCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTcc                                                                                                                                               <  2:367322/148‑1 (MQ=17)
    cATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCa                                                                                                                                              >  2:1353210/1‑148 (MQ=32)
    cATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCa                                                                                                                                              <  2:1159684/148‑1 (MQ=17)
     aTGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCat                                                                                                                                             >  1:3211279/1‑148 (MQ=17)
      tGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTACata                                                                                                                                            <  1:967730/148‑1 (MQ=17)
       taCAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCataa                                                                                                                                           >  2:2933441/2‑148 (MQ=17)
        aCAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCataat                                                                                                                                          >  2:2047558/1‑148 (MQ=17)
        aCAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCataat                                                                                                                                          >  2:381050/1‑148 (MQ=17)
           gTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGc                                                                                                                                       <  1:92304/148‑1 (MQ=17)
              aCGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGa                                                                                                                                    <  1:2406359/148‑1 (MQ=17)
               cGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGaa                                                                                                                                   <  1:2686210/148‑1 (MQ=17)
                    tCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACaaaa                                                                                                                              >  2:1421078/1‑148 (MQ=14)
                    tCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACaaaa                                                                                                                              >  2:364458/1‑148 (MQ=14)
                     cAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAaaag                                                                                                                             >  2:2728342/1‑148 (MQ=11)
                     cAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAaaag                                                                                                                             <  1:2912877/148‑1 (MQ=11)
                     cAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAaaag                                                                                                                             >  2:928366/1‑148 (MQ=11)
                             agagTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGcc                                                                                                                     >  1:2002971/1‑148 (MQ=11)
                             agagTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGcc                                                                                                                     <  2:1762927/148‑1 (MQ=11)
                               agTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAg                                                                                                                   <  1:210938/148‑1 (MQ=17)
                                     ggCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGa                                                                                                             <  2:237760/148‑1 (MQ=21)
                                     ggCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGa                                                                                                             <  1:1215164/148‑1 (MQ=21)
                                       cAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACg                                                                                                           >  2:1400821/1‑148 (MQ=21)
                                         gagaAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGt                                                                                                         >  1:1630161/1‑148 (MQ=21)
                                           gaagaaATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGACATGACGGTCa                                                                                                       >  2:1195132/1‑148 (MQ=21)
                                           gaagaaATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCa                                                                                                       >  2:792883/1‑148 (MQ=21)
                                            aagaaATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGTAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGACAGAAATGACGGTCAg                                                                                                      >  2:101581/1‑148 (MQ=255)
                                             agtaaTGCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCata                                                                                                     <  1:2232145/145‑3 (MQ=19)
                                              gaaATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGat                                                                                                    <  2:3137488/148‑1 (MQ=21)
                                                 aTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtatt                                                                                                 <  1:89045/148‑1 (MQ=21)
                                                 aTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtatt                                                                                                 <  2:603886/148‑1 (MQ=21)
                                                 aTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtatt                                                                                                 <  2:1501505/148‑1 (MQ=21)
                                                   ccGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtattat                                                                                               <  2:1044870/148‑1 (MQ=21)
                                                    cGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATtata                                                                                              >  2:2973329/1‑148 (MQ=21)
                                                     gTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTatat                                                                                             <  2:463893/148‑1 (MQ=21)
                                                      tCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATc                                                                                            >  2:140146/1‑148 (MQ=21)
                                                       cTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCa                                                                                           >  2:557617/1‑148 (MQ=21)
                                                       cTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCa                                                                                           >  2:1902695/1‑148 (MQ=21)
                                                        ttCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAg                                                                                          >  2:2139910/1‑148 (MQ=21)
                                                          cATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGtt                                                                                        <  1:1157575/148‑1 (MQ=21)
                                                           atatTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTc                                                                                       <  2:1254174/148‑1 (MQ=34)
                                                           atatTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTc                                                                                       <  2:1911082/148‑1 (MQ=34)
                                                             atTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTctg                                                                                     <  1:2888374/148‑1 (MQ=21)
                                                             atTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTctg                                                                                     <  1:3164692/148‑1 (MQ=21)
                                                               tGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGATCCAGAAATGACGGTCAGATATTATATCAGTTctgct                                                                                   >  2:940157/1‑148 (MQ=21)
                                                               tGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTctgct                                                                                   >  1:1936824/1‑148 (MQ=21)
                                                               tGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTctgct                                                                                   >  1:1494332/1‑148 (MQ=21)
                                                                gTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTctgctg                                                                                  <  2:742200/148‑1 (MQ=21)
                                                                gTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTctgctg                                                                                  <  1:442337/148‑1 (MQ=21)
                                                                 tttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGGCGGTCAGATATTATATCAGTTCTGCTGa                                                                                 <  1:3022501/148‑1 (MQ=21)
                                                                 tttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGa                                                                                 >  1:2560045/1‑148 (MQ=21)
                                                                  ttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGAt                                                                                <  2:1426186/148‑1 (MQ=21)
                                                                  ttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGAt                                                                                <  1:1413283/148‑1 (MQ=21)
                                                                      gATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTaa                                                                            <  1:1552491/148‑1 (MQ=21)
                                                                         gTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCg                                                                         >  1:1664852/1‑148 (MQ=34)
                                                                          tCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGc                                                                        >  2:2578126/1‑148 (MQ=21)
                                                                          tCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGc                                                                        >  1:1443745/1‑148 (MQ=21)
                                                                            ccTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTg                                                                      <  2:1564380/148‑1 (MQ=25)
                                                                                   aaCTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAAtt                                                               >  1:1244168/1‑148 (MQ=255)
                                                                                   aaCTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAAtt                                                               <  1:1872448/148‑1 (MQ=255)
                                                                                    aCTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTc                                                              >  2:572054/1‑148 (MQ=255)
                                                                                      ttATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGc                                                            <  2:1863558/148‑1 (MQ=255)
                                                                                      ttATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGc                                                            >  1:2817713/1‑148 (MQ=255)
                                                                                       tATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGcc                                                           <  1:1366768/148‑1 (MQ=255)
                                                                                       tATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGcc                                                           >  1:1568247/1‑148 (MQ=255)
                                                                                        aTTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCca                                                          >  1:2398688/1‑148 (MQ=255)
                                                                                        aTTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCca                                                          >  1:470880/1‑148 (MQ=255)
                                                                                          tGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCcaca                                                        >  1:664116/1‑148 (MQ=255)
                                                                                             tttCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCg                                                     >  1:474378/1‑148 (MQ=255)
                                                                                              ttCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGa                                                    >  1:1959413/1‑148 (MQ=255)
                                                                                               tCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGAt                                                   >  2:1901006/1‑148 (MQ=255)
                                                                                                cACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATc                                                  <  1:2264309/148‑1 (MQ=255)
                                                                                                  cGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCg                                                >  1:2741641/1‑148 (MQ=255)
                                                                                                    tttGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGaa                                              <  2:646273/148‑1 (MQ=255)
                                                                                                    tttGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGaa                                              <  1:2047558/148‑1 (MQ=255)
                                                                                                        aaTGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCa                                          <  1:3255733/148‑1 (MQ=255)
                                                                                                             aaaGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGc                                     >  1:172197/1‑148 (MQ=255)
                                                                                                             aaaGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGc                                     >  1:1087624/1‑148 (MQ=255)
                                                                                                                 ggCTGAAGAACTTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGt                                 <  2:3206952/148‑1 (MQ=255)
                                                                                                                 ggCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGt                                 <  2:2881163/148‑1 (MQ=255)
                                                                                                                 ggCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGt                                 >  1:1093343/1‑148 (MQ=255)
                                                                                                                 ggCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGt                                 <  2:1828703/148‑1 (MQ=255)
                                                                                                                   cTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTgg                               <  1:2466507/148‑1 (MQ=255)
                                                                                                                    tGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGga                              >  2:2922221/1‑148 (MQ=255)
                                                                                                                    tGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGga                              >  2:1273929/1‑148 (MQ=255)
                                                                                                                      aagaaATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGgaga                            >  1:2317967/1‑148 (MQ=255)
                                                                                                                      aagaaATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGgaga                            >  2:1649736/1‑148 (MQ=255)
                                                                                                                       agaaATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGaa                           >  1:523643/1‑148 (MQ=255)
                                                                                                                        gaaATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAAt                          <  1:721938/148‑1 (MQ=255)
                                                                                                                         aaaTTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATa                         <  2:1990708/148‑1 (MQ=255)
                                                                                                                           aTTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAg                       <  1:578201/148‑1 (MQ=255)
                                                                                                                           aTTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGGGGAGCATAAg                       >  1:1686242/1‑148 (MQ=255)
                                                                                                                            ttATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGc                      <  1:339450/148‑1 (MQ=255)
                                                                                                                                   tGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTg               <  1:1998975/148‑1 (MQ=255)
                                                                                                                                    ggCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTgg              <  1:3193664/148‑1 (MQ=255)
                                                                                                                                      cAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCg            >  1:981988/1‑148 (MQ=255)
                                                                                                                                         tctcCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCt         >  1:276819/1‑148 (MQ=255)
                                                                                                                                         tctcCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCt         >  1:2541380/1‑148 (MQ=255)
                                                                                                                                           tcCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTgg       <  2:3060716/148‑1 (MQ=255)
                                                                                                                                           tcCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTgg       <  2:909101/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTG  >  CP009273/3618391‑3618682

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: