Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913941,7130TG42.9% 42.9 / 12.0 35V252G (GTG→GGG) dmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (5/15);  new base G (15/0);  total (20/15)
Fisher's exact test for biased strand distribution p-value = 4.87e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTC  >  NC_000913/941575‑941856
                                                                                                                                          |                                                                                                                                               
ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATc                                                                                                                                      <  2:707682/150‑1 (MQ=255)
       aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAg                                                                                                                              <  2:393445/151‑1 (MQ=255)
       aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTaa                                                                                                                                 >  2:681194/1‑148 (MQ=255)
         ttGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACaag                                                                                                                             >  2:479471/1‑148 (MQ=255)
                  gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGat                                                                                                                   <  1:286248/151‑2 (MQ=255)
                      gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCt                                                                                                               <  1:320975/151‑1 (MQ=255)
                      gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCt                                                                                                               <  2:351656/151‑1 (MQ=255)
                             aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGGACGCCCGAAATATAATGccc                                                                                                        >  1:798553/1‑151 (MQ=255)
                                tACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTTACATTTTTAAAAAGCACGCCAGaaaaattatgtcc                                                                                                        >  2:503811/1‑138 (MQ=255)
                                     cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTTTTATCTTTAACAGGCACGCCCGAAATCTAAAGccc                                                                                                        >  1:459530/1‑143 (MQ=255)
                                            gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCg                                                                                         <  1:559173/151‑1 (MQ=255)
                                                  ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcg                                                                                    <  1:214112/150‑1 (MQ=255)
                                                       cAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTTAACAGGCACGCCCGAAATCTAATGCCCGCCTGATCATCATCGATCcgc                                                                                   >  1:161259/1‑146 (MQ=255)
                                                          cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGTTCCTATTATTTTGTAGAGGCACGCCAGCAATCTAATTGCCGCGCAGTCATCACCGATCCGCGCcataa                                                                            >  2:693560/1‑149 (MQ=255)
                                                                 gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTTAAAAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGGTCCGCGCTATACCGacac                                                                      >  2:623693/1‑149 (MQ=255)
                                                                  atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGGAATATTAACCTGAAAAAGAACGCAAGAAAAGTGAAGGCCCCATGAACCTCATCCATCCCCGCTGTAccccgaaacg                                                                   >  1:4347/1‑143 (MQ=255)
                                                                    atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGt                                                                  <  1:680360/150‑1 (MQ=255)
                                                                    atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCGGTg                                                                 >  2:103327/1‑151 (MQ=255)
                                                                      cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGACTTACTATCTTGAAAAAGCACGCCCGAAATCTAAAGCCCGCATGGTCATCATCGATCCGCGCTATACCGACACCGGTg                                                                 >  2:645667/1‑149 (MQ=255)
                                                                       gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTTCTATCTTTAACAAGCACGCCCGAAAACTAATGCCCGCATGATCATCCTCGATCCGCGCTATACCGACACCGGTg                                                                 >  1:800379/1‑148 (MQ=255)
                                                                        aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCgg                                                             <  2:571778/151‑1 (MQ=255)
                                                                                        aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAggg                                              >  2:421880/1‑148 (MQ=255)
                                                                                         gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGa                                            <  2:511192/151‑1 (MQ=255)
                                                                                             tGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGTGACTTACTATTTTGAAAAGGCACGCCCGAAATCTAATGTGCGCCTGATCATCATCGATGCGGGCTATATCCACACCCGGGCCGGGGGCGTAGATGAGTGGGTccc                                        >  2:278616/1‑151 (MQ=255)
                                                                                                ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCCGGGCCGGGCGCGGAGATGAGTGGATccc                                        >  1:72979/1‑148 (MQ=255)
                                                                                                   gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGACGTTCTATCTTGAACAGGCACGCCCGAAAACTAATGGCCCCATGATCATTATCCGTCCGCGCTATAACGACCCCCGGGGCGGGCGCGGAGAGGAGTGGGTCCCTAtt                                    >  2:476495/1‑149 (MQ=255)
                                                                                                     aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGt                                 <  2:432303/150‑1 (MQ=255)
                                                                                                         aCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCCCACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCggg                            >  1:474863/1‑151 (MQ=255)
                                                                                                              ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAAc                         <  2:244796/149‑1 (MQ=255)
                                                                                                              ggCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAAGCACGCCAGAAAACTAATGCCCGCCTGAACATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGGGGATCCCTATTCCTCCggg                            >  1:61635/1‑146 (MQ=255)
                                                                                                               gCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACCCCAGAAATCTAAAGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGGAGATGAGGGGGACCCCATTTGTCCgggga                          >  2:528461/1‑145 (MQ=255)
                                                                                                                   aaCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGcc                  <  2:444333/151‑1 (MQ=255)
                                                                                                                              gTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTa        <  2:226741/150‑1 (MQ=255)
                                                                                                                                ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAAc      <  2:200947/150‑1 (MQ=255)
                                                                                                                                   ggTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTc  >  1:420547/1‑151 (MQ=255)
                                                                                                                                          |                                                                                                                                               
CCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTC  >  NC_000913/941575‑941856

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