Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,107 T→C 57.7% H855H (CAT→CAC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,1070TC57.7% 17.6 / 42.2 59H855H (CAT→CACnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (15/10);  new base C (20/14);  total (35/24)
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.75e-01

TGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGA  >  minE/834054‑834174
                                                     |                                                                   
tgCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTg                                                        >  1:523068/1‑67 (MQ=255)
tgCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGaa                                                      >  1:2262736/1‑69 (MQ=255)
tgCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGaa                                                      >  1:2202709/1‑69 (MQ=255)
                   aaGAAACTGACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTg                                                   >  1:1035166/1‑53 (MQ=255)
                   aaGAAACTGACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTg                                                   >  1:2621894/1‑53 (MQ=255)
                   aaGAAACTGACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTg                                                   >  1:2447446/1‑53 (MQ=255)
                   aaGAAACTGACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTg                                                   >  1:2177150/1‑53 (MQ=255)
                   aaGAAACTGACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTg                                                   >  1:2841081/1‑53 (MQ=255)
                   aaGAAACTGACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTg                                                   >  1:2936416/1‑53 (MQ=255)
                   aaGAAACTGACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTg                                                   >  1:194809/1‑53 (MQ=255)
                   aaGAAACTGACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTg                                                   >  1:1941700/1‑53 (MQ=255)
                   aaGAAACTGACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTg                                                   >  1:3138525/1‑53 (MQ=255)
                   aaGAAACTGACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTg                                                   >  1:1836885/1‑53 (MQ=255)
                   aaGAAACTGACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTg                                                   >  1:1816368/1‑53 (MQ=255)
                   aaGAAACTGACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTg                                                   >  1:3194174/1‑53 (MQ=255)
                   aaGAAACTGACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTg                                                   >  1:1664688/1‑53 (MQ=255)
                   aaGAAACTGACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTg                                                   >  1:3229071/1‑53 (MQ=255)
                   aaGAAACTGACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTg                                                   >  1:1570834/1‑53 (MQ=255)
                   aaGAAACTGACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTg                                                   >  1:1080198/1‑53 (MQ=255)
                   aaGAAACTGACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTg                                                   >  1:908802/1‑53 (MQ=255)
                   aaGAAACTGACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTg                                                   >  1:591271/1‑53 (MQ=255)
                   aaGAAACTGACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTg                                                   >  1:837319/1‑53 (MQ=255)
                   aaGAAACTGACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTg                                                   >  1:1419337/1‑53 (MQ=255)
                      aaaCCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGa                                <  1:978939/69‑1 (MQ=255)
                      aaaCCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGa                                <  1:278906/69‑1 (MQ=255)
                       aaCCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGAtgt                                  >  1:530613/1‑66 (MQ=255)
                          cGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGa                            <  1:308508/69‑1 (MQ=255)
                          cGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGa                            <  1:2193603/69‑1 (MQ=255)
                           gACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTGGCGCAACCGCTGGATGTGAAAGAc                           <  1:1612252/69‑1 (MQ=255)
                           gACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTGGCGCAACCGCTGGATGTGAAAGAc                           <  1:2405481/69‑1 (MQ=255)
                           gACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTGGCGCAACCGCTGGATGTGAAAGAc                           <  1:2710147/69‑1 (MQ=255)
                           gACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTGGCGCAACCGCTGGATGTGAAAGAc                           <  1:2276909/69‑1 (MQ=255)
                           gACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTGGCGCAACCGCTGGATGTGAAAGAc                           <  1:1520932/69‑1 (MQ=255)
                           gACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTGGCGCAACCGCTGGATGTGAAAGAc                           <  1:317010/69‑1 (MQ=255)
                           gACATCGTCACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTGGCGCAACCGCTGGATGTGAAAGAc                           <  1:1833031/69‑1 (MQ=255)
                           gACATCGTCACGCTGCCTAACCACCACGACTCTGCCGCTGAACTGGCGCAACCGCTGGATGTGAAAGAc                           <  1:1899343/69‑1 (MQ=255)
                                   cACGCTGCCTATCCAGCACGACTCTGCCGCTGAACTGGCGCAACCGCTGGATGTGAAAGAc                           <  1:1549700/61‑1 (MQ=255)
                                        tGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCg              <  1:383702/69‑1 (MQ=255)
                                        tGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCg              <  1:270032/69‑1 (MQ=255)
                                        tGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCg              <  1:105807/69‑1 (MQ=255)
                                        tGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCg              <  1:1932621/69‑1 (MQ=255)
                                             aTCCAGCACGACTCTGCCGCTGAACTGGCGCAACCGCTGGATGTGAATGAc                           <  1:902562/51‑1 (MQ=38)
                                             aTCCAGCACGACTCTGCCGCTGAACTGGCGCAACCGCTGGATGTGAAAGAc                           <  1:3328005/51‑1 (MQ=39)
                                             aTCCACCACGACTCTGCCGCTGAACTGGCGCAACCGCTGGATGTGAAAGAc                           <  1:262588/51‑1 (MQ=38)
                                             aTCCACCACGACTCTGCCGCTGAACTGGCGCAACCGCTGGATGTGAAAGAc                           <  1:3148449/51‑1 (MQ=38)
                                               ccAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGc               <  1:1548990/61‑1 (MQ=255)
                                               ccAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGc               >  1:3125758/1‑61 (MQ=255)
                                               ccAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGc               <  1:2003030/61‑1 (MQ=255)
                                                cAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGAc                           >  1:2445180/1‑48 (MQ=255)
                                                cAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGAc      >  1:860090/1‑69 (MQ=255)
                                                cAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGAc      >  1:134236/1‑69 (MQ=255)
                                                cAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGAc      >  1:3139990/1‑69 (MQ=255)
                                                cAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGAc      >  1:672852/1‑69 (MQ=255)
                                                cAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGAc      >  1:1495621/1‑69 (MQ=255)
                                                cAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGAc      >  1:2301947/1‑69 (MQ=255)
                                                cAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGAc      >  1:352087/1‑69 (MQ=255)
                                                cAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGAc      >  1:418746/1‑69 (MQ=255)
                                                  gCACGACTCTGCCGCTGAACTGGCGCAACCGCTGGATGTGAAAGAc                           <  1:2173587/46‑1 (MQ=38)
                                                    aTGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGa  >  1:3012109/1‑69 (MQ=255)
                                                     tGACTCTGGCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCg              >  1:1426121/1‑56 (MQ=255)
                                                     tGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCg              >  1:1293495/1‑56 (MQ=255)
                                                     tGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCg              >  1:1497909/1‑56 (MQ=255)
                                                     tGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCg              >  1:526202/1‑56 (MQ=255)
                                                     tGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCg              >  1:1526702/1‑56 (MQ=255)
                                                     tGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCg              >  1:1640857/1‑56 (MQ=255)
                                                     tGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCg              >  1:1731022/1‑56 (MQ=255)
                                                     tGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCg              >  1:1835201/1‑56 (MQ=255)
                                                     tGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCg              >  1:3083598/1‑56 (MQ=255)
                                                     tGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCg              >  1:3076164/1‑56 (MQ=255)
                                                     tGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCg              >  1:2015773/1‑56 (MQ=255)
                                                     tGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGCCg              >  1:2201807/1‑56 (MQ=255)
                                                     |                                                                   
TGCGTCGGCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGA  >  minE/834054‑834174

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: