Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,767,000 G→T 100% D454Y (GAC→TAC)  mtlA → fused mannitol‑specific PTS enzymes: IIA components/IIB components/IIC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,767,0000GT98.4% 639.0 / 0.9 202D454Y (GAC→TAC) mtlAfused mannitol‑specific PTS enzymes: IIA components/IIB components/IIC components
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (83/115);  total (85/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 202 total aligned reads displayed.

GTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATC  >  CP009273/3766857‑3767143
                                                                                                                                               |                                                                                                                                               
gTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAg                                                                                                                                             >  2:804146/1‑148 (MQ=255)
gTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAg                                                                                                                                             >  2:3048139/1‑148 (MQ=255)
 tCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGc                                                                                                                                            >  2:679160/1‑148 (MQ=255)
 tCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGc                                                                                                                                            >  2:3004973/1‑148 (MQ=255)
 tCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGc                                                                                                                                            <  1:739626/148‑1 (MQ=255)
    gTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGc                                                                                                                                         >  2:1716434/1‑148 (MQ=255)
    gTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGc                                                                                                                                         >  1:3439637/1‑148 (MQ=255)
    gTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGc                                                                                                                                         <  1:2224878/148‑1 (MQ=255)
      cGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCt                                                                                                                                       >  2:90714/1‑148 (MQ=255)
      cGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCt                                                                                                                                       >  2:108343/1‑148 (MQ=255)
      cGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCt                                                                                                                                       >  2:108237/1‑148 (MQ=255)
       gCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTg                                                                                                                                      <  1:1027423/148‑1 (MQ=255)
           aTTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTaca                                                                                                                                  >  2:1493710/1‑148 (MQ=255)
             ttCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACAcc                                                                                                                                >  1:2834114/1‑148 (MQ=255)
                 gTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGcc                                                                                                                            >  1:2106543/1‑148 (MQ=255)
                 gTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGcc                                                                                                                            <  1:2269172/148‑1 (MQ=255)
                 gTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGcc                                                                                                                            <  1:940016/148‑1 (MQ=255)
                   tACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTg                                                                                                                          >  1:1562519/1‑148 (MQ=255)
                        aCAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGa                                                                                                                     <  2:503839/148‑1 (MQ=255)
                            cgcgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGt                                                                                                                 >  2:3811796/1‑148 (MQ=255)
                            cgcgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGt                                                                                                                 >  2:207154/1‑148 (MQ=255)
                            cgcgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGt                                                                                                                 >  1:2829138/1‑148 (MQ=255)
                             gcgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTc                                                                                                                <  2:793049/148‑1 (MQ=255)
                             gcgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTc                                                                                                                <  1:3717217/148‑1 (MQ=255)
                                aTCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTgg                                                                                                             >  1:330463/1‑148 (MQ=255)
                                aTCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTgg                                                                                                             <  2:3317401/148‑1 (MQ=255)
                                   aacaacCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTtg                                                                                                          <  2:2454726/148‑1 (MQ=255)
                                       acCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTtgctgc                                                                                                      >  1:2086963/1‑148 (MQ=255)
                                               ccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGcc                                                                                              <  2:418626/148‑1 (MQ=255)
                                                 aGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCcac                                                                                            <  1:1147492/148‑1 (MQ=255)
                                                  gATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCcaca                                                                                           <  1:2678475/148‑1 (MQ=255)
                                                  gATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCcaca                                                                                           >  2:2106255/1‑148 (MQ=255)
                                                    tgtgGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACg                                                                                         <  1:2473090/148‑1 (MQ=255)
                                                    tgtgGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACg                                                                                         <  2:253069/148‑1 (MQ=255)
                                                        gACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCaa                                                                                     <  1:1244739/148‑1 (MQ=255)
                                                        gACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCaa                                                                                     >  2:3502457/1‑148 (MQ=255)
                                                        gACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCaa                                                                                     >  2:1736526/1‑148 (MQ=255)
                                                         aCCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCaaa                                                                                    <  1:540491/148‑1 (MQ=255)
                                                          ccTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAAc                                                                                   >  2:2156170/1‑148 (MQ=255)
                                                          ccTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAAc                                                                                   >  2:227234/1‑148 (MQ=255)
                                                            tCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGa                                                                                 >  1:3717479/1‑148 (MQ=255)
                                                            tCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGa                                                                                 >  1:2804312/1‑148 (MQ=255)
                                                            tCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGa                                                                                 >  1:3379417/1‑148 (MQ=255)
                                                              gTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAag                                                                               <  1:3450583/148‑1 (MQ=255)
                                                              gTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAag                                                                               <  2:507315/148‑1 (MQ=255)
                                                              gTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAag                                                                               <  1:2390734/148‑1 (MQ=255)
                                                               tcatcaCTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAaga                                                                              <  2:3735022/148‑1 (MQ=255)
                                                                catcaCTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAagag                                                                             <  1:858821/148‑1 (MQ=255)
                                                                catcaCTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAagag                                                                             <  1:3274722/148‑1 (MQ=255)
                                                                  tcactcacCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGaa                                                                           >  1:1057823/1‑148 (MQ=255)
                                                                  tcactcacCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGaa                                                                           <  1:1541476/148‑1 (MQ=255)
                                                                     ctcacCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGt                                                                        >  1:2734119/1‑148 (MQ=255)
                                                                        acCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTaaa                                                                     <  2:2973452/148‑1 (MQ=255)
                                                                        acCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTaaa                                                                     >  1:2771864/1‑148 (MQ=255)
                                                                        acCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTaaa                                                                     <  2:1726870/148‑1 (MQ=255)
                                                                          cGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAg                                                                   <  2:2829138/148‑1 (MQ=255)
                                                                          cGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAg                                                                   >  2:2986420/1‑148 (MQ=255)
                                                                          cGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAg                                                                   <  1:3783351/148‑1 (MQ=255)
                                                                           gTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGa                                                                  <  1:135063/148‑1 (MQ=255)
                                                                           gTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGa                                                                  <  1:3085923/148‑1 (MQ=255)
                                                                           gTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATAGTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGa                                                                  <  2:497896/148‑1 (MQ=255)
                                                                            tgacctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGAc                                                                 <  2:752398/148‑1 (MQ=255)
                                                                              acctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAg                                                               >  1:2268562/1‑148 (MQ=255)
                                                                              acctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAg                                                               >  1:526506/1‑148 (MQ=255)
                                                                                     cGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCACACGAAGAGACAGTAAAAGACAGCCTGaaa                                                        >  2:1087201/1‑148 (MQ=255)
                                                                                     cGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGaaa                                                        >  2:1087580/1‑148 (MQ=255)
                                                                                     cGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGaaa                                                        <  1:544844/148‑1 (MQ=255)
                                                                                      gAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAg                                                       <  2:1998808/148‑1 (MQ=255)
                                                                                       aaCGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGa                                                      >  1:2219985/1‑148 (MQ=255)
                                                                                       aaCGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGa                                                      <  2:2924755/148‑1 (MQ=255)
                                                                                        ccgcgcgATGCGCCAGGGTCCGCAGGCCCAGCAGAGGACGCCGACCCCCGTCCTCTCCAGCGGCCCCTACACCAGGCTGACCGATCGTCTGGGTGCTGCCCAACGCCACACCGCAATCGAAGAGAAAGTAAAAGACACCCTGAAAGAc                                                     <  2:2582543/147‑1 (MQ=255)
                                                                                        aCGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGAc                                                     >  2:123699/1‑148 (MQ=255)
                                                                                           cgcTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGc                                                  <  2:2736456/148‑1 (MQ=255)
                                                                                                  cgcCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACg                                           <  1:3639121/148‑1 (MQ=255)
                                                                                                       ggTTCCGCAGGCACAGCATGTTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTcc                                      >  1:2288019/1‑148 (MQ=255)
                                                                                                       ggTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTcc                                      <  1:1246946/148‑1 (MQ=255)
                                                                                                       ggTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTcc                                      >  1:819355/1‑148 (MQ=255)
                                                                                                           ccGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTg                                  <  1:3003545/148‑1 (MQ=255)
                                                                                                           ccGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTg                                  <  1:679160/148‑1 (MQ=255)
                                                                                                            cGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGc                                 <  1:1006995/148‑1 (MQ=255)
                                                                                                            cGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGc                                 <  1:108259/148‑1 (MQ=255)
                                                                                                              cAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTa                               <  2:1631420/148‑1 (MQ=255)
                                                                                                              cAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTa                               <  2:3526800/148‑1 (MQ=255)
                                                                                                               aGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTaa                              <  1:1320074/148‑1 (MQ=255)
                                                                                                                ggCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAAc                             >  2:1118515/1‑148 (MQ=255)
                                                                                                                ggCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAAc                             >  2:2628932/1‑148 (MQ=255)
                                                                                                                 gCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAAcc                            >  1:1407252/1‑148 (MQ=255)
                                                                                                                 gCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAAcc                            >  1:1160048/1‑148 (MQ=255)
                                                                                                                       cATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCa                      >  1:3180345/1‑148 (MQ=255)
                                                                                                                          atTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGc                   >  2:1184413/1‑148 (MQ=255)
                                                                                                                          atTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGc                   <  1:2000781/148‑1 (MQ=255)
                                                                                                                          atTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGAGCTAACCTGTTCAAGc                   >  2:3188569/1‑148 (MQ=255)
                                                                                                                          atTGCGCTGACCAACTTCCTCGACCGCGGCCTGTACCCCCGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTATCCTGTTCAAGc                   <  2:330463/148‑1 (MQ=255)
                                                                                                                           tttCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCATCGCCACACGGCAAACGTAGAGAAAGTAAAAGACAGCCTTTAATTCAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCt                  <  2:1920469/148‑1 (MQ=255)
                                                                                                                            ttCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTa                 >  2:1407182/1‑148 (MQ=255)
                                                                                                                            ttCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTa                 >  2:688108/1‑148 (MQ=255)
                                                                                                                              cGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAgg               >  1:166523/1‑148 (MQ=255)
                                                                                                                                   aCCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCgg          <  1:3804597/148‑1 (MQ=255)
                                                                                                                                          tCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACAt   <  1:2361410/148‑1 (MQ=255)
                                                                                                                                           ccTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATc  <  1:1512972/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
GTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATC  >  CP009273/3766857‑3767143

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: