Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,131,4400TC43.2% 50.2 / 10.1 37G323G (GGA→GGG) yghQputative inner membrane polysaccharide flippase
Reads supporting (aligned to +/- strand):  ref base T (19/2);  new base C (0/16);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 1.19e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTT  >  NC_000913/3131292‑3131583
                                                                                                                                                     |                                                                                                                                               
aaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTcc                                                                                                                                               <  2:278959/151‑1 (MQ=255)
 aatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTcc                                                                                                                                               >  1:1045568/1‑150 (MQ=255)
        gAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCa                                                                                                                                       <  1:238589/151‑1 (MQ=255)
                cGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCg                                                                                                                               >  1:911733/1‑151 (MQ=255)
                cGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCg                                                                                                                               >  1:73689/1‑151 (MQ=255)
                               cGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTaa                                                                                                                 >  1:613499/1‑150 (MQ=255)
                                  ttATTCCCCCCTTCCCAAAAGCAGGTCATCCGCTTCGAGATACTAAACGCAAAC‑CACCAGAGAAAGGAGCGGCTTGCCGACGATCACCCCCGCGAGCCCCACAAGAATACCGTTCCCGCCCGCCAGTAAACCGGACTTCACACCTaac                                                                                                                <  2:454126/148‑1 (MQ=255)
                                         ccccccccaaaaaTCAGTTCATCCGCTTCAAAATACTAACCCCAAAA‑CACCAGAAAAATGAGCGCCTTGCCCACAATCACCCCCGCGCGCGCCACAAAAACACCTATCCCCCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTc                                                                                                      <  2:919598/141‑1 (MQ=255)
                                           ccTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCtg                                                                                                    >  1:861966/1‑151 (MQ=255)
                                           ccTCCCCAAAATCGGGTCATCCGCTTCGAGAT‑CCTAACCCCAAACCCCCAGCAAAAGGGGCGGCTGGCCGACGACCACCCCCGGGGGCGCCACCAGAATCCCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCtg                                                                                                    <  2:757226/151‑1 (MQ=255)
                                                   aaaTCAGGTCATACGCTTCGAGATACTTAACGCCAAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTgcg                                                                                            >  1:830598/1‑151 (MQ=255)
                                                   aaaTAAGGTCAACCGCTTCGAGATACTAAACCCAAAA‑CCCCAGCAAAATGGGCGGCTTCCCGACGAAAAACCCCGCGAGCCCCACAAAAACACCTACCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTggtggt                                                                                               <  2:861818/148‑1 (MQ=255)
                                                             atccGCTTCAAAATACTTACCGCCAAC‑CACCAGAGAAAGGGGCGCCTCCCCGACGACCCACACCGCGACCGCCACCAAACTCCCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGa                                                                                        <  1:346729/142‑1 (MQ=255)
                                                             aTACGCTTCGAGATACTTAACGCCAAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTaa                                                                                   >  2:875078/1‑150 (MQ=255)
                                                                    tCGAGATACTTAACGCCAAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATc                                                                            >  2:52606/1‑150 (MQ=255)
                                                                      gagaTACTTAACGCCAAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCAcc                                                                         >  2:588363/1‑151 (MQ=255)
                                                                           aCTTAACGCCAAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTgg                                                                    >  1:424050/1‑151 (MQ=255)
                                                                               accccaaaa‑CACCAGCAAAAGGAGCGGCTGGCCGACCAACAACCCCGGGCCCGCCACAAAAACCCCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTaa                                                                <  2:555596/145‑1 (MQ=255)
                                                                                 ccccAAC‑CCCCAGAAAAAGGGGCGGCTGGCCGCCAACCAACCCGCCGCGCCCCCACAAAATCCCTATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTg                                                                     <  1:5824/142‑1 (MQ=255)
                                                                                       accacCAGAGAAAGGAGCGGCTCCCCGACGACCAACCCCGCGAGCCCCACAAAAACACCGACCCCCCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACt                                                           <  1:1012944/147‑1 (MQ=255)
                                                                                           ccAGCAAAAGGGGCGGCTGGCCGACAATCACCCCCGCGGCCCCCACAAAAATCCCTATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTa                                                       <  2:588652/149‑1 (MQ=255)
                                                                                                gAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGc                                                 >  1:322742/1‑150 (MQ=255)
                                                                                                 aaaTGGGCGGCTTGCCGACGACCAACCCCGCGAGCGCCACAAGAATCCCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAg                                                  <  2:1035288/148‑1 (MQ=255)
                                                                                                 aaaTGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCa                                                >  2:555526/1‑150 (MQ=255)
                                                                                                     gAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAAccc                                           >  1:644875/1‑151 (MQ=255)
                                                                                                                cgacgaTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTcc                                >  2:334384/1‑151 (MQ=255)
                                                                                                                 gacgaTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcg                               >  1:319893/1‑151 (MQ=255)
                                                                                                                    caccacccccGCGCGCCCCACCAAAATCCCGTCCCCCCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcggcg                            <  1:607677/143‑1 (MQ=255)
                                                                                                                       cacccccGCGACCCCCACCAAAACCCCGTTCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAAccc                                           <  2:27780/128‑1 (MQ=255)
                                                                                                                       cAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGgcggcg                         >  1:41940/1‑151 (MQ=255)
                                                                                                                          cccggcgACCCCCACCAAAATCCCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTc                       <  2:277012/146‑1 (MQ=255)
                                                                                                                             cgcggccccccccaAAATACCGATCCCCCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCg                      <  1:465573/137‑1 (MQ=255)
                                                                                                                              ggatcccccccaacaaTACCGTCCCCCCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCga                     <  2:860199/134‑1 (MQ=255)
                                                                                                                                    gcCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTc             >  2:1038240/1‑150 (MQ=255)
                                                                                                                                    cccacaaGAAACCCTATCCCTCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCa            <  2:670235/150‑1 (MQ=255)
                                                                                                                                      caccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATg          >  1:1017275/1‑151 (MQ=255)
                                                                                                                                                ccGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCCGGTGTTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGCCATGGCGAtttt  >  1:258551/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
AAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTT  >  NC_000913/3131292‑3131583

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑