Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 941,253 C→T 8.1% H99Y (CAC→TAC)  dmsA → dimethyl sulfoxide reductase, anaerobic, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913941,2530CT 8.1% 147.3 / 3.1 49H99Y (CAC→TAC) dmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base C (21/24);  new base T (2/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

CAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAGGCAAATTCGAGCGCATTAGCTGGGAAGAAGCCTACGACA  >  NC_000913/941110‑941394
                                                                                                                                               |                                                                                                                                             
caacaaAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAg                                                                                                                                          <  6:70825/149‑1 (MQ=255)
 aacaaAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAgg                                                                                                                                         >  5:204256/1‑149 (MQ=255)
  acaaAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGt                                                                                                                                        >  6:10181/1‑149 (MQ=255)
      aaTCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTcgc                                                                                                                                    <  1:231427/149‑1 (MQ=255)
      aaTCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTcgc                                                                                                                                    <  3:71098/149‑1 (MQ=255)
                          tGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTc                                                                                                                >  2:45622/1‑149 (MQ=255)
                           ggAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTc                                                                                                                >  6:29864/1‑148 (MQ=255)
                                gcCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAAAATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCg                                                                                                          >  6:168060/1‑149 (MQ=255)
                                      aCAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCg                                                                                                    >  1:231088/1‑149 (MQ=255)
                                       gaGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTCGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGGCGCCgt                                                                                                   <  6:168865/148‑1 (MQ=255)
                                        aGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCgtg                                                                                                  >  3:192378/1‑149 (MQ=255)
                                         gTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCgtgt                                                                                                 >  3:45929/1‑149 (MQ=255)
                                           tAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTc                                                                                                <  7:138844/148‑1 (MQ=255)
                                             aCTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTa                                                                                              <  5:79459/148‑1 (MQ=255)
                                             aCTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTa                                                                                              <  3:96656/148‑1 (MQ=255)
                                             aCTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTa                                                                                              <  7:62282/148‑1 (MQ=255)
                                             aCTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAAATACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGACGTGTCTa                                                                                              >  2:43835/1‑148 (MQ=255)
                                                 tgGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAg                                                                                                                                          >  7:45791/1‑100 (MQ=255)
                                                 tgGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAg                                                                                                                                          <  8:45791/100‑1 (MQ=255)
                                                   gTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCg                                                                                       >  2:222941/1‑149 (MQ=255)
                                                         gCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGc                                                                                 >  8:82212/1‑149 (MQ=255)
                                                              ccGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGa                                                                             >  7:30253/1‑148 (MQ=255)
                                                               cGCTACGTCTGCACGTCGTGGTCGGGGAAATCATATATGGCGAAACTGACTATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCTTGGGCGTTCCATGCGTCGCCGTGGCTACAATCCGGACCGCCTGaaa                                                                           <  5:10181/149‑1 (MQ=255)
                                                                gCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGTACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGaaa                                                                           >  8:67184/1‑148 (MQ=255)
                                                                   aCGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATc                                                                       >  5:28433/1‑149 (MQ=255)
                                                                    cGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAatat                                                                        <  1:222941/147‑1 (MQ=255)
                                                                      tATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGa                                                                    <  5:95526/149‑1 (MQ=255)
                                                                      tATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGa                                                                    >  4:12354/1‑149 (MQ=255)
                                                                        tGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGTTGACAATTACGTCTGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGAt                                                                   <  2:173480/148‑1 (MQ=255)
                                                                        tGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATg                                                                  <  2:231088/149‑1 (MQ=255)
                                                                            cgtcgtGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGTACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAAc                                                              <  7:67184/149‑1 (MQ=255)
                                                                              tcgtGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGa                                                            <  1:114139/149‑1 (MQ=255)
                                                                                                 atatGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAg                                         <  4:45929/149‑1 (MQ=255)
                                                                                                     gTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCCCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAAAACGCACCGCCTGAAATATCCGATGAAACGAGTCGTGGCGCCCGGTCAAggaa                                      >  6:2440/1‑146 (MQ=255)
                                                                                                      tCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGCGTCTACAATCCGGACCGCCTGAAATATACGATGAAACGAGTCGGGGCGCGCGGTGAAGGCaaa                                    >  5:180081/1‑149 (MQ=255)
                                                                                                        gAAACGGACAATACCGGCGATCACAAGTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATTCGTCGCCGTGTCTACAATCCGGAACGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAGGCAAAtt                                  >  6:202305/1‑149 (MQ=255)
                                                                                                         aaaCGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAGGCAAATTc                                 <  5:168060/149‑1 (MQ=255)
                                                                                                           aCGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAGGCAAATTCGa                               <  4:192378/149‑1 (MQ=255)
                                                                                                                cAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAGGCAAATTCGAGCGCa                          <  1:45622/149‑1 (MQ=255)
                                                                                                                   tACCGGCGATGACAATTACGACGGCCTGTACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAGGCAAATTCGAGCGCAt                         >  7:155697/1‑147 (MQ=255)
                                                                                                                   tACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGGGGGCGTTCCATGCGTCGCCGTGGCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAGGCAAATTCGAGCGCATTa                       <  5:29864/149‑1 (MQ=255)
                                                                                                                        gCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAGGCAAATTCGAGCGCATTAGCTgg                  <  8:30253/149‑1 (MQ=255)
                                                                                                                               cAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACACTCCGGACCGCCTGAAATATCCGATGACACGCGTCGGGGCGCGCGGTCAAGGCAAATTCGAGCGCATTAGCTGGGaagaag           >  5:68932/1‑149 (MQ=255)
                                                                                                                                  ttACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAatat                                                                        <  6:198212/85‑1 (MQ=255)
                                                                                                                                  ttACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAatat                                                                        >  5:198212/1‑85 (MQ=255)
                                                                                                                                      gacgGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAGGCAAATTCGAGCGCATTAGCTGGGAAGAAGCCTACGa    <  6:28433/149‑1 (MQ=255)
                                                                                                                                       acgGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGCCCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGGGAAGACAAATTCGAGCGCATTAGCTGGGAAGAAGCcgacga    >  3:17506/1‑148 (MQ=255)
                                                                                                                                       acgGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAGGCAAATTCGAGCGCATTAGCTGGGAAGAAGCCTACGAc   <  6:204256/149‑1 (MQ=255)
                                                                                                                                        cgGCCTGTACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAGGCAAATTCGAGCGCATTAGCTGGGAAGAAGCCTACGAca  <  8:155697/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                             
CAACAAAATCAGACGAAAAGGTTATCTGGAGCGCCTGTACAGTTAACTGTGGTAGTCGCTGCCCGCTACGTATGCACGTCGTGGACGGTGAAATCAAATATGTCGAAACGGACAATACCGGCGATGACAATTACGACGGCCTGCACCAGGTTCGCGCCTGCCTGCGTGGGCGTTCCATGCGTCGCCGTGTCTACAATCCGGACCGCCTGAAATATCCGATGAAACGAGTCGGGGCGCGCGGTGAAGGCAAATTCGAGCGCATTAGCTGGGAAGAAGCCTACGACA  >  NC_000913/941110‑941394

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