Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 932,920 T→C 63.7% A139A (GCT→GCC fdnG → formate dehydrogenase‑N, alpha subunit, nitrate‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE932,9200TC63.7% 33.9 / 30.6 55A139A (GCT→GCCfdnGformate dehydrogenase‑N, alpha subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base T (8/12);  new base C (23/12);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 9.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGG  >  minE/932885‑932988
                                   |                                                                    
gaagCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAg                                     <  1:1131121/69‑1 (MQ=255)
gaagCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAg                                     <  1:1389734/69‑1 (MQ=255)
gaagCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAg                                     <  1:910470/69‑1 (MQ=255)
gaagCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAg                                     <  1:1930266/69‑1 (MQ=255)
           ccGTATTGCGAAGCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGt                             <  1:1040921/66‑1 (MQ=255)
           ccGTATTGCGAAGCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGt                             <  1:3369292/66‑1 (MQ=255)
           ccGTATTGCGAAGCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGt                             <  1:3214593/66‑1 (MQ=255)
           ccGTATTGCGAAGCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGt                             <  1:2842310/66‑1 (MQ=255)
           ccGTATTGCGAAGCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGt                             <  1:2499608/66‑1 (MQ=255)
           ccGTATTGCGAAGCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGt                             <  1:2358301/66‑1 (MQ=255)
           ccGTATTGCGAAGCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGt                             <  1:2074575/66‑1 (MQ=255)
           ccGTATTGCGAAGCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGt                             <  1:2048451/66‑1 (MQ=255)
           ccGTATTGCGAAGCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGt                             <  1:790887/66‑1 (MQ=255)
           ccGTATTGCGAAGCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGt                             <  1:1693344/66‑1 (MQ=255)
           ccGTATTGCGAAGCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGt                             <  1:879677/66‑1 (MQ=255)
           ccGTATTGCGAAGCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGt                             <  1:863015/66‑1 (MQ=255)
               aTTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTa                      <  1:1196379/69‑1 (MQ=255)
               aTTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTa                      <  1:3416980/69‑1 (MQ=255)
               aTTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTa                      <  1:2634076/69‑1 (MQ=255)
               aTTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTa                      >  1:2469241/1‑69 (MQ=255)
               aTTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTa                      <  1:2409907/69‑1 (MQ=255)
               aTTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTa                      <  1:491382/69‑1 (MQ=255)
                ttGCGAAGCTGATGAAAGCTGCCCGTGACGCTAACTTTATTGAAAAg                                           >  1:382574/1‑47 (MQ=255)
                ttGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTaa                     <  1:2410621/69‑1 (MQ=255)
                ttGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTaa                     >  1:3383686/1‑69 (MQ=255)
                ttGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTaa                     >  1:2917339/1‑69 (MQ=255)
                ttGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTaa                     <  1:1349736/69‑1 (MQ=255)
                ttGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTaa                     >  1:1350451/1‑69 (MQ=255)
                ttGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTaa                     >  1:1479330/1‑69 (MQ=255)
                ttGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTaa                     >  1:990405/1‑69 (MQ=255)
                ttGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTaa                     <  1:1969604/69‑1 (MQ=255)
                ttGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTaa                     >  1:856104/1‑69 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGa                                          >  1:664485/1‑41 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:434080/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:526166/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:687689/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:904682/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:3431626/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:809056/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:1128960/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:316251/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:3161279/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:2526276/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:2413681/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:2312340/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:2234901/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:2046107/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:2042261/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:1892956/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:1761751/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:1691285/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:1680135/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:1614503/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGtt               >  1:158754/1‑68 (MQ=255)
                       gCTGATGAAAGCCGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGt                >  1:2460530/1‑67 (MQ=255)
                                   tGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCgg  >  1:219546/1‑69 (MQ=255)
                                   |                                                                    
GAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGG  >  minE/932885‑932988

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: