Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 935,197 C→G 49.5% T898T (ACC→ACG fdnG → formate dehydrogenase‑N, alpha subunit, nitrate‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE935,1970CG49.5% ‑2.3 / 119.1 95T898T (ACC→ACGfdnGformate dehydrogenase‑N, alpha subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base C (34/14);  new base G (27/20);  total (61/34)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

AACAAGACGCGCTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTGGACCAAGCACGCATTGCTCAACGCAATTGCTCAGCCGGAACAGTTTG  >  minE/935130‑935261
                                                                   |                                                                
aaCAAGACGCGCTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGCCCTATCGTCTGACCg                                                                 >  1:196294/1‑69 (MQ=255)
aaCAAGACGCGCTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCg                                                                 >  1:3155182/1‑69 (MQ=255)
aaCAAGACGCGCTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCg                                                                 >  1:2488139/1‑69 (MQ=255)
aaCAAGACGCGCTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCg                                                                 >  1:2610946/1‑69 (MQ=255)
aaCAAGACGCGCTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCg                                                                 >  1:2753512/1‑69 (MQ=255)
aaCAAGACGCGCTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCg                                                                 >  1:304368/1‑69 (MQ=255)
aaCAAGACGCGCTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCg                                                                 >  1:126644/1‑69 (MQ=255)
aaCAAGACGCGCTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCg                                                                 >  1:3293589/1‑69 (MQ=255)
aaCAAGACGCGCTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCg                                                                 >  1:39840/1‑69 (MQ=255)
aaCAAGACGCGCTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCg                                                                 >  1:583565/1‑69 (MQ=255)
aaCAAGACGCGCTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCg                                                                 >  1:1128779/1‑69 (MQ=255)
    aGACGCGCTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCa                                                             <  1:3004884/69‑1 (MQ=255)
    aGACGCGCTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCa                                                             <  1:2863952/69‑1 (MQ=255)
     gACGCGCTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCAc                                                            <  1:1159332/69‑1 (MQ=255)
     gACGCGCTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCAc                                                            <  1:3411538/69‑1 (MQ=255)
     gACGCGCTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCAc                                                            <  1:886017/69‑1 (MQ=255)
     gACGCGCTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCAc                                                            <  1:1927041/69‑1 (MQ=255)
        gcgcTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTc                                                         <  1:2690176/69‑1 (MQ=255)
        gcgcTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTc                                                         <  1:1018229/69‑1 (MQ=255)
                  tGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCg                                                                 >  1:1637233/1‑51 (MQ=255)
                  tGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTgg                                                >  1:2348531/1‑68 (MQ=255)
                  tGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTgg                                                >  1:3185968/1‑68 (MQ=255)
                  tGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTgg                                                >  1:2680829/1‑68 (MQ=255)
                  tGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTgg                                                >  1:705029/1‑68 (MQ=255)
                  tGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTgg                                                >  1:2166546/1‑68 (MQ=255)
                  tGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTgg                                                >  1:109469/1‑68 (MQ=255)
                  tGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTgg                                                >  1:1114005/1‑68 (MQ=255)
                  tGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTgg                                                >  1:1218512/1‑68 (MQ=255)
                  tGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTgg                                                >  1:1236116/1‑68 (MQ=255)
                  tGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTgg                                                >  1:1400206/1‑68 (MQ=255)
                  tGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCAAACCTgg                                                >  1:1481178/1‑68 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAGCACTTCCACACCTGGACCAAGc                                         <  1:1980794/69‑1 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAGCACTTCCACACCTGGACCAAGc                                         <  1:973874/69‑1 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAGCACTTCCACACCTGGACCAAGc                                         <  1:875122/69‑1 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAGCACTTCCACACCTGGACCAAGc                                         <  1:1119242/69‑1 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAGCACTTCCACACCTGGACCAAGc                                         <  1:4721/69‑1 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAGCACTTCCACACCTGGACCAAGc                                         <  1:366570/69‑1 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAGCACTTCCACACCTGGACCAAGc                                         <  1:2592750/69‑1 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAGCACTTCCACACCTGGACCAAGc                                         <  1:1006481/69‑1 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAGCACTTCCACACCTGGACCAAGc                                         <  1:2274636/69‑1 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAGCACTTCCACACCTGGACCAAGc                                         <  1:195942/69‑1 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAGCACTTCCACACCTGGACCAAGc                                         <  1:1984540/69‑1 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAGCACTTCCACACCTGGACCAAGc                                         <  1:217364/69‑1 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAGCACTTCCACACCTGGACCAAGc                                         <  1:1995174/69‑1 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAGCACTTCCACACCTGGACCAAGc                                         <  1:2081494/69‑1 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAGCACTTCCACACCTGGACCAAGc                                         <  1:2112532/69‑1 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAGCACTTCCACACCTGGACCAAGc                                         <  1:2134397/69‑1 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTGGACCAAg                                          >  1:2899731/1‑68 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTGGACCAAg                                          >  1:3291961/1‑68 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTGGACCAAg                                          >  1:1412137/1‑68 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTGGACCAAGc                                         >  1:2345637/1‑69 (MQ=255)
                        aaaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTGGACCAAGc                                         >  1:2955324/1‑69 (MQ=255)
                         aaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAGCACTTCCACACCTGGACCAAGc                                         <  1:1255190/68‑1 (MQ=255)
                         aaaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAGCACTTCCACACCTGGACCAAGc                                         <  1:2757531/68‑1 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGa                                                                >  1:3138583/1‑44 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGa                                                                >  1:2438041/1‑44 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:2052993/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:2022447/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:3320258/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:1204126/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:428951/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:2071533/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:471810/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:3267883/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:2017224/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:2015296/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:700092/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:1308513/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:2012865/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:2240210/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:2450336/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:2181632/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:2253847/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:2330751/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:2381577/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:1677333/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:2544055/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:1548760/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:2884466/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:1408334/1‑45 (MQ=255)
                          aaaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACGGAg                                                               >  1:1409815/1‑45 (MQ=255)
                           aaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTGGACCAAGCACg                                      <  1:1847764/69‑1 (MQ=255)
                           aaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTGGACCAAGCACg                                      <  1:806662/69‑1 (MQ=255)
                           aaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTGGACCAAGCACg                                      <  1:867353/69‑1 (MQ=255)
                           aaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTGGACCAAGCACg                                      <  1:2318257/69‑1 (MQ=255)
                           aaGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTGGACCAAGCACg                                      <  1:901830/69‑1 (MQ=255)
                            agagCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTGGACCAAGCACg                                      <  1:730825/68‑1 (MQ=255)
                                     cccGTATGTGGGTACGACCTATCGTCTGACGGAGCACTTCCACACCTGGACCAAGc                                         <  1:2429758/56‑1 (MQ=255)
                                     cccGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTGGACCAAGCACGCATTGCTCaa                            >  1:408476/1‑69 (MQ=255)
                                                      ccTATCGTCTGACGGAGCACTTCCACACCTGGACCAAGc                                         <  1:2990790/39‑1 (MQ=255)
                                                          tCGTCTGACCGAGCACTTCCACACCTGGACCAAGCACGCATTGCTCAACGCAATTGCTCAGCCGGAACa       >  1:1378535/1‑69 (MQ=255)
                                                          tCGTCTGACCGAGCACTTCCACACCTGGACCAAGCACGCATTGCTCAACGCAATTGCTCAGCCGGAACa       >  1:729654/1‑69 (MQ=255)
                                                          tCGTCTGACCGAGCACTTCCACACCTGGACCAAGCACGCATTGCTCAACGCAATTGCTCAGCCGGAACa       >  1:2551608/1‑69 (MQ=255)
                                                          tCGTCTGACCGAGCACTTCCACACCTGGACCAAGCACCCATTGCTCAACGCAATTGCTCAGCCGGAACa       >  1:2466432/1‑69 (MQ=255)
                                                               tGACCGAGCACTTCCACACCTGGACCAAGCACGCATTGCTCAACGCAATTGCTCAGCCGGAACAGTTtg  >  1:2852839/1‑69 (MQ=255)
                                                                   |                                                                
AACAAGACGCGCTGCGGATGGGTAAAAAAGAGCAGTTCCCGTATGTGGGTACGACCTATCGTCTGACCGAGCACTTCCACACCTGGACCAAGCACGCATTGCTCAACGCAATTGCTCAGCCGGAACAGTTTG  >  minE/935130‑935261

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