Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,072,314 C→A 6.5% 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,3140CA 6.5% 1015.7 / 32.2 272E126* (GAA→TAA) cmtBputative mannitol‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base C (151/102);  new base A (7/12);  total (158/114)
Fisher's exact test for biased strand distribution p-value = 5.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 289 total aligned reads displayed.

GTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTT  >  CP009273/3072166‑3072463
                                                                                                                                                    |                                                                                                                                                     
gTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttct                                                                                                                                                      <  1:650354/150‑1 (MQ=255)
 tCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttctt                                                                                                                                                     <  2:2349060/150‑1 (MQ=255)
  cAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttcttc                                                                                                                                                    >  2:223177/1‑150 (MQ=255)
   aaaaaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttattct                                                                                                                                                   <  2:4271998/150‑1 (MQ=255)
    aaaaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttcttctt                                                                                                                                                  <  1:4484392/150‑1 (MQ=255)
      aaGCACCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTACATAATATTATCCGGTAGAGGAGAGAGTAATTATCCGCGGCTTATAATGTCCGATAATTTTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCg                                                                                                                                                <  2:973005/150‑1 (MQ=255)
       aGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGc                                                                                                                                               >  1:3095486/1‑150 (MQ=255)
        gCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGca                                                                                                                                              <  1:3755804/150‑1 (MQ=255)
        gCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGca                                                                                                                                              <  2:647269/150‑1 (MQ=255)
           cccAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAg                                                                                                                                           >  1:2103149/1‑150 (MQ=255)
            ccAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGt                                                                                                                                          >  2:4198608/1‑150 (MQ=255)
             cAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTa                                                                                                                                         >  1:774530/1‑150 (MQ=255)
              aaaaGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTaa                                                                                                                                        >  2:12603/1‑150 (MQ=255)
               aaaGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAAt                                                                                                                                       <  1:2589786/150‑1 (MQ=255)
                 aGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATTATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTc                                                                                                                                     <  1:1386754/150‑1 (MQ=255)
                   ccTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCAc                                                                                                                                   >  2:534522/1‑150 (MQ=255)
                   ccTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGCACAGTAATTCAc                                                                                                                                   >  2:756474/1‑150 (MQ=255)
                    cTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACt                                                                                                                                  <  1:954505/150‑1 (MQ=255)
                     tGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACtt                                                                                                                                 <  1:4605196/150‑1 (MQ=255)
                          cTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGc                                                                                                                            >  2:2965443/1‑150 (MQ=255)
                           tttGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGcc                                                                                                                           <  1:458228/150‑1 (MQ=255)
                                   aGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAgc                                                                                                                   <  2:3432608/150‑1 (MQ=255)
                                    gCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAgcg                                                                                                                  <  1:1496890/150‑1 (MQ=255)
                                     cAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAgcgc                                                                                                                 >  2:3420485/1‑150 (MQ=255)
                                      aGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGcc                                                                                                                >  1:2288267/1‑150 (MQ=255)
                                           tGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAAtgt                                                                                                           >  2:284318/1‑150 (MQ=255)
                                           tGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAAtgt                                                                                                           <  2:4064828/150‑1 (MQ=255)
                                              tttCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCg                                                                                                        <  2:3149890/150‑1 (MQ=255)
                                                 ccATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAat                                                                                                     >  2:1479565/1‑150 (MQ=255)
                                                 ccATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAat                                                                                                     >  2:1218199/1‑150 (MQ=255)
                                                   ataataTTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCGTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAatcg                                                                                                   <  2:1230702/150‑1 (MQ=255)
                                                        attatCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAatcggcatc                                                                                              >  1:929344/1‑150 (MQ=255)
                                                         ttatCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAatcggcatcg                                                                                             >  1:3407597/1‑150 (MQ=255)
                                                          tatCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAatcggcatcgg                                                                                            >  2:32251/1‑150 (MQ=255)
                                                                  agagGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCagag                                                                                    >  1:2945806/1‑150 (MQ=255)
                                                                   gaggagAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCagaga                                                                                   >  1:3074356/1‑150 (MQ=255)
                                                                    aggagAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCagagag                                                                                  >  1:4308559/1‑150 (MQ=255)
                                                                     ggagAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGt                                                                                 >  2:1733534/1‑150 (MQ=255)
                                                                     ggagAGAGTAATTATCCGCGGCTGAAAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGt                                                                                 <  2:1378589/150‑1 (MQ=255)
                                                                      gagagaGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTc                                                                                >  1:3230681/1‑150 (MQ=255)
                                                                      gagagaGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTc                                                                                >  1:4673953/1‑150 (MQ=255)
                                                                            gTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGa                                                                          >  2:116152/1‑150 (MQ=255)
                                                                             tAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAg                                                                         <  2:3060126/150‑1 (MQ=255)
                                                                             tAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAg                                                                         <  1:1870144/150‑1 (MQ=255)
                                                                              aaTTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGt                                                                        <  1:4317200/150‑1 (MQ=255)
                                                                              aaTTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGt                                                                        <  1:1278279/150‑1 (MQ=255)
                                                                              aaTTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGt                                                                        >  2:2776907/1‑150 (MQ=255)
                                                                               aTTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGta                                                                       >  1:1025771/1‑150 (MQ=255)
                                                                                ttATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGAGGAGtaa                                                                      >  2:3793906/1‑150 (MQ=255)
                                                                                   tCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGtaattt                                                                   <  2:1266578/150‑1 (MQ=255)
                                                                                       cgGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGtaatttaatt                                                               <  1:2738407/150‑1 (MQ=255)
                                                                                       cgGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGtaatttaatt                                                               >  2:649223/1‑150 (MQ=255)
                                                                                        ggCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTg                                                              >  1:634776/1‑150 (MQ=255)
                                                                                           tGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGtt                                                           >  2:374173/1‑150 (MQ=255)
                                                                                           tGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGtt                                                           <  1:3023183/150‑1 (MQ=255)
                                                                                           tGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGtt                                                           >  1:3880221/1‑150 (MQ=255)
                                                                                            gATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTc                                                          >  1:4714138/1‑150 (MQ=255)
                                                                                             aTAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGttag                                                         >  1:3996256/1‑148 (MQ=255)
                                                                                             aTAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTAACTTAACGCCTGAATATCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCg                                                         >  1:1592036/1‑150 (MQ=255)
                                                                                              tAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGt                                                        >  2:454269/1‑150 (MQ=255)
                                                                                              tAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGt                                                        <  2:1878098/150‑1 (MQ=255)
                                                                                               aaTGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTc                                                       >  2:331894/1‑150 (MQ=255)
                                                                                                aTGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTTGTTCGTCa                                                      >  1:1659786/1‑150 (MQ=255)
                                                                                                aTGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGACa                                                      <  1:2911673/150‑1 (MQ=255)
                                                                                                aTGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGACa                                                      <  1:1822682/150‑1 (MQ=255)
                                                                                                  gTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCAtt                                                    <  2:51313/150‑1 (MQ=255)
                                                                                                    ccGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTcc                                                  >  1:1003455/1‑150 (MQ=255)
                                                                                                     cGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTccc                                                 >  1:3064048/1‑150 (MQ=255)
                                                                                                     cGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTccc                                                 >  1:394049/1‑150 (MQ=255)
                                                                                                         aaTTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGa                                             >  1:241167/1‑150 (MQ=255)
                                                                                                         aaTTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGa                                             >  1:2578910/1‑150 (MQ=255)
                                                                                                          aTTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGaa                                            <  2:3578441/150‑1 (MQ=255)
                                                                                                          aTTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGaa                                            <  2:2941058/150‑1 (MQ=255)
                                                                                                           ttGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGaaa                                           <  1:1107972/150‑1 (MQ=255)
                                                                                                            tGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGaaaa                                          <  1:2639377/150‑1 (MQ=255)
                                                                                                             gTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAAt                                         <  1:3822383/150‑1 (MQ=255)
                                                                                                                 ttCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTTGTTCGTCATTCCCCGGAAAATaaaa                                     >  1:2805017/1‑149 (MQ=255)
                                                                                                                 ttCTGATGATGATTTTAAGAGTTGTTCTAGTATTTCTTCTTCGCAGAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGCTGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAac                                     <  1:2027200/150‑1 (MQ=255)
                                                                                                                  tCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAaca                                    <  1:3726739/150‑1 (MQ=255)
                                                                                                                   cTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAacac                                   >  2:131087/1‑150 (MQ=255)
                                                                                                                    tgatgaTGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACAcc                                  <  1:4379240/150‑1 (MQ=255)
                                                                                                                     gatgatGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCt                                 <  1:309310/150‑1 (MQ=255)
                                                                                                                      atgatgCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCtt                                <  2:3953832/150‑1 (MQ=255)
                                                                                                                        gatgCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGt                              >  1:2795855/1‑150 (MQ=255)
                                                                                                                        gatgCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGt                              >  1:567482/1‑150 (MQ=255)
                                                                                                                         atgCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATATCGCCAATGTGCGAATCGGCATCGGCAGCATAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGtt                             >  2:1022916/1‑150 (MQ=255)
                                                                                                                            cTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCaa                          <  2:1565413/150‑1 (MQ=255)
                                                                                                                             tGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAg                         >  1:4115498/1‑150 (MQ=255)
                                                                                                                             tGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCATGTTCAAg                         >  1:1235913/1‑150 (MQ=255)
                                                                                                                               ttAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTa                       <  2:3915391/150‑1 (MQ=255)
                                                                                                                                   gagTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGt                   <  2:2105339/150‑1 (MQ=255)
                                                                                                                                   gagTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGt                   <  2:2618664/150‑1 (MQ=255)
                                                                                                                                         ttCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGa             >  1:1971966/1‑150 (MQ=255)
                                                                                                                                         ttCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGa             >  2:2514983/1‑150 (MQ=255)
                                                                                                                                         ttCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGa             <  1:2047074/150‑1 (MQ=255)
                                                                                                                                            gAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACAt          <  1:2513385/150‑1 (MQ=255)
                                                                                                                                             aGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATc         >  2:4741466/1‑150 (MQ=255)
                                                                                                                                                aTTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCg      >  1:4032170/1‑150 (MQ=255)
                                                                                                                                                  ttcttcttcGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGt    >  2:4026780/1‑150 (MQ=255)
                                                                                                                                                   tcttcttcGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGtt   >  1:346186/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGGTTT  >  CP009273/3072166‑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: