Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,738,431 C→A 90.6% P90T (CCA→ACA)  yiaN → 2,3‑diketo‑L‑gulonate TRAP transporter large permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,738,4310CA90.6% 1048.3 / 91.7 337P90T (CCA→ACA) yiaN2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
Reads supporting (aligned to +/- strand):  ref base C (13/19);  new base A (184/121);  total (197/140)
Fisher's exact test for biased strand distribution p-value = 3.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 337 total aligned reads displayed.

AGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGC  >  CP009273/3738285‑3738578
                                                                                                                                                  |                                                                                                                                                   
aGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACa                                                                                                                                                   <  1:2604131/149‑1 (MQ=255)
aGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACa                                                                                                                                                   >  2:3525965/1‑149 (MQ=255)
 gATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAg                                                                                                                                                  >  1:1200976/1‑149 (MQ=255)
 gATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAg                                                                                                                                                  >  2:1784965/1‑149 (MQ=255)
   tCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAAAGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGc                                                                                                                                                >  1:4797930/1‑149 (MQ=255)
    cATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCg                                                                                                                                               <  2:1313143/149‑1 (MQ=255)
     aTGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCgg                                                                                                                                              >  1:2030599/1‑149 (MQ=255)
     aTGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCgg                                                                                                                                              <  1:3746247/149‑1 (MQ=255)
     aTGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCgg                                                                                                                                              >  1:3053726/1‑149 (MQ=255)
       ggCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCgggc                                                                                                                                            >  1:2985618/1‑149 (MQ=255)
        gCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCgggct                                                                                                                                           <  2:3105241/149‑1 (MQ=255)
             aaCGCTGGTGAACGGAGCCGATAGCTTCACCCTGCTGGCGATTCCGTTCTTTGTGCTGTCGGGTGAAAACCTGAATGCGGGCGGCGTGTCAAAACGCAGTGTTGACCTGCCGAAGATACTGGTGGTGCATACAACAAGCgggcttggat                                                                                                                                      >  2:1393361/1‑147 (MQ=255)
                 cTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGt                                                                                                                                  <  2:5124439/149‑1 (MQ=255)
                 cTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGt                                                                                                                                  <  2:4578040/149‑1 (MQ=255)
                  tGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTc                                                                                                                                 <  1:4201543/149‑1 (MQ=255)
                  tGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTc                                                                                                                                 <  1:855785/149‑1 (MQ=255)
                   ggTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCg                                                                                                                                >  2:3855559/1‑149 (MQ=255)
                   ggTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCg                                                                                                                                >  2:795574/1‑149 (MQ=255)
                    gTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCgg                                                                                                                               <  1:2110853/149‑1 (MQ=255)
                      gAACTGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGAGGGTGATATCATGAATGCGTGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCg                                                                                                                             <  1:1838808/149‑1 (MQ=255)
                      gAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTGTTTGTGCTGGCGGGTGAAATAATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCg                                                                                                                             >  1:4506957/1‑149 (MQ=255)
                      gAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCg                                                                                                                             >  2:273836/1‑149 (MQ=255)
                       aaCGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGt                                                                                                                            >  2:2977242/1‑149 (MQ=255)
                        aCGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTg                                                                                                                           >  2:299947/1‑149 (MQ=255)
                         cGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGc                                                                                                                          <  1:4579455/149‑1 (MQ=255)
                           gAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATTAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTg                                                                                                                        >  1:5044197/1‑149 (MQ=255)
                           gAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTg                                                                                                                        >  2:4294628/1‑149 (MQ=255)
                           gAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAAACATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCGGGGCTACGTCGGCGTGCTg                                                                                                                        >  1:3772866/1‑149 (MQ=255)
                             gCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGc                                                                                                                      <  2:5233048/149‑1 (MQ=255)
                             gCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGc                                                                                                                      <  2:4933111/149‑1 (MQ=255)
                                gATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGc                                                                                                                   <  1:2476615/149‑1 (MQ=255)
                                gATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGc                                                                                                                   <  2:5008113/149‑1 (MQ=255)
                                gATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGc                                                                                                                   >  1:4069415/1‑149 (MQ=255)
                                 aTAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCt                                                                                                                  >  1:1721919/1‑149 (MQ=255)
                                 aTAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCt                                                                                                                  >  1:1668838/1‑149 (MQ=255)
                                  tAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGTGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTa                                                                                                                 <  1:3652965/149‑1 (MQ=255)
                                   aGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTAt                                                                                                                >  1:1000848/1‑149 (MQ=255)
                                    gCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATg                                                                                                               <  2:1662579/149‑1 (MQ=255)
                                     cTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGa                                                                                                              <  1:5126270/149‑1 (MQ=255)
                                      ttCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGAt                                                                                                             <  2:1445418/149‑1 (MQ=255)
                                         tcCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTAt                                                                                                          >  1:2008570/1‑149 (MQ=255)
                                           ccTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATgg                                                                                                        <  1:3610429/149‑1 (MQ=255)
                                                tgGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGc                                                                                                   >  1:442478/1‑149 (MQ=255)
                                                tgGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGc                                                                                                   <  1:2803849/149‑1 (MQ=255)
                                                 ggCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGcc                                                                                                  >  1:1581669/1‑149 (MQ=255)
                                                  gCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCt                                                                                                 >  2:2772373/1‑149 (MQ=255)
                                                      ttCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCt                                                                                             <  2:448871/149‑1 (MQ=255)
                                                         cGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGt                                                                                          >  2:66798/1‑149 (MQ=255)
                                                          gTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGtt                                                                                         >  2:1984961/1‑149 (MQ=255)
                                                          gTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGtt                                                                                         >  2:2622115/1‑149 (MQ=255)
                                                           ttCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTc                                                                                        >  1:1887771/1‑149 (MQ=255)
                                                            tCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCt                                                                                       >  2:4110704/1‑149 (MQ=255)
                                                                  tgCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGTTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGtt                                                                                 >  1:3997214/1‑149 (MQ=255)
                                                                  tgCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGtt                                                                                 >  2:4262724/1‑149 (MQ=255)
                                                                     tGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCg                                                                              >  1:4963465/1‑149 (MQ=255)
                                                                      ggCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCgg                                                                             >  1:714195/1‑149 (MQ=255)
                                                                        cGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGAt                                                                           <  1:2850623/149‑1 (MQ=255)
                                                                         gggTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATa                                                                          >  2:1223577/1‑149 (MQ=255)
                                                                         gggTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATa                                                                          >  1:2643426/1‑149 (MQ=255)
                                                                         gggTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATa                                                                          >  2:1038158/1‑149 (MQ=255)
                                                                          ggTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATAc                                                                         >  1:3242675/1‑149 (MQ=255)
                                                                              aaaTCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATAccgcc                                                                     >  2:930049/1‑149 (MQ=255)
                                                                              aaaTCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATAccgcc                                                                     >  2:4034996/1‑149 (MQ=255)
                                                                                aTCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATAccgccgc                                                                   <  1:129240/149‑1 (MQ=255)
                                                                                   aTGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGt                                                                >  1:1444009/1‑149 (MQ=255)
                                                                                    tGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGtt                                                               >  2:1505373/1‑149 (MQ=255)
                                                                                    tGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGtt                                                               >  2:4902717/1‑149 (MQ=255)
                                                                                      aaTGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGc                                                             >  1:2023531/1‑149 (MQ=255)
                                                                                      aaTGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGc                                                             <  1:2310133/149‑1 (MQ=255)
                                                                                         gcgggcggCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGc                                                          >  2:3363550/1‑149 (MQ=255)
                                                                                            ggcggcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGctct                                                       >  2:3732179/1‑149 (MQ=255)
                                                                                             gcggcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTa                                                      <  2:4859415/149‑1 (MQ=255)
                                                                                              cggcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTAc                                                     <  2:3640188/149‑1 (MQ=255)
                                                                                               ggcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACt                                                    <  1:912085/149‑1 (MQ=255)
                                                                                                gcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTCTGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTg                                                   >  2:4154454/1‑149 (MQ=255)
                                                                                                gcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTg                                                   >  2:2018712/1‑149 (MQ=255)
                                                                                                gcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTg                                                   >  1:874367/1‑149 (MQ=255)
                                                                                                        aaaCGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCgat                                           >  2:2607850/1‑149 (MQ=255)
                                                                                                           cGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCgatgat                                        <  1:760032/149‑1 (MQ=255)
                                                                                                            gCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGatgatg                                       >  1:4725993/1‑149 (MQ=255)
                                                                                                            gCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGatgatg                                       >  1:3113510/1‑149 (MQ=255)
                                                                                                               ttgttgACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATgcgc                                    >  2:4371628/1‑149 (MQ=255)
                                                                                                                 gttgACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAg                                  >  1:5454/1‑149 (MQ=255)
                                                                                                                     aCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGcc                              <  2:4575534/149‑1 (MQ=255)
                                                                                                                       cTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCaa                            >  2:3835689/1‑149 (MQ=255)
                                                                                                                         gCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACt                          <  2:2515838/149‑1 (MQ=255)
                                                                                                                          ccGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTa                         <  2:4473656/149‑1 (MQ=255)
                                                                                                                            gATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATc                       <  2:624159/149‑1 (MQ=255)
                                                                                                                              tGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCg                     >  2:1174760/1‑149 (MQ=255)
                                                                                                                                 aaCTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTc                  >  2:2115460/1‑149 (MQ=255)
                                                                                                                                 aaCTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTc                  >  2:1939399/1‑149 (MQ=255)
                                                                                                                                 aaCTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTc                  >  1:4816986/1‑149 (MQ=255)
                                                                                                                                   cTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCaa                <  1:3076714/149‑1 (MQ=255)
                                                                                                                                      gtggGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCTATGATGCGCAGCGCCAACTATCCGGTCAACCg             <  1:347287/149‑1 (MQ=255)
                                                                                                                                      gtggGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCg             <  2:2339785/149‑1 (MQ=255)
                                                                                                                                       tggGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCgg            >  2:4232394/1‑149 (MQ=255)
                                                                                                                                        ggggCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCggg           >  2:3747246/1‑149 (MQ=255)
                                                                                                                                          ggCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGggcg         <  1:1091859/149‑1 (MQ=255)
                                                                                                                                           gCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGggcgg        >  1:1137782/1‑149 (MQ=255)
                                                                                                                                             aTAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGgcggcg      >  2:4775556/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGC  >  CP009273/3738285‑3738578

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: