Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,105,662 A→T 37.8% W497R (TGG→AGG)  glpK ← glycerol kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,105,6620AT37.8% 292.3 / inf 304W497R (TGG→AGG) glpKglycerol kinase
Reads supporting (aligned to +/- strand):  ref base A (94/95);  new base T (61/54);  total (155/149)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 306 total aligned reads displayed.

CTCTAAGCAGTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGA  >  CP009273/4105514‑4105811
                                                                                                                                                    |                                                                                                                                                     
 tctAAGCAGTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACg                                                                                                                                                     >  1:2199647/1‑150 (MQ=255)
    aaGCAGTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCa                                                                                                                                                  >  1:258597/1‑150 (MQ=255)
     aGCAGTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCAt                                                                                                                                                 >  1:4111687/1‑150 (MQ=255)
     aGCAGTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCAt                                                                                                                                                 >  1:2169225/1‑150 (MQ=255)
      gCAGTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATc                                                                                                                                                <  1:959006/150‑1 (MQ=255)
      gCAGTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATc                                                                                                                                                <  1:3749708/150‑1 (MQ=255)
         gTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATcgcg                                                                                                                                             <  1:3383345/150‑1 (MQ=255)
           caGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATcgcgcg                                                                                                                                           >  2:3991123/2‑150 (MQ=255)
           aaGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATcgcgcg                                                                                                                                           <  1:759265/150‑1 (MQ=255)
           aaGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATcgcgcg                                                                                                                                           <  2:3046059/150‑1 (MQ=255)
               aaTGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTa                                                                                                                                       <  2:2978210/150‑1 (MQ=255)
               aaTGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTa                                                                                                                                       <  2:141903/150‑1 (MQ=255)
               aaTGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTa                                                                                                                                       <  2:2015830/150‑1 (MQ=255)
               aaTGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTa                                                                                                                                       >  2:2754225/1‑150 (MQ=255)
                aTGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTaa                                                                                                                                      >  2:1646063/1‑150 (MQ=255)
                aTGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTaa                                                                                                                                      <  1:530260/150‑1 (MQ=255)
                 tGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAAc                                                                                                                                     <  2:467656/150‑1 (MQ=255)
                      cgAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCtt                                                                                                                                <  1:1618783/150‑1 (MQ=255)
                         aGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCttttt                                                                                                                             >  2:1158116/1‑150 (MQ=255)
                          gtgtAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCtttttt                                                                                                                            <  2:2531600/150‑1 (MQ=255)
                             tAGCACCGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCa                                                                                                                         >  2:3935094/1‑150 (MQ=255)
                             tAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCa                                                                                                                         >  2:3240040/1‑150 (MQ=255)
                             tAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCa                                                                                                                         <  2:359190/150‑1 (MQ=255)
                                 acaGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCt                                                                                                                     <  1:1796781/150‑1 (MQ=255)
                                    ggggAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCg                                                                                                                  >  1:2301499/1‑150 (MQ=255)
                                    ggggAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCg                                                                                                                  >  1:2418873/1‑150 (MQ=255)
                                        aaGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAAc                                                                                                              <  1:1822271/150‑1 (MQ=255)
                                            gagggaTTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTa                                                                                                          <  2:1525780/150‑1 (MQ=255)
                                              gggaTTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAAt                                                                                                        >  1:1554980/1‑150 (MQ=255)
                                              gggaTTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAAt                                                                                                        <  2:2823820/150‑1 (MQ=255)
                                               ggaTTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAAtt                                                                                                       <  1:3133298/150‑1 (MQ=255)
                                                gaTTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTa                                                                                                      <  2:1137096/150‑1 (MQ=255)
                                                 attattGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTAc                                                                                                     <  2:3774758/150‑1 (MQ=255)
                                                 attattGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTAc                                                                                                     <  1:2529357/150‑1 (MQ=255)
                                                  ttattGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACg                                                                                                    <  1:4020843/150‑1 (MQ=255)
                                                   tattGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGc                                                                                                   <  1:3991959/150‑1 (MQ=255)
                                                      tGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCa                                                                                                <  1:2510623/150‑1 (MQ=255)
                                                       gATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAg                                                                                               <  2:3355988/150‑1 (MQ=255)
                                                       gATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAg                                                                                               >  2:3026509/1‑150 (MQ=255)
                                                       gATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAg                                                                                               <  1:1214327/150‑1 (MQ=255)
                                                        aTGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGt                                                                                              <  1:125456/150‑1 (MQ=255)
                                                         tgtgtgCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTCATTACGCTCAGTg                                                                                             >  2:745528/1‑150 (MQ=255)
                                                            gtgCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGtt                                                                                          <  2:1351954/150‑1 (MQ=255)
                                                             tgCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTCATTACGCTCAGTGGttt                                                                                         >  2:2430934/1‑150 (MQ=255)
                                                               cGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCg                                                                                       <  1:2188552/150‑1 (MQ=255)
                                                                ggggTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGa                                                                                      >  2:1243556/1‑150 (MQ=255)
                                                                 gggTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGAt                                                                                     <  1:3942677/150‑1 (MQ=255)
                                                                  ggTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCACAGTGGTTTCGATg                                                                                    <  2:2204056/150‑1 (MQ=255)
                                                                  ggTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATg                                                                                    <  2:1577202/150‑1 (MQ=255)
                                                                  ggTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATg                                                                                    >  2:4021027/1‑150 (MQ=255)
                                                                  ggTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATg                                                                                    >  2:2113281/1‑150 (MQ=255)
                                                                     tGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCt                                                                                 <  1:3413924/150‑1 (MQ=255)
                                                                       cccccGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTgg                                                                               >  1:906514/1‑150 (MQ=255)
                                                                           cGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACgg                                                                           <  2:2587125/150‑1 (MQ=255)
                                                                           cGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACgg                                                                           >  1:2129517/1‑150 (MQ=255)
                                                                               cGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACt                                                                       <  1:2694428/150‑1 (MQ=255)
                                                                                gTTTCGGGACTACCGGATGCGGCATATACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTc                                                                      <  1:3830029/150‑1 (MQ=255)
                                                                                   tCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTcgcg                                                                   <  2:222914/150‑1 (MQ=255)
                                                                                     gggACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGATTTCGATGCCTGGACGGAACTCGCGCt                                                                 <  2:1019299/150‑1 (MQ=255)
                                                                                      ggACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTCTTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTc                                                                >  2:878207/1‑150 (MQ=255)
                                                                                      ggACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTc                                                                <  1:1430416/150‑1 (MQ=255)
                                                                                        aCTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCaa                                                              >  1:2416120/1‑150 (MQ=255)
                                                                                          tACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATc                                                            >  1:2800928/1‑150 (MQ=255)
                                                                                           aCCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCa                                                           >  1:297646/1‑150 (MQ=255)
                                                                                           aCCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCa                                                           <  2:1870750/150‑1 (MQ=255)
                                                                                           aCCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGctcactca                                                           >  1:1745256/1‑150 (MQ=255)
                                                                                             cGGATGCGGCATAAACGCTTCAGTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCAcc                                                         >  2:1459654/1‑150 (MQ=255)
                                                                                              ggATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCg                                                        >  1:323505/1‑150 (MQ=255)
                                                                                                aTGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCACAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCt                                                      <  1:1110317/150‑1 (MQ=255)
                                                                                                 tGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCtt                                                     >  1:2544837/1‑150 (MQ=255)
                                                                                                  gCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCttt                                                    <  2:1855247/150‑1 (MQ=255)
                                                                                                   cGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTtc                                                   <  1:1385487/150‑1 (MQ=255)
                                                                                                    ggCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTtct                                                  <  2:549186/150‑1 (MQ=255)
                                                                                                      cATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTctct                                                >  1:3574151/1‑150 (MQ=255)
                                                                                                       aTAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCtt                                               <  2:1732077/150‑1 (MQ=255)
                                                                                                        tAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTg                                              <  1:3727683/150‑1 (MQ=255)
                                                                                                          aaCGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCa                                            <  1:4114544/150‑1 (MQ=255)
                                                                                                               ttCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCtcg                                       >  1:1924356/1‑150 (MQ=255)
                                                                                                                tCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCtcgt                                      <  2:434008/150‑1 (MQ=255)
                                                                                                                tCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCtcgt                                      >  2:3142720/1‑150 (MQ=255)
                                                                                                                  aTTCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCtcgtcg                                    >  1:4043260/1‑150 (MQ=255)
                                                                                                                  aTTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCtcgtcg                                    >  2:4045231/1‑150 (MQ=255)
                                                                                                                   ttCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACACTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGa                                   >  2:1384667/1‑150 (MQ=255)
                                                                                                                   ttCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGa                                   <  2:2965277/150‑1 (MQ=255)
                                                                                                                    tCGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAg                                  >  2:638802/1‑150 (MQ=255)
                                                                                                                     cGGCATTTACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAgg                                 <  1:4170792/150‑1 (MQ=255)
                                                                                                                         aTTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCt                             <  2:102964/150‑1 (MQ=255)
                                                                                                                         aTTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCt                             <  1:2760735/150‑1 (MQ=255)
                                                                                                                          tttACATTATTCGTCGTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTg                            >  1:3079278/1‑150 (MQ=255)
                                                                                                                            tACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGcc                          <  2:3694131/150‑1 (MQ=255)
                                                                                                                            tACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGcc                          >  2:3302726/1‑150 (MQ=255)
                                                                                                                              cATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAg                        >  1:3974283/1‑150 (MQ=255)
                                                                                                                                    tcgtcgTGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCa                  >  2:3656600/1‑150 (MQ=255)
                                                                                                                                       tcgtGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAAcc               >  1:519161/1‑150 (MQ=255)
                                                                                                                                        cgtGTTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCg              <  1:3781825/150‑1 (MQ=255)
                                                                                                                                           gtTCTTCCCTCGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCa           >  2:3566234/1‑150 (MQ=255)
                                                                                                                                           gtTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCa           >  2:2370321/1‑150 (MQ=255)
                                                                                                                                             tcttcCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAgg         >  1:2034828/1‑150 (MQ=255)
                                                                                                                                                tcCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCt      <  1:1550824/150‑1 (MQ=255)
                                                                                                                                                  ccACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGc    <  1:824907/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTCTAAGCAGTAAGGAATGGCGCGAAGTGTAGCACAGGGGAAGGGAGGGATTATTGATGTGTGCGGGGTTGCCCCCGCACGTTTCGGGACTACCGGATGCGGCATAAACGCTTCATTCGGCATTTACATTATTCGTCGTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGA  >  CP009273/4105514‑4105811

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: