Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,664,769 (G)5→6 coding (264/615 nt) dmsD → redox enzyme maturation protein DmsD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,664,7641.G92.9% 42.7 / ‑2.7 14coding (259/615 nt)dmsDredox enzyme maturation protein DmsD
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (7/6);  total (8/6)
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.89e-01

GGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTG  >  NC_000913/1664626‑1664898
                                                                                                                                           |                                                                                                                                       
ggCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTg                                                                                                                                        >  1:400715/1‑140 (MQ=255)
     ggAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTgggggg                                                                                                                                  >  2:439109/1‑139 (MQ=255)
       aaaCTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTACCGTCTCCGCCGTggggggtc                                                                                                                                <  1:452617/140‑4 (MQ=255)
         aCTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTggggggtcgg                                                                                                                              <  2:193783/140‑6 (MQ=255)
                            aaGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGGTCGGTCTGGCTGGATCGCGAATc                                                                                                            >  2:370736/1‑140 (MQ=255)
                                 tCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGGTCGGTCTGGCTGGATcg                                                                                                                  >  1:280698/1‑129 (MQ=255)
                                 tCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGGTCGGTCTGGCTGGATcg                                                                                                                  <  2:280698/129‑1 (MQ=255)
                                         gccgcTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGc                                                                                               >  1:244526/1‑140 (MQ=255)
                                              tGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTc                                                                                          <  1:370736/140‑1 (MQ=255)
                                                          ttCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACt                                                                              >  1:326149/1‑140 (MQ=255)
                                                                    gtgtGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTgg                                                                    >  1:64701/1‑140 (MQ=255)
                                                                                         ggCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTACCGTCTCCGCCGTGGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAg                                               >  1:491279/1‑140 (MQ=255)
                                                                                              ggCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCAtt                                                  <  2:160560/132‑1 (MQ=255)
                                                                                              ggCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCAtt                                                  >  1:160560/1‑132 (MQ=255)
                                                                                              ggCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGGCTCCGCCGTGGGCGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGa                                          >  2:336923/1‑140 (MQ=255)
                                                                                                                 cccGTGGGCACTACCGTCTCCGCCGTGGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTTAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCg                       <  2:491279/140‑1 (MQ=255)
                                                                                                                    gTGGGCACTGCCGTCTCCGCCGTGGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGaa                    <  2:326149/140‑1 (MQ=255)
                                                                                                                                      gccgTGGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCActg  <  2:99887/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                       
GGCTGGGAAACTCAGTGGCCTCTACCAGAAGCGTCATTAGCGCCGCTGGTGACGGCTTTTCAGACACAGTGTGAAGAGACTCACGCCCAGGCCTGGCAGCGTTTGTTTGTCGGCCCGTGGGCACTGCCGTCTCCGCCGTGGGGGTCGGTCTGGCTGGATCGCGAATCTGTGCTGTTTGGCGATTCAACATTGGCACTTCGTCAGTGGATGCGCGAGAAAGGCATTCAGTTTGAAATGAAGCAAAACGAACCGGAAGATCATTTTGGTTCACTG  >  NC_000913/1664626‑1664898

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: