Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913941,7130TG48.8% 53.6 / 10.5 41V252G (GTG→GGG) dmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (2/19);  new base G (20/0);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 9.40e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAG‑TGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941569‑941862
                                                                                                                                                |                                                                                                                                                       
actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACtt                                                                                                                                                    <  1:825874/150‑1 (MQ=255)
  tcctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtt                                                                                                                                                    >  1:513586/1‑148 (MQ=255)
     tccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTAc                                                                                                                                                  >  1:439393/1‑147 (MQ=255)
       cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCtt                                                                                                                                            <  1:704070/151‑1 (MQ=255)
             aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAg                                                                                                                                      <  2:94985/151‑1 (MQ=255)
               ttGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACa                                                                                                                                       >  2:500667/1‑148 (MQ=255)
                 gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCAc                                                                                                                                  <  1:13456/151‑1 (MQ=255)
                 gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCAc                                                                                                                                  <  1:501839/151‑1 (MQ=255)
                     aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTAATATTTTTAACACGCACGcc                                                                                                                               >  1:955504/1‑150 (MQ=255)
                        gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGcc                                                                                                                               >  1:753672/1‑147 (MQ=255)
                            gAACTATACCTACGGCTGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATc                                                                                                                       <  2:145001/150‑1 (MQ=255)
                                   aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAAAGGCACGCCAG‑AAATCTAATGccc                                                                                                               >  1:335621/1‑151 (MQ=255)
                                       aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCAt                                                                                                           <  1:901691/151‑1 (MQ=255)
                                       aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACCTATTATCTTTAAAAGGCCCGACCG‑AAAAATAATTccc                                                                                                               >  1:603650/1‑147 (MQ=255)
                                         ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGa                                                                                                         <  1:808158/151‑1 (MQ=255)
                                               tGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAAGCACGCCAG‑AAATCTAAAGccc                                                                                                               >  2:284591/1‑139 (MQ=255)
                                                 ggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCAACCCCG‑AAATCTAATGCCCGCCTGatcat                                                                                                     >  2:763855/1‑147 (MQ=255)
                                                    aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGAt                                                                                              <  2:374076/151‑1 (MQ=255)
                                                           aaCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTTAAAAAGGACGCCAG‑AAAACTAATGCCCGCATGATCATCATCGATcc                                                                                            >  1:198442/1‑146 (MQ=255)
                                                                cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAcc                                                                                  <  1:622258/151‑1 (MQ=255)
                                                                        atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCgg                                                                          <  1:987448/151‑1 (MQ=255)
                                                                          atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAAGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAACGACACCGGTg                                                                        >  1:843237/1‑151 (MQ=255)
                                                                             gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCg                                                                     <  1:457990/151‑1 (MQ=255)
                                                                                aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgc                                                                  <  2:800412/151‑1 (MQ=255)
                                                                                aaCAGTAAGCTGGTAGTGCTGTTTAGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATTTTTAAAAGGCAACCCAG‑AAAACTAAAGCCCGCCTGATTCTCATTGATCCCCGCTATACCGAAAACGGGGCCggg                                                                   >  2:951284/1‑150 (MQ=255)
                                                                                 aCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAAGCACGCCAG‑AAAACTAATGCCCGCATGATCATTATCGATCCGCGCTATACCGACACCGGGGCCggg                                                                   >  1:961014/1‑149 (MQ=255)
                                                                                  cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgcgc                                                                <  2:524934/151‑1 (MQ=255)
                                                                                  cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCggg                                                                   >  1:32579/1‑148 (MQ=255)
                                                                                   aGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGcgcg                                                               <  2:713778/151‑1 (MQ=255)
                                                                                          tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATTTTTAACAGGCACGCCAG‑AAAACTAATGCCCGCATGATCATCATCGAACCGCGCTATACCGACACCCGGGCCGGGGGCGGAGATg                                                         >  2:761110/1‑150 (MQ=255)
                                                                                           ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATg                                                         <  2:899471/149‑1 (MQ=255)
                                                                                            gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCCCGCTATACCGACAACGGTGCCGGGcgcgga                                                             >  1:528891/1‑142 (MQ=255)
                                                                                             tAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTATATCTTGAACAGGCACGCCAG‑AAAACTAATACCCGCCTGATCATCATTCATTCGCGCTTTACCGGCACCGGGGGCGGGGGCGGAGAAGAG‑Tg                                                    >  2:203098/1‑151 (MQ=255)
                                                                                                          taataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTATTATCTTGAACAGGCACGCCAG‑AAAACTAAAGACCGCATGATGATCATCGATCCGCGCCATACCGACAACGGGGTCGGGGGGGGAGAGGAG‑TGGGTCCCTATTc                                         >  2:327889/1‑149 (MQ=255)
                                                                                                            ataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTGAATAGGGACGCCAGAAAACCTAAT‑CCACCAAGCACCTTATTTTTGCCCGCTATACCGTCACCCGGGCCGAGCGGGAAGATGAGATGGTTCCC‑ATTCGt                                       >  2:582090/1‑149 (MQ=255)
                                                                                                            ataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCCG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAG‑TGGGTCCCTATTCGTcc                                     >  1:320970/1‑151 (MQ=255)
                                                                                                                     gCGAAACGTGAATGAGTGGCGGTGGGGTGGCTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGTTCCGCGCTCTACCGACACCGGTGCCGGGCGCGAAGATGAG‑TGGAGCCCTATTCGTCCGGGAACAGa                            <  2:580646/151‑1 (MQ=255)
                                                                                                                             cgAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAG‑TGGATCCCTATTCGTCCGGGAACAGATGCCGc                      >  1:444601/1‑149 (MQ=255)
                                                                                                                                  gAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAG‑TGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGtt               <  2:657210/151‑1 (MQ=255)
                                                                                                                                      ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAG‑TGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg           >  2:408142/1‑151 (MQ=255)
                                                                                                                                        cGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAG‑TGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt         <  2:32578/151‑1 (MQ=255)
                                                                                                                                                tGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAG‑TGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGt  >  1:637874/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                       
ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAG‑AAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAG‑TGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941569‑941862

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