Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,062 C→T 51.8% G840G (GGC→GGT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,0620CT51.8% ‑2.2 / 72.8 81G840G (GGC→GGTnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (18/21);  new base T (26/16);  total (44/37)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

AAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAA  >  minE/833998‑834122
                                                                |                                                            
aaGCGAAAAGCGACTGTGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCAt                                                           <  1:2179751/68‑1 (MQ=255)
aaGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCTAAGAAATTCTCCGAAGTGTGCGTTGGTCAt                                                           <  1:373424/68‑1 (MQ=255)
aaGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCAt                                                           <  1:1972254/68‑1 (MQ=255)
aaGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCAt                                                           <  1:2002614/68‑1 (MQ=255)
aaGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCAt                                                           <  1:1636959/68‑1 (MQ=255)
aaGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCAt                                                           <  1:2572561/68‑1 (MQ=255)
aaGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCAt                                                           <  1:3048937/68‑1 (MQ=255)
aaGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCAt                                                           <  1:2517925/68‑1 (MQ=255)
aaGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCAt                                                           <  1:144576/68‑1 (MQ=255)
aaGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCAt                                                           <  1:3411540/68‑1 (MQ=255)
aaGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCAt                                                           <  1:442802/68‑1 (MQ=255)
aaGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCAt                                                           <  1:484687/68‑1 (MQ=255)
aaGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCAt                                                           <  1:2459416/68‑1 (MQ=255)
aaGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCAt                                                           <  1:827457/68‑1 (MQ=255)
aaGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCAt                                                           <  1:2323190/68‑1 (MQ=255)
aaGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAACAAATTCTCCGAAGTGTGCGTTGGTCAt                                                           <  1:2515061/68‑1 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATc                                                          >  1:3280213/1‑48 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAAGTGACATc                                      >  1:3391382/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:2616834/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:2341592/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:2800266/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:2908151/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:2880597/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:289969/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:2900156/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:938117/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:3155589/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:519947/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:808160/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:90610/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:934978/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:1057876/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:2088851/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:1173125/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:1196444/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:1371294/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:1462947/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:1897517/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:1973169/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:204195/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTTGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:2045779/1‑68 (MQ=255)
                     tCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTAGGTCATCTGGGTAAAGAAACTGACATc                                      >  1:1064542/1‑68 (MQ=255)
                      cTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGt                                    <  1:2550409/69‑1 (MQ=255)
                      cTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGt                                    <  1:1068048/69‑1 (MQ=255)
                      cTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGt                                    <  1:671584/69‑1 (MQ=255)
                      cTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGt                                    <  1:667372/69‑1 (MQ=255)
                      cTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGt                                    <  1:1740803/69‑1 (MQ=255)
                      cTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGGCGGCCATCTGGGTAAAGAAACCGACATCGt                                    <  1:3277857/69‑1 (MQ=255)
                                   gcgAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTa                         >  1:2111572/1‑67 (MQ=255)
                                   gcgAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTa                         >  1:792549/1‑67 (MQ=255)
                                   gcgAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTa                         >  1:599474/1‑67 (MQ=255)
                                   gcgAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTa                         >  1:1314562/1‑67 (MQ=255)
                                   gcgAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTa                         >  1:478556/1‑67 (MQ=255)
                                   gcgAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTa                         >  1:2060349/1‑67 (MQ=255)
                                   gcgAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTa                         >  1:1490452/1‑67 (MQ=255)
                                   gcgAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTa                         >  1:1581622/1‑67 (MQ=255)
                                   gcgAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTa                         >  1:2038018/1‑67 (MQ=255)
                                   gcgAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTa                         >  1:1596059/1‑67 (MQ=255)
                                   gcgAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTa                         >  1:2710378/1‑67 (MQ=255)
                                   gcgAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTa                         >  1:2457491/1‑67 (MQ=255)
                                   gcgAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTa                         >  1:2855102/1‑67 (MQ=255)
                                   gcgAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTa                         >  1:2857191/1‑67 (MQ=255)
                                   gcgAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTa                         >  1:2447417/1‑67 (MQ=255)
                                      aagaaATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAg                    <  1:204585/69‑1 (MQ=255)
                                      aagaaATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAg                    <  1:2965654/69‑1 (MQ=255)
                                      aagaaATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAg                    <  1:1333334/69‑1 (MQ=255)
                                      aagaaATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAg                    <  1:1578441/69‑1 (MQ=255)
                                                        tgCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTg    >  1:523068/1‑67 (MQ=255)
                                                        tgCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGaa  >  1:2262736/1‑69 (MQ=255)
                                                        tgCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGaa  >  1:2202709/1‑69 (MQ=255)
                                                          cGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGcc                           <  1:45671/42‑1 (MQ=255)
                                                          cGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGcc                           <  1:1349511/42‑1 (MQ=255)
                                                          cGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGcc                           <  1:2087097/42‑1 (MQ=255)
                                                          cGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGcc                           <  1:1223076/42‑1 (MQ=255)
                                                          cGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGcc                           <  1:3407953/42‑1 (MQ=255)
                                                          cGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGcc                           <  1:639227/42‑1 (MQ=255)
                                                          cGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGcc                           <  1:3406009/42‑1 (MQ=255)
                                                          cGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGcc                           <  1:1620945/42‑1 (MQ=255)
                                                          cGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGcc                           <  1:294348/42‑1 (MQ=255)
                                                          cGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGcc                           <  1:111287/42‑1 (MQ=255)
                                                          cGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGcc                           <  1:2055849/42‑1 (MQ=255)
                                                                |                                                            
AAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAA  >  minE/833998‑834122

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: