Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,072,314 C→A 10.8% E126* (GAA→TAA)  cmtB ← putative mannitol‑specific enzyme IIA component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,072,3140CA10.8% 892.9 / 54.0 256E126* (GAA→TAA) cmtBputative mannitol‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base C (139/89);  new base A (13/15);  total (152/104)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 271 total aligned reads displayed.

AGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTT  >  CP009273/3072165‑3072463
                                                                                                                                                     |                                                                                                                                                     
aGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttc                                                                                                                                                       >  2:1956899/1‑150 (MQ=255)
aGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttc                                                                                                                                                       >  2:3817125/1‑150 (MQ=255)
    aaaaaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttcttct                                                                                                                                                   >  2:1233428/1‑150 (MQ=255)
     aaaaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttcttctt                                                                                                                                                  <  1:1803115/150‑1 (MQ=255)
     aaaaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttcttctt                                                                                                                                                  <  1:1803339/150‑1 (MQ=255)
      aaaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttcttcttc                                                                                                                                                 <  2:4068980/150‑1 (MQ=255)
       aaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCg                                                                                                                                                <  1:2487747/150‑1 (MQ=255)
          cccccAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGcac                                                                                                                                             <  1:1274827/150‑1 (MQ=255)
           ccccAAAAGCCTGGTCATTTGCAAGATCAGACTTGTTTTCAATAATATTATCATGTAGAGGAGAGAGTAATTATCCGCGGCTGATATTGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGcaca                                                                                                                                            <  2:3251602/150‑1 (MQ=255)
            cccAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAg                                                                                                                                           >  2:454814/1‑150 (MQ=255)
             ccAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGt                                                                                                                                          >  2:1841677/1‑150 (MQ=255)
              cAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTa                                                                                                                                         <  2:384361/150‑1 (MQ=255)
              cAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTa                                                                                                                                         <  2:1676589/150‑1 (MQ=255)
                  aGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTc                                                                                                                                     <  2:397710/150‑1 (MQ=255)
                    ccTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCAc                                                                                                                                   >  1:3409413/1‑150 (MQ=255)
                     cTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACt                                                                                                                                  >  1:4353147/1‑150 (MQ=255)
                      tGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACtt                                                                                                                                 <  1:610727/150‑1 (MQ=255)
                       ggACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTa                                                                                                                                <  2:2741775/150‑1 (MQ=255)
                          ccTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACg                                                                                                                             <  1:391856/150‑1 (MQ=255)
                          ccTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACg                                                                                                                             <  1:391879/150‑1 (MQ=255)
                           cTTTTCACGAGCAGACTTGATTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGTCTGATAATGTCCTCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTACCGc                                                                                                                            <  1:4252985/150‑1 (MQ=255)
                           cTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGc                                                                                                                            <  1:2573847/150‑1 (MQ=255)
                              tGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTg                                                                                                                         <  2:843161/150‑1 (MQ=255)
                              tGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTg                                                                                                                         <  2:3123906/150‑1 (MQ=255)
                              tGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTg                                                                                                                         >  2:3307789/1‑150 (MQ=255)
                                cACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGaa                                                                                                                       <  1:1354767/150‑1 (MQ=255)
                                cACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGaa                                                                                                                       >  2:433407/1‑150 (MQ=255)
                                  cGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATa                                                                                                                     >  2:1078350/1‑150 (MQ=255)
                                    aGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAgc                                                                                                                   >  1:3474888/1‑150 (MQ=255)
                                    aGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAgc                                                                                                                   >  1:3475251/1‑150 (MQ=255)
                                     gCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAgcg                                                                                                                  <  1:589531/150‑1 (MQ=255)
                                         aCTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCaa                                                                                                              >  1:177150/1‑150 (MQ=255)
                                           ttGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAAtg                                                                                                            <  1:464450/150‑1 (MQ=255)
                                              ttttCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGc                                                                                                         >  1:3239336/1‑150 (MQ=255)
                                                ttCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGa                                                                                                       <  1:2721844/150‑1 (MQ=255)
                                                    ataataTTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAatcg                                                                                                   <  1:1906343/150‑1 (MQ=255)
                                                    ataataTTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAatcg                                                                                                   >  2:3649214/1‑150 (MQ=255)
                                                      aataTTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAatcggc                                                                                                 >  2:449824/1‑150 (MQ=255)
                                                       atatTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAatcggca                                                                                                >  1:3912990/1‑150 (MQ=255)
                                                         attatCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAatcggcatc                                                                                              <  1:977978/150‑1 (MQ=255)
                                                           tatCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAatcggcatcgg                                                                                            <  2:4147918/150‑1 (MQ=255)
                                                            atCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAatcggcatcggc                                                                                           >  1:113625/1‑150 (MQ=255)
                                                            atCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGAGCGAatcggcatcggc                                                                                           >  1:3738292/1‑150 (MQ=255)
                                                              ccTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGgcag                                                                                         >  2:3975343/1‑150 (MQ=255)
                                                               cTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGgcagc                                                                                        >  1:2843726/1‑150 (MQ=255)
                                                                tGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGgcagca                                                                                       <  1:1395490/150‑1 (MQ=255)
                                                                 gTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGcagcag                                                                                      >  1:1578432/1‑150 (MQ=255)
                                                                 gTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGcagcag                                                                                      <  1:1809531/150‑1 (MQ=255)
                                                                   agagGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGACTCGGCATCGGCAGCagag                                                                                    >  2:2433122/1‑150 (MQ=255)
                                                                   agagGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCagag                                                                                    <  1:3650026/150‑1 (MQ=255)
                                                                    gaggagAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCagaga                                                                                   >  1:3613093/1‑150 (MQ=255)
                                                                        agagagTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTcc                                                                               <  1:2164386/150‑1 (MQ=255)
                                                                        agagagTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTcc                                                                               <  2:2188189/150‑1 (MQ=255)
                                                                         gagagTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCg                                                                              >  2:2667095/1‑150 (MQ=255)
                                                                          agagTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGa                                                                             >  2:2955355/1‑150 (MQ=255)
                                                                            agTAATTATCCGCGGCTTATAATGTCCGCTAATTGTTTTTCTGCTGATGCTGTTAAGCGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATg                                                                           <  2:4411187/150‑1 (MQ=255)
                                                                            agTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATg                                                                           <  1:3227984/150‑1 (MQ=255)
                                                                               aaTTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGt                                                                        <  1:1021708/150‑1 (MQ=255)
                                                                                 ttATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGtaa                                                                      <  2:707696/150‑1 (MQ=255)
                                                                                   aTCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGtaatt                                                                    >  1:2917131/1‑150 (MQ=255)
                                                                                   aTCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGtaatt                                                                    <  2:2561106/150‑1 (MQ=255)
                                                                                      cgcgGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGtaatttaa                                                                 >  2:677537/1‑150 (MQ=255)
                                                                                       gcgGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGtaatttaat                                                                <  2:3705310/150‑1 (MQ=255)
                                                                                        cgGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGtaatttaatt                                                               <  2:754762/150‑1 (MQ=255)
                                                                                        cgGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGtaatttaatt                                                               >  1:4206277/1‑150 (MQ=255)
                                                                                         ggCTGATAATGTCCGCTATTTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTg                                                              <  2:2113239/150‑1 (MQ=255)
                                                                                          gCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTgg                                                             >  2:3003728/1‑150 (MQ=255)
                                                                                            tGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGtt                                                           >  2:3505571/1‑150 (MQ=255)
                                                                                            tGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGtt                                                           >  2:3446872/1‑150 (MQ=255)
                                                                                            tGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGtt                                                           >  2:4079190/1‑150 (MQ=255)
                                                                                             gATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTc                                                          >  2:4018453/1‑150 (MQ=255)
                                                                                             gATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTc                                                          >  1:1216110/1‑150 (MQ=255)
                                                                                              aTAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCg                                                         >  1:289824/1‑150 (MQ=255)
                                                                                              aTAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCg                                                         >  1:3470158/1‑150 (MQ=255)
                                                                                               tAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGt                                                        >  2:725422/1‑150 (MQ=255)
                                                                                                aaTGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTc                                                       <  2:2616412/150‑1 (MQ=255)
                                                                                                      cGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTccc                                                 >  1:2479661/1‑150 (MQ=255)
                                                                                                         tAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCTAGTATTTCTTATTCGCACAGTAATTCACTTAACGCCTGAATAGCTCCAATTTGCGAATCTGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCgg                                              >  2:23101/1‑150 (MQ=255)
                                                                                                              gTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATATCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAAt                                         >  2:604365/1‑150 (MQ=255)
                                                                                                              gTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAAt                                         >  2:2697673/1‑150 (MQ=255)
                                                                                                              gTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAAt                                         >  2:2889538/1‑150 (MQ=255)
                                                                                                                   tCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAaca                                    <  1:3394475/150‑1 (MQ=255)
                                                                                                                       atgatgCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATTTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCtt                                <  2:1472927/150‑1 (MQ=255)
                                                                                                                        tgatgCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGTGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTg                               >  1:3477328/1‑150 (MQ=255)
                                                                                                                        tgatgCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTg                               >  1:988521/1‑150 (MQ=255)
                                                                                                                             cTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCaa                          >  1:2444918/1‑150 (MQ=255)
                                                                                                                                ttAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTa                       <  1:4489849/150‑1 (MQ=255)
                                                                                                                                     agTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGACAATAAACACCTTGTTCAAGTAATGTc                  >  2:3769373/1‑150 (MQ=255)
                                                                                                                                     agTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTc                  >  2:4212236/1‑150 (MQ=255)
                                                                                                                                     agTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTc                  >  2:1921391/1‑150 (MQ=255)
                                                                                                                                      gttgttCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCa                 <  2:2784309/150‑1 (MQ=255)
                                                                                                                                      gttgttCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCa                 <  2:2783756/150‑1 (MQ=255)
                                                                                                                                       ttgttCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCaa                >  2:1909588/1‑150 (MQ=255)
                                                                                                                                        tgttCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCaaa               >  2:2728826/1‑150 (MQ=255)
                                                                                                                                         gttCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAg              >  1:2685385/1‑150 (MQ=255)
                                                                                                                                         gttCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAg              >  1:4116830/1‑150 (MQ=255)
                                                                                                                                             gAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACAt          >  1:2044712/1‑150 (MQ=255)
                                                                                                                                              aGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATc         <  1:1483390/150‑1 (MQ=255)
                                                                                                                                               gTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATcc        <  1:1206720/150‑1 (MQ=255)
                                                                                                                                                     cttcttcGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGttt  <  2:987278/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTT  >  CP009273/3072165‑3072463

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: