Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913941,7130TG41.3% 75.6 / 11.8 46V252G (GTG→GGG) dmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (4/23);  new base G (19/0);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 4.26e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGC  >  NC_000913/941565‑941860
                                                                                                                                                    |                                                                                                                                                   
gACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGa                                                                                                                                                   <  1:76183/151‑1 (MQ=255)
  ctactCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAc                                                                                                                                                  >  1:1012155/1‑150 (MQ=255)
           cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCtt                                                                                                                                        <  1:337101/151‑1 (MQ=255)
             cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGa                                                                                                                                      <  1:313642/151‑1 (MQ=255)
                  aTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAgg                                                                                                                                 <  2:103337/151‑1 (MQ=255)
                  aTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACtttat                                                                                                                                         >  2:827384/1‑139 (MQ=255)
                     gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCa                                                                                                                               <  2:599575/150‑1 (MQ=255)
                     gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCAc                                                                                                                              <  1:842744/151‑1 (MQ=255)
                        gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGGAGTACTTTTCTGAAAAGGCACGcc                                                                                                                           >  1:687570/1‑151 (MQ=255)
                              ttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaaa                                                                                                                      <  1:534606/150‑1 (MQ=255)
                              ttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACGTATTAAATTGaaaacggccccccgaaa                                                                                                                      >  2:986991/1‑135 (MQ=255)
                               tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTTCTATCTTTTAAAAGGACGCCCGaaa                                                                                                                      >  1:535221/1‑149 (MQ=255)
                                     ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGc                                                                                                              <  2:952206/151‑1 (MQ=255)
                                       aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTTCTATCTTTAACAAGCACGCCCGAAAACTAATGccc                                                                                                            >  2:1125654/1‑151 (MQ=255)
                                            cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATg                                                                                                       <  2:633404/151‑1 (MQ=255)
                                             ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTTAACAAGCACGCCAGAAAACTAATGCCCGCCTg                                                                                                       >  2:90414/1‑150 (MQ=255)
                                              gcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTAAATTTAAAAGGCAACGCTCAAAACTAATTCCCGCag                                                                                                        >  2:614459/1‑147 (MQ=255)
                                               cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAAGCACGCCAGAAAACTAATGCCCGCATg                                                                                                       >  2:1125040/1‑148 (MQ=255)
                                                     ggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcatc                                                                                              <  2:578722/151‑1 (MQ=255)
                                                       cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGCACGCCAGAAATCTAATGCCCGCATGatcatcatc                                                                                              >  2:94290/1‑149 (MQ=255)
                                                       cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGAGGTAGTAAATTTAACCGAAACGCCCGAAATCTTATGGCCGCAAGataatgtact                                                                                             >  1:1132106/1‑145 (MQ=255)
                                                            ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcgc                                                                                       <  1:109697/151‑1 (MQ=255)
                                                            ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcgc                                                                                       <  1:440691/151‑1 (MQ=255)
                                                              cAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGGACGCCCGAAAACTAATGCCCGCATGATCATCATCGATcc                                                                                         >  1:36635/1‑147 (MQ=255)
                                                                aCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGGACGCCAGAAAACTAATAGGCGCCTGatcatcatttgcacgtgc                                                                                     >  1:118230/1‑139 (MQ=255)
                                                                        tCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTACTATTTTGAACAAACACGCCCGAAAACTAATGCCCGCCTGATCATCTCCGATCCCCGCTATAccca                                                                             >  1:968727/1‑147 (MQ=255)
                                                                           gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCg                                                                        <  1:18077/151‑1 (MQ=255)
                                                                           gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAGTTACTATTTTGAAAAGGGACCGCCGAAATCTAAAGCACGCATGAACATAAACGAACCCCCCTAAACCGACAccc                                                                        >  2:1076625/1‑150 (MQ=255)
                                                                                   aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGGACGCCAGAAAACTAAAGCCCCCCTGATCATCATCGATCCGCGCTATACCCACACCGGTGCCggg                                                                >  2:887608/1‑151 (MQ=255)
                                                                                      cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgcgc                                                             <  1:177544/151‑1 (MQ=255)
                                                                                       aGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGcgcg                                                            <  2:543521/151‑1 (MQ=255)
                                                                                               ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCCGTGCCGGGCGCGAAGa                                                        >  2:462129/1‑147 (MQ=255)
                                                                                                   gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGTGAAGTTCAATCTTTGAAAGGGACGGCAGAAAATTAAATCCCGCATGGACATCAACGTTACGCGCTTTACCGGCACCGGGGGCGGGcgcggaaatagggggg                                                >  2:75834/1‑139 (MQ=255)
                                                                                                         tttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTAATCTCTTGAACAAGCACGCCAGAAATCTAAAGGCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGAAGAGTGGATcc                                             >  2:969/1‑148 (MQ=255)
                                                                                                                      cTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCCTTAAAAGGCACGCCAGAAAACTTATGCCCGCATGATCATCATCGATTCGCGCCATACCGACCCCGGGGCCGGGCGCGAAGATGAGGGGATCCCCATTCGTCCgggga                              >  1:650900/1‑148 (MQ=255)
                                                                                                                        ggCGAAACGCGAATGAGTGGCGGTGGGGGGGATTACTATTTTGAAAAGGCACGCCAGAAAATTCATGCCCGCAAGATCATCATCGATCCGCGCTATAACGACACCCGGGGCGGGGGCGGAGATTAGTTGGACCCTTTTTGTCCGGGaa                              >  2:772901/1‑148 (MQ=255)
                                                                                                                           gAAACGCGAATGAGTGGCGGTGGGGTGGCTTACTATCTTGAAAAAGCACGCCAGAAATCTAATGCCCCCATGATCATCATCGATCCGCGCTATACCCACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGAt                         >  2:885419/1‑150 (MQ=255)
                                                                                                                            aaaCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCTCGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGc                       <  2:1027608/151‑1 (MQ=255)
                                                                                                                            aaaCGCGAATGAGTGGCGGTGGGGGGGACTACCATTTTTAAAAAGCACCCCAGAAATCTAATGCCCGCATGATCCTCATCGATCCCCGCTATACCGACACCCGGGCCGGGCGCGGAGATGAGTGGGTTCCTATTTGTCCCGGGACa                            >  1:610371/1‑146 (MQ=255)
                                                                                                                              aCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCg                     <  1:482990/151‑1 (MQ=255)
                                                                                                                                  gAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACt                 <  2:486290/151‑1 (MQ=255)
                                                                                                                                  gAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACt                 <  2:1330/151‑1 (MQ=255)
                                                                                                                                  gAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACt                 <  1:918447/151‑1 (MQ=255)
                                                                                                                                      gAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGtt             <  1:525492/151‑1 (MQ=255)
                                                                                                                                       aGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTa            <  2:808071/151‑1 (MQ=255)
                                                                                                                                                 gggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc  <  1:723981/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGC  >  NC_000913/941565‑941860

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: ‑