Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,119 T→G 63.5% A279A (GCT→GCG narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE836,1190TG63.5% 33.6 / 27.6 41A279A (GCT→GCGnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (7/8);  new base G (6/20);  total (13/28)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

CCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCTTTACCAGCGTCAGCTGGACGTGTTCCTC  >  minE/836052‑836185
                                                                   |                                                                  
ccGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTa                                                                   <  1:3390407/69‑1 (MQ=255)
ccGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTa                                                                   <  1:2424121/69‑1 (MQ=255)
     tgttCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGCTGTACGATGCCGACGCGAtt                                                                 <  1:2871531/63‑1 (MQ=255)
     tgCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGaa                                                              <  1:1413445/69‑1 (MQ=255)
     tgCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGaa                                                              <  1:1417569/69‑1 (MQ=255)
     tgCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGaa                                                              <  1:1438522/69‑1 (MQ=255)
     tgCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGaa                                                              <  1:3364354/69‑1 (MQ=255)
        tCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGaa                                                              <  1:1148584/66‑1 (MQ=255)
                     tcggtcgTATCCGTTATCTTGGCGTGCTGCTGTACGATGCCGACGCGATTGAACGTGcggcaagta                                                 >  1:1500408/1‑64 (MQ=255)
                     tcggtcgTATCCGTTATCTTGGCGTGCTGCTGTACGATGCCGACGCGATTGAACGTGcggcaagta                                                 >  1:2654348/1‑64 (MQ=255)
                     tcggtcgTATCCGTTATCTTGGCGTGCTGCTGTACGATGCCGACGCGATTGAACGTGcggcaagta                                                 >  1:243255/1‑64 (MQ=255)
                     tcggtcgTATCCGTTATCTTGGCGTGCTGCTGTACGATGCCGACGCGATTGAACGTGcggcaagta                                                 >  1:1218498/1‑64 (MQ=255)
                     tcggtcgTATCCGTTATCTTGGCGTGCTGCTGTACGATGCCGACGCGATTGAACGTGcggcaagta                                                 >  1:1151431/1‑64 (MQ=255)
                     tcggtcgTATCCGTTATCTTGGCGTGCTGCTGTACGATGCCGACGCGATTGAACGTGcggcaagta                                                 >  1:3183048/1‑64 (MQ=255)
                                cGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACcgagaacga                                       >  1:778388/1‑65 (MQ=255)
                                    aTCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTg                                                          >  1:416477/1‑42 (MQ=255)
                                       ttGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCtt                             <  1:925299/68‑1 (MQ=255)
                                                 gctgtACGATGCCGACGCGATTGAACGTGCGGCAAGTACTGAGAACGAGAAAGATCTCTACCAGCGTCa                  <  1:238120/67‑1 (MQ=255)
                                                 gctgtACGATGCCGACGCGATTGAACGTGCGGCAAGTACCGAGAACGAGAAAGATCTCTACCAGCGTCa                  <  1:3351754/67‑1 (MQ=255)
                                                 gctgtACGATGCCGACGCGATTGAACGTGCGGCAAGTACCGAGAACGAGAAAGATCTCTACCAGCGTCa                  <  1:546074/67‑1 (MQ=255)
                                                 gctgtACGATGCCGACGCGATTGAACGTGCGGCAAGTACCGAGAACGAGAAAGATCTCTACCAGCGTCa                  <  1:344941/67‑1 (MQ=255)
                                                 gctgtACGATGCCGACGCGATTGAACGTGCGGCAAGTACCGAGAACGAGAAAGATCTCTACCAGCGTCa                  <  1:522549/67‑1 (MQ=255)
                                                 gctgtACGATGCCGACGCGATTGAACGTGCGGCAAGTACCGAGAACGAGAAAGATCTCTACCAGCGTCa                  <  1:2791760/67‑1 (MQ=255)
                                                 gctgtACGATGCCGACGCGATTGAACGTGCGGCAAGTACCGAGAACGAGAAAGATCTCTACCAGCGTCa                  <  1:650158/67‑1 (MQ=255)
                                                 gctgtACGATGCCGACGCGATTGAACGTGCGGCAAGTACCGAGAACGAGAAAGATCTCTACCAGCGTCa                  <  1:989059/67‑1 (MQ=255)
                                                 gctgtACGATGCCGACGCGATTGAACGTGCGGCAAGTACCGAGAACGAGAAAGATCTCTACCAGCGTCa                  <  1:3302751/67‑1 (MQ=255)
                                                 gctgtACGATGCCGACGCGATTGAACGTGCGGCAAGTACCGAGAACGAGAAAGATCTCTACCAGCGTCa                  <  1:1055057/67‑1 (MQ=255)
                                                 gctgtACGATGCCGACGCGATTGAACGTGCGGCAAGTACCGAGAACGAGAAAGATCTCTACCAGCGTCa                  <  1:2465319/67‑1 (MQ=255)
                                                 gctgtACGATGCCGACGCGATTGAACGTGCGGCAAGTACCGAGAACGAGAAAGATCTCTACCAGCGTCa                  <  1:1847760/67‑1 (MQ=255)
                                                 gctgtACGATGCCGACGCGATTGAACGTGCGGCAAGTACCGAGAACGAGAAAGATCTCTACCAGCGTCa                  <  1:16617/67‑1 (MQ=255)
                                                 gctgtACGATGCCGACGCGATTGAACGTGCGGCAAGTACCGAGAACGAGAAAGATCTCTACCAGCGTCa                  <  1:1483394/67‑1 (MQ=255)
                                                 gctgtACGATGCCGACGCGATTGAACGTGCGGCAAGTACCGAGAACGAGAAAGATCTCTACCAGCGTCa                  <  1:1380268/67‑1 (MQ=255)
                                                 gctgtACGATGCCGACGCGATTGAACGTGCGGCAAGTACCGAGAACGAGAAAGATCTCTACCAGCGTCa                  <  1:1081539/67‑1 (MQ=255)
                                                 gctgtACGATGCCGACGCGATTGAACGTGCGGCAAGTACCGAGAACGAGAAAGAGCTCTACCAGCGTCa                  <  1:471790/67‑1 (MQ=255)
                                                  ctgtACGATGCCGACGCGATTGAACGTGCGGCAAGTACCGAGAACGAGAAAGATCTCTACCAGCGTCa                  <  1:2974790/67‑1 (MQ=255)
                                                  ctgtACGATGCCGACGCGATTGAACGTGCGGCAAGTACCGAGAACGAGAAAGATCTCTACCAGCGTCa                  <  1:3382707/67‑1 (MQ=255)
                                                         aTGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCTTTACCAGCGTCAGCTGGACg          >  1:2892679/1‑69 (MQ=255)
                                                         aTGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCTTTACCAGCGTCAGCTGGACg          >  1:2343760/1‑69 (MQ=255)
                                                         aTGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCTTTACCAGCGTCAGCTGGACg          >  1:557210/1‑69 (MQ=255)
                                                         aTGCCGAAGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCTTTACCAGCGTCAGCTGGACg          >  1:1729600/1‑69 (MQ=255)
                                                                 gCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCTTTACCAGCGTCAGCTGGACGTGTTCCTc  >  1:3386408/1‑69 (MQ=255)
                                                                   |                                                                  
CCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCTTTACCAGCGTCAGCTGGACGTGTTCCTC  >  minE/836052‑836185

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