Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,234,3440TG46.3% 47.7 / 16.5 41S348A (TCT→GCT) sstTserine/threonine transporter SstT
Reads supporting (aligned to +/- strand):  ref base T (5/17);  new base G (19/0);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 3.15e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATC  >  NZ_CP009273/3234255‑3234430
                                                                                         |                                                                                      
tGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGt                                                                                        >  2:68055/1‑90 (MQ=255)
       tCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTctgct                                                                                 >  1:498424/1‑90 (MQ=255)
           gTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTctgctgctg                                                                             <  2:337381/90‑1 (MQ=255)
           gTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTctgctgctg                                                                             <  2:498424/90‑1 (MQ=255)
           gTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTctgctgctg                                                                             >  2:278205/1‑90 (MQ=255)
           gTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTctgctgctg                                                                             >  1:337381/1‑90 (MQ=255)
           gTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGGGGCGCATCCGGCGTGGCGGGGGGGGCCctgcggctg                                                                             >  2:495411/1‑90 (MQ=255)
             cGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTCTGCTGCTGAt                                                                           >  2:160919/1‑90 (MQ=255)
               aTCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATcc                                                                         >  1:263663/1‑90 (MQ=255)
                    cccACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTg                                                                    <  1:160919/90‑1 (MQ=255)
                       aCGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGcc                                                                 <  2:263663/90‑1 (MQ=255)
                       aCGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGcc                                                                 <  2:382590/90‑1 (MQ=255)
                       aCGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGcc                                                                 <  1:807091/90‑1 (MQ=255)
                       aCGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGcc                                                                 <  1:68055/90‑1 (MQ=255)
                        cGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTCTGCTGCTGATCCCCCTTGCCt                                                                >  1:862000/1‑90 (MQ=255)
                             ctgctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGtc                                                                                       <  2:857903/62‑1 (MQ=255)
                             ctgctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGtc                                                                                       <  1:76029/62‑1 (MQ=255)
                             ctgctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGtc                                                                                       >  1:857903/1‑62 (MQ=255)
                             ctgctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGcc                                                                 <  2:931355/84‑1 (MQ=255)
                             ctgctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTCTGCTGCTGATCCCACTGGcc                                                                 >  1:931355/1‑84 (MQ=255)
                                    tGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCAc                                                                      <  2:811021/72‑1 (MQ=255)
                                    tGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCAc                                                                      <  2:423269/72‑1 (MQ=255)
                                    tGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTCTGCTGCTGATCCCAc                                                                      >  1:811021/1‑72 (MQ=255)
                                    tGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTCTGCTGCTGATCCCAc                                                                      >  1:423269/1‑72 (MQ=255)
                                     gAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCgg                                                   <  2:833162/90‑1 (MQ=255)
                                       gCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCCCTGGCCTGTAATATGTTTggga                                                 >  1:497778/1‑88 (MQ=255)
                                         gtggtgGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCAc                                                                      >  1:811492/1‑67 (MQ=255)
                                         gtggtgGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCAc                                                                      <  2:811492/67‑1 (MQ=255)
                                         gtggtgGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGGTCTGCTGCTGATCCCCCCGGCCTGTTATATGTTTggggtt                                               >  2:720821/1‑86 (MQ=255)
                                         gtggtgGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTCTGCTGCTGTTTCCCCCGGCCTGTAATATGTTCggggtt                                               >  2:883122/1‑86 (MQ=255)
                                         gtggtgGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTTGGGAtt                                               >  1:721302/1‑90 (MQ=255)
                                            gtggCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGtc                                                                                       <  1:378640/47‑1 (MQ=255)
                                            gtggCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGtc                                                                                       <  1:1013962/47‑1 (MQ=255)
                                            gtggCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTCTGCTGCTGATCCCCCTGGCCTGTAATATGTTTGGGATTTCg                                            >  2:189051/1‑90 (MQ=255)
                                             tggCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGa                                           <  1:883122/90‑1 (MQ=255)
                                                 ttCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTCTGCTGCTGATCCCCCTGGCCTGTAATATGTTCGGGATTTTGAACGa                                       >  2:258891/1‑90 (MQ=255)
                                                      tgtgtgCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGGATTTCGAACGATATCg                                  >  1:714421/1‑90 (MQ=255)
                                                        tgtgCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGcc                                                                 <  2:666454/57‑1 (MQ=255)
                                                        tgtgCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGcc                                                                 <  1:647762/57‑1 (MQ=255)
                                                        tgtgCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTCTGCTGCTGATCCCCCTGGcc                                                                 >  2:647762/1‑57 (MQ=255)
                                                        tgtgCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTCTGCTGCTGATCCCACTGGcc                                                                 >  1:666454/1‑57 (MQ=255)
                                                               gtgGCGCATCCGGCGTGGCGGGGGGGGCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAgg                         >  1:952645/1‑90 (MQ=255)
                                                                 ggCGCATCCGGCGTGGCGGGGGGGGCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAggtg                       >  1:426083/1‑90 (MQ=255)
                                                                   cgcATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGGATTTCGAACGATATCGCCATGCAggtggt                     >  1:419344/1‑90 (MQ=255)
                                                                      aTCCGGCGTGGCGGGGGGGGCTCTGCTGCTGATCCCACTGGCCTTTAATATGTTCGGGATTTTGAACGATATCGCCCTGCAGGGGGGTGc                  >  2:122816/1‑90 (MQ=255)
                                                                                    ggggggCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGGGGGTGCCGTCGGCTTTttca    >  1:139466/1‑90 (MQ=255)
                                                                                      gggTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTatcatc  <  2:524416/90‑1 (MQ=255)
                                                                                      gggTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTatcatc  <  2:95453/90‑1 (MQ=255)
                                                                                         |                                                                                      
TGGGTATTCCGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATC  >  NZ_CP009273/3234255‑3234430

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: