Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 937,004 A→G 37.5% V188V (GTA→GTG fdnI → formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE937,0040AG37.5% 79.5 / 92.9 141V188V (GTA→GTGfdnIformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base A (60/28);  new base G (30/23);  total (90/51)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

Only 100 of 163 total aligned reads displayed.

CATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAA  >  minE/936945‑937069
                                                           |                                                                 
cATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGc                                                          >  1:3105946/1‑69 (MQ=255)
cATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGc                                                          >  1:3170329/1‑69 (MQ=255)
cATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGc                                                          >  1:270125/1‑69 (MQ=255)
cATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGc                                                          >  1:638417/1‑69 (MQ=255)
cATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGc                                                          >  1:2357464/1‑69 (MQ=255)
cATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGc                                                          >  1:2330464/1‑69 (MQ=255)
cATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGc                                                          >  1:1712234/1‑69 (MQ=255)
cATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGc                                                          >  1:137240/1‑69 (MQ=255)
    tGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTgg                                                       >  1:2471052/1‑68 (MQ=255)
           ggCATTTTGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCgaaga                                                >  1:140437/1‑68 (MQ=255)
           ggCATTCTGGGTAAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTGAgccgtc                                                            >  1:3056682/1‑56 (MQ=255)
           ggCATTCTGGGTAAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTGAgccgtc                                                            >  1:1922723/1‑56 (MQ=255)
           ggCATTCTGGGTAAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTGAgccgtc                                                            >  1:2039018/1‑56 (MQ=255)
           ggCATTCTGGGTAAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTGAgccgtc                                                            >  1:413847/1‑56 (MQ=255)
           ggCATTCTGGGTAAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTGAgccgtc                                                            >  1:19823/1‑56 (MQ=255)
           ggCATTCTGGGTAAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTGAgccgtc                                                            >  1:3352026/1‑56 (MQ=255)
           ggCATTCTGGGTAAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTGAgccgtc                                                            >  1:1934809/1‑56 (MQ=255)
           ggCATTCTGGGTAAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTGAgccgtc                                                            >  1:1525232/1‑56 (MQ=255)
           ggCATTCTGGGTAAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTGAgccgtc                                                            >  1:1475329/1‑56 (MQ=255)
           ggCATTCTGGGTAAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTGAgccgtc                                                            >  1:478255/1‑56 (MQ=255)
           ggCATTCTGGGTAAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTGAgccgtc                                                            >  1:2180486/1‑56 (MQ=255)
           ggCATTCTGGGTAAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTGAgccgtc                                                            >  1:1324504/1‑56 (MQ=255)
           ggCATTCTGGGTAAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTGAgccgtc                                                            >  1:1154402/1‑56 (MQ=255)
           ggCATTCTGGGTAAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTGAgccgtc                                                            >  1:2651092/1‑56 (MQ=255)
           ggCATTCTGGGTAAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTGAgccgtc                                                            >  1:1007212/1‑56 (MQ=255)
           ggCATTCTGGGTAAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTGAgccgtc                                                            >  1:975815/1‑56 (MQ=255)
           ggCATTCTGGGTAAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTGAgacgtc                                                            >  1:2920305/1‑56 (MQ=255)
               ttttGGGTGAAAGGATCGATTTAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAactc                                           <  1:3233799/69‑3 (MQ=255)
               ttttGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAacac                                           <  1:935737/69‑1 (MQ=255)
               ttttGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAacac                                           <  1:803549/69‑1 (MQ=255)
               ttttGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAacac                                           <  1:27956/69‑1 (MQ=255)
               ttttGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAacac                                           <  1:868010/69‑1 (MQ=255)
               ttttGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAacac                                           <  1:2352691/69‑1 (MQ=255)
               ttttGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAacac                                           <  1:1392384/69‑1 (MQ=255)
               ttttGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAacac                                           <  1:1567370/69‑1 (MQ=255)
                tttGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAacac                                           <  1:303553/68‑1 (MQ=255)
                     gTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCcg                                     >  1:588774/1‑69 (MQ=255)
                     gTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCcg                                     >  1:2910739/1‑69 (MQ=255)
                           ggATCGATTAAAGGGATGATCGAAGGGAAGGTAAgtc                                                               >  1:378747/1‑37 (MQ=255)
                           ggATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAacac                                           >  1:3147728/1‑57 (MQ=255)
                           ggATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTgg                               <  1:3086615/69‑1 (MQ=255)
                           ggATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTgg                               >  1:1561471/1‑69 (MQ=255)
                           ggATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTgg                               >  1:2973107/1‑69 (MQ=255)
                           ggATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTgg                               >  1:227236/1‑69 (MQ=255)
                           ggATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTgg                               >  1:531296/1‑69 (MQ=255)
                           ggATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTgg                               >  1:446697/1‑69 (MQ=255)
                           ggATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTgg                               >  1:2884920/1‑69 (MQ=255)
                           ggATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTgg                               >  1:797576/1‑69 (MQ=255)
                                  ttAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACAcc                                          >  1:2650669/1‑51 (MQ=255)
                                  ttAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACAcc                                          >  1:1257100/1‑51 (MQ=255)
                                  ttAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCa                                         >  1:3412845/1‑52 (MQ=255)
                                  ttAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCa                                         >  1:2608088/1‑52 (MQ=255)
                                  ttAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCa                                         >  1:1238660/1‑52 (MQ=255)
                                  ttAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCa                                         >  1:1830769/1‑52 (MQ=255)
                                  ttAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCa                                         >  1:2569736/1‑52 (MQ=255)
                                  ttAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCa                                         >  1:1955143/1‑52 (MQ=255)
                                  ttAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCa                                         <  1:1052522/52‑1 (MQ=255)
                                  ttAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCa                                         >  1:206580/1‑52 (MQ=255)
                                  ttAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCa                                         >  1:2242485/1‑52 (MQ=255)
                                    aaaGGGATGATCGAAGGGAAGGTGAGCCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGaa                      <  1:129321/69‑1 (MQ=255)
                                    aaaGGGATGATCGAAGGGAAGGTGAGCCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGaa                      <  1:2749816/69‑1 (MQ=255)
                                    aaaGGGATGATCGAAGGGAAGGTGAGCCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGaa                      <  1:3161966/69‑1 (MQ=255)
                                    aaaGGGATGATCGAAGGGAAGGTGAGCCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGaa                      <  1:2047419/69‑1 (MQ=255)
                                    aaaGGGATGATCGAAGGGAAGGTGAGCCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGaa                      <  1:685308/69‑1 (MQ=255)
                                     aaGGGATGATCGAAGGGAAGGTGAGCCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGaa                      <  1:596273/68‑1 (MQ=255)
                                      aGGGATGGTCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGaaa                     <  1:1437216/68‑1 (MQ=255)
                                      aGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGaaa                     <  1:1161646/68‑1 (MQ=255)
                                      aGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGaaa                     <  1:1281166/68‑1 (MQ=255)
                                      aGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGaaa                     <  1:909047/68‑1 (MQ=255)
                                      aGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGaaa                     <  1:1825054/68‑1 (MQ=255)
                                       gggATGATCGAAGGGAAGGTGAGCCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATc                   <  1:1132521/69‑1 (MQ=255)
                                       gggATGATCGAAGGGAAGGTGAGCCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATc                   <  1:1894339/69‑1 (MQ=255)
                                       gggATGATCGAAGGGAAGGTGAGCCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATc                   <  1:640764/69‑1 (MQ=255)
                                        ggATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCg                  <  1:17976/69‑1 (MQ=255)
                                        ggATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCg                  <  1:2818495/69‑1 (MQ=255)
                                         gatgatCGAAGGGAAGGTAAGTCGTCGCTTGGCGAAGAAACACCATCCGCGCTGGt                              <  1:1022099/56‑1 (MQ=255)
                                          atgatCGAAGGGAAGGTGAGCCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGaa                      <  1:2306617/63‑1 (MQ=255)
                                              tCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTTGTATCTTGAAATCGAGaa              >  1:1674222/1‑67 (MQ=255)
                                              tCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGaa              >  1:1043764/1‑67 (MQ=255)
                                              tCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGaa              >  1:2963635/1‑67 (MQ=255)
                                              tCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGaa              >  1:1971857/1‑67 (MQ=255)
                                              tCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGaa              >  1:2369049/1‑67 (MQ=255)
                                              tCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGaa              >  1:2223220/1‑67 (MQ=255)
                                                   gggaaggTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGaa       >  1:94259/1‑69 (MQ=255)
                                                     gaaggTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAg      >  1:341436/1‑68 (MQ=255)
                                                          tAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGaa  >  1:750213/1‑67 (MQ=255)
                                                           aaGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTAt                            >  1:2436180/1‑40 (MQ=255)
                                                           aaGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAAt                    >  1:1622733/1‑48 (MQ=255)
                                                           aaGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGaa  >  1:3248613/1‑66 (MQ=255)
                                                           aaGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGaa  >  1:48311/1‑66 (MQ=255)
                                                           aaGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGaa  >  1:363142/1‑66 (MQ=255)
                                                           aaGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGaa  >  1:2090233/1‑66 (MQ=255)
                                                           aaGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGaa  >  1:3035103/1‑66 (MQ=255)
                                                           aaGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGaa  >  1:3007965/1‑66 (MQ=255)
                                                           aaGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGaa  >  1:2922942/1‑66 (MQ=255)
                                                           aaGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGaa  >  1:2876266/1‑66 (MQ=255)
                                                           aaGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGaa  >  1:2841681/1‑66 (MQ=255)
                                                           aaGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGaa  >  1:245188/1‑66 (MQ=255)
                                                           aaGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGaa  >  1:2726360/1‑66 (MQ=255)
                                                           aaGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGaa  >  1:2493233/1‑66 (MQ=255)
                                                           |                                                                 
CATATGTATATGGCATTTTGGGTGAAAGGATCGATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAA  >  minE/936945‑937069

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: