Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,285,540 G→A 43.7% S120S (TCG→TCA narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,285,5400GA43.7% 1.3 / 14.5 16S120S (TCG→TCAnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (6/3);  new base A (4/3);  total (10/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 = 1.00e+00

TCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAACCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATCGAAAAAGGGCC  >  W3110S.gb/1285470‑1285605
                                                                      |                                                                 
tCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAACCTGCATACCGCGCCGGAAGGCAGCAAATCg                                                                   >  1:2202960/1‑71 (MQ=255)
      attattACGAGCCGTTCGATTTTGACTATCAGAACCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCg                                                             <  1:1381371/71‑1 (MQ=255)
           taCGAGCCGTTCGATTTTGATTATCAGAACCTGCATACCGCGCCGGAAGGCAGCAAATCACAGCCGATTGc                                                        >  1:730696/1‑71 (MQ=255)
           taCGAGCCGTTCGATTTTGATTATCAGAACCTGCATACCGCGCCGGAAGGCAGCAAATCACAGCCGATTGc                                                        >  1:3354106/1‑71 (MQ=255)
           taCGAGCCGTTCGATTTTGATTATCAGAACCTGCATACCGCGCCGGAAGGCAGCAAATCACAGCCGATTGc                                                        >  1:2859320/1‑71 (MQ=255)
           taCGAGCCGTTCGATTTTGATTATCAGAACCTGCATACCGCGCCGGAAGGCAGCAAATCACAGCCGATTGc                                                        >  1:2348953/1‑71 (MQ=255)
           taCGAGCCGTTCGATTTTGACTATCAGAACCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGc                                                        >  1:1371823/1‑71 (MQ=255)
            aCGAGCCGTTCGATTTTGACTATCAGAACCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGcc                                                       >  1:1332936/1‑71 (MQ=255)
                          tttGACTATCAGAACCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCt                                         >  1:1188258/1‑71 (MQ=255)
                                            cATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGAt                       >  1:2083865/1‑71 (MQ=255)
                                                 cgcgCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTg                                        >  1:19861/1‑49 (MQ=255)
                                                    gcCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGaaaa               <  1:1623103/71‑1 (MQ=255)
                                                            ggagcaAATCACAGCCGATTGCCCGTCCGCGTTCACTGATTACCGGCGAACGGAt                       <  1:3722613/53‑1 (MQ=37)
                                                            gcagcaAATCACAGCCGATTGCCCGTCCGCGTTCACTGATTACCGGCGAACGGAt                       <  1:2101052/55‑1 (MQ=37)
                                                            gcagcaAATCACAGCCGATTGCCCGTCCGCGTTCACTGATTACCGGCGAACGGAt                       <  1:2669818/55‑1 (MQ=37)
                                                             cagcaAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATCGAAAAAg      <  1:1132446/71‑1 (MQ=255)
                                                                  aaTCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATCGAAAAAGGGcc  >  1:219576/1‑70 (MQ=255)
                                                                      |                                                                 
TCGACGATTATTACGAGCCGTTCGATTTTGACTATCAGAACCTGCATACCGCGCCGGAAGGCAGCAAATCGCAGCCGATTGCCCGTCCGCGTTCGCTGATTACCGGGGAACGGATGGCGAAAATCGAAAAAGGGCC  >  W3110S.gb/1285470‑1285605

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