Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,072,314 C→A 14.2% 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,3140CA14.2% 824.8 / 88.7 260E126* (GAA→TAA) cmtBputative mannitol‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base C (106/117);  new base A (19/18);  total (125/135)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 282 total aligned reads displayed.

AGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGG  >  CP009273/3072165‑3072460
                                                                                                                                                     |                                                                                                                                                  
aGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttc                                                                                                                                                    >  2:798490/1‑150 (MQ=255)
 gTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttct                                                                                                                                                   <  2:2867985/150‑1 (MQ=255)
 gTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttct                                                                                                                                                   >  2:2388187/1‑150 (MQ=255)
  tCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGtatttatt                                                                                                                                                  <  2:315218/150‑1 (MQ=255)
   cAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttcttc                                                                                                                                                 >  1:409712/1‑150 (MQ=255)
    aaaaaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttcttct                                                                                                                                                <  1:1723585/150‑1 (MQ=255)
    aaaaaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttcttct                                                                                                                                                <  2:3277207/150‑1 (MQ=255)
     aaaaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttcttctt                                                                                                                                               <  2:698081/150‑1 (MQ=255)
     aaaaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttcttctt                                                                                                                                               <  1:991410/150‑1 (MQ=255)
     aaaaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGAGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttcttctt                                                                                                                                               <  1:3816076/150‑1 (MQ=255)
      aaaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATttcttcttc                                                                                                                                              <  2:2172517/150‑1 (MQ=255)
       aaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCg                                                                                                                                             <  1:483954/150‑1 (MQ=255)
       aaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCg                                                                                                                                             <  1:3052229/150‑1 (MQ=255)
       aaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCg                                                                                                                                             >  2:1315687/1‑150 (MQ=255)
       aaGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCAATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCg                                                                                                                                             <  2:1606624/150‑1 (MQ=255)
        aGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGc                                                                                                                                            >  2:3807024/1‑150 (MQ=255)
        aGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGc                                                                                                                                            >  1:1939518/1‑150 (MQ=255)
            cccAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAg                                                                                                                                        <  1:4015705/150‑1 (MQ=255)
              cAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTa                                                                                                                                      <  2:864528/150‑1 (MQ=255)
                      tGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACtt                                                                                                                              <  1:2391792/150‑1 (MQ=255)
                      tGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACtt                                                                                                                              <  1:1789492/150‑1 (MQ=255)
                      tGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACtt                                                                                                                              <  2:1561993/150‑1 (MQ=255)
                       ggACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTATTATTCGAACAGTAATTAACTTa                                                                                                                             >  2:3508841/1‑150 (MQ=255)
                        gACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGTGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTaa                                                                                                                            <  1:4072670/150‑1 (MQ=255)
                        gACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGCACAGTAATTCACTTaa                                                                                                                            >  2:2069702/1‑150 (MQ=255)
                         aCCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAAc                                                                                                                           >  2:3980546/1‑150 (MQ=255)
                         aCCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAAc                                                                                                                           >  1:4318251/1‑150 (MQ=255)
                             ttGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCt                                                                                                                       <  1:3476371/150‑1 (MQ=255)
                                  cGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATa                                                                                                                  <  2:1855372/150‑1 (MQ=255)
                                   gAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAg                                                                                                                 >  2:3253591/1‑150 (MQ=255)
                                     gCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAgcg                                                                                                               <  1:3922253/150‑1 (MQ=255)
                                           ttGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAAtg                                                                                                         <  2:595172/150‑1 (MQ=255)
                                                   cATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAatc                                                                                                 <  2:4119180/150‑1 (MQ=255)
                                                   cATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAatc                                                                                                 <  1:862258/150‑1 (MQ=255)
                                                   cATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAatc                                                                                                 <  1:1011060/150‑1 (MQ=255)
                                                     taataTTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAatcgg                                                                                               <  2:3019066/150‑1 (MQ=255)
                                                        tattatCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAatcggcat                                                                                            >  2:3564225/1‑150 (MQ=255)
                                                          ttatCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAatcggcatcg                                                                                          >  2:163983/1‑150 (MQ=255)
                                                           tatCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAatcggcatcgg                                                                                         <  1:3349559/150‑1 (MQ=255)
                                                              ccTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGgcag                                                                                      >  1:500062/1‑150 (MQ=255)
                                                              ccTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGgcag                                                                                      <  1:522102/150‑1 (MQ=255)
                                                              ccTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGgcag                                                                                      >  2:2302867/1‑150 (MQ=255)
                                                               cTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGgcagc                                                                                     >  2:2073087/1‑150 (MQ=255)
                                                                tGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGgcagca                                                                                    >  1:781028/1‑150 (MQ=255)
                                                                 gTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGcagcag                                                                                   >  1:1705633/1‑150 (MQ=255)
                                                                   agagGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCagag                                                                                 >  1:1079474/1‑150 (MQ=255)
                                                                    gaggagAGAGTAATTATCCGCGGCTGAGAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCagaga                                                                                <  2:4413803/150‑1 (MQ=255)
                                                                      ggagAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGt                                                                              <  2:4058579/150‑1 (MQ=255)
                                                                        agagagTAATTATCCGCGGCTGATAAGGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTcc                                                                            <  2:1910301/150‑1 (MQ=255)
                                                                          agagTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGa                                                                          <  2:3244929/150‑1 (MQ=255)
                                                                           gagTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGAt                                                                         <  1:1402298/150‑1 (MQ=255)
                                                                           gagTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGAt                                                                         <  1:4330390/150‑1 (MQ=255)
                                                                            agTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATg                                                                        <  1:1668667/150‑1 (MQ=255)
                                                                            agTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATg                                                                        <  1:2624562/150‑1 (MQ=255)
                                                                                   aTCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGtaatt                                                                 <  2:1726386/150‑1 (MQ=255)
                                                                                    tCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGtaattt                                                                <  1:2946887/150‑1 (MQ=255)
                                                                                     ccGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGtaattta                                                               >  1:2860162/1‑150 (MQ=255)
                                                                                       gcgGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGtaatttaat                                                             <  1:3864342/150‑1 (MQ=255)
                                                                                       gcgGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGtaatttaaa                                                             <  2:2532968/150‑2 (MQ=255)
                                                                                        cgGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGtaatttaatt                                                            <  2:1088725/150‑1 (MQ=255)
                                                                                        cgGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGtaatttaatt                                                            >  2:3602127/1‑150 (MQ=255)
                                                                                        cgGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGtaatttaatt                                                            <  2:1385645/150‑1 (MQ=255)
                                                                                         ggCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTg                                                           <  2:4275434/150‑1 (MQ=255)
                                                                                          gCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTgg                                                          >  2:1191396/1‑150 (MQ=255)
                                                                                           cTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGt                                                         >  1:3680408/1‑150 (MQ=255)
                                                                                             gATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTc                                                       <  2:363316/150‑1 (MQ=255)
                                                                                              aTAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCg                                                      >  2:1610253/1‑150 (MQ=255)
                                                                                              aTAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCg                                                      >  2:3419355/1‑150 (MQ=255)
                                                                                                aaTGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTc                                                    >  2:2155233/1‑150 (MQ=255)
                                                                                                 aTGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCa                                                   >  2:944320/1‑150 (MQ=255)
                                                                                                  tGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCAt                                                  <  2:2232071/150‑1 (MQ=255)
                                                                                                  tGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCAt                                                  <  1:645332/150‑1 (MQ=255)
                                                                                                     ccGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTcc                                               <  1:2762101/150‑1 (MQ=255)
                                                                                                       gCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTcccc                                             >  1:4417672/1‑150 (MQ=255)
                                                                                                                ttttCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATaa                                    >  1:4338682/1‑150 (MQ=255)
                                                                                                                   tCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAaca                                 >  2:935105/1‑150 (MQ=255)
                                                                                                                       atgatgCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCtt                             <  2:2032135/150‑1 (MQ=255)
                                                                                                                       atgatgCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCtt                             >  2:4347285/1‑150 (MQ=255)
                                                                                                                       atgatgCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCtt                             <  1:1191269/150‑1 (MQ=255)
                                                                                                                         gatgCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGt                           <  1:3254183/150‑1 (MQ=255)
                                                                                                                          atgCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGtt                          >  2:1595833/1‑150 (MQ=255)
                                                                                                                           tgCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTc                         <  2:911424/150‑1 (MQ=255)
                                                                                                                           tgCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTc                         <  1:2137556/150‑1 (MQ=255)
                                                                                                                           tgCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTc                         <  2:623184/150‑1 (MQ=255)
                                                                                                                            gCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCa                        <  2:375944/150‑1 (MQ=255)
                                                                                                                            gCTGTTAAGAGTTGTTCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCa                        <  2:105554/150‑1 (MQ=255)
                                                                                                                              tGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAg                      >  1:3630185/1‑150 (MQ=255)
                                                                                                                                 tAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTaa                   <  2:3877066/150‑1 (MQ=255)
                                                                                                                                  aaGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAAt                  >  2:3228736/1‑150 (MQ=255)
                                                                                                                                   agagTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATg                 <  1:2189685/150‑1 (MQ=255)
                                                                                                                                   agagTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATg                 <  2:2096149/150‑1 (MQ=255)
                                                                                                                                    gagTTGTTCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGt                <  1:2012592/150‑1 (MQ=255)
                                                                                                                                     agTTGTTCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGACAATAAACACCTTGTTCAAGTAATGTc               >  2:2637986/1‑150 (MQ=255)
                                                                                                                                      gttgttCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCACTGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCa              <  2:4211565/150‑1 (MQ=255)
                                                                                                                                      gttgttCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCa              <  2:4043418/150‑1 (MQ=255)
                                                                                                                                       ttgttCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCaa             >  2:4227251/1‑150 (MQ=255)
                                                                                                                                        tgttCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCaaa            >  1:2648722/1‑150 (MQ=255)
                                                                                                                                         gttCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAg           >  1:3348946/1‑150 (MQ=255)
                                                                                                                                          ttCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGa          <  1:3310084/150‑1 (MQ=255)
                                                                                                                                          ttCGAGTATTTATTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGa          >  1:1898768/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AGTCAAAAAGCCCCCAAAAGCCTGGACCTTTGCACGAGCAGACTTGTTTTCCATAATATTATCCTGTAGAGGAGAGAGTAATTATCCGCGGCTGATAATGTCCGCTAATTGTTTTTCTGATGATGCTGTTAAGAGTTGTTCGAGTATTTCTTCTTCGCACAGTAATTCACTTAACGCCTGAATAGCGCCAATGTGCGAATCGGCATCGGCAGCAGAGAGTCCGATGAGTAATTTAATTGGTTCGTCATTCCCCGGAAAATAAACACCTTGTTCAAGTAATGTCAAAGACATCCCGG  >  CP009273/3072165‑3072460

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: