Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,728 C→G 12.7% F67L (TTC→TTG narI → nitrate reductase 1, gamma (cytochrome b(NR)) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE837,7280CG12.7% 96.3 / 4.9 39F67L (TTC→TTGnarInitrate reductase 1, gamma (cytochrome b(NR)) subunit
Reads supporting (aligned to +/- strand):  ref base C (21/13);  new base G (2/3);  total (23/16)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAA  >  minE/837659‑837798
                                                                     |                                                                      
tCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCActtct                                                                       >  1:2943459/1‑71 (MQ=255)
tCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCActtct                                                                       <  1:3077969/71‑1 (MQ=255)
    aaaGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCgg                                                                   >  1:2214863/1‑71 (MQ=255)
     aaGGGATGAACCTGGCGTCGAACCTGTTCCATATAGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCgg                                                                   >  1:3945069/1‑70 (MQ=255)
       gggATGAACCTGGCGTCGCACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCgg                                                                   >  1:338379/1‑68 (MQ=255)
         gATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGc                                                              <  1:2548172/71‑1 (MQ=255)
          aTGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt                                                             >  1:498578/1‑71 (MQ=255)
          aTGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt                                                             >  1:174865/1‑71 (MQ=255)
               ccTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgc                                                        >  1:2300660/1‑71 (MQ=255)
                cTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgc                                                        <  1:3764346/70‑1 (MQ=255)
                  ggCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgccg                                                      >  1:1505350/1‑70 (MQ=255)
                   gCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgccg                                                      <  1:2446114/69‑1 (MQ=255)
                      tCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTg                                                 <  1:2020350/71‑1 (MQ=255)
                        gAATCTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTgg                                                <  1:999754/70‑1 (MQ=255)
                        gAATCTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTGCGGCATGTTAACGCCGCACTg                                                 <  1:1437546/69‑1 (MQ=255)
                        gAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGa                                               >  1:1273085/1‑71 (MQ=255)
                          aCCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATg                                             >  1:1463209/1‑71 (MQ=255)
                           tcTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGATg                                             >  1:2560582/2‑70 (MQ=255)
                             tGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTgg                                                <  1:1092000/65‑1 (MQ=255)
                               ttCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGa                                        <  1:1064298/71‑1 (MQ=255)
                                tCCATTTCGGGATTCTGGGGATTTTTGTCGGTCAGTTCTTCGGTATGCTGACGCCGCACTGGATGTATGa                                        <  1:2305298/70‑1 (MQ=255)
                                tCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCg                                                                    >  1:3939225/1‑42 (MQ=255)
                                tCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATg                                         >  1:3013147/1‑69 (MQ=255)
                                     atCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGa                                                           <  1:2855573/46‑1 (MQ=255)
                                               tGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGATGTa                                           >  1:723555/1‑52 (MQ=255)
                                                 gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCt                               >  1:1730880/1‑62 (MQ=255)
                                                 gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCt                               >  1:210283/1‑62 (MQ=255)
                                                 gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGa                      <  1:766058/71‑1 (MQ=255)
                                                  ggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGc                                     >  1:941868/1‑55 (MQ=255)
                                                    aTTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGt                   <  1:1054019/71‑1 (MQ=255)
                                                    aTTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGt                   <  1:994830/71‑1 (MQ=255)
                                                     tttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGt                   >  1:3052113/1‑70 (MQ=255)
                                                     tttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGt                   >  1:3873386/1‑70 (MQ=255)
                                                        ttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGaaa               >  1:2970484/1‑71 (MQ=255)
                                                        ttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGaa                <  1:2830759/70‑1 (MQ=255)
                                                         tGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCt                                   >  1:3248328/1‑50 (MQ=255)
                                                         tGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCt                                   >  1:3079998/1‑50 (MQ=255)
                                                           tCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGaaa               >  1:2326049/1‑68 (MQ=255)
                                                               tCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGaa                                       <  1:403390/40‑1 (MQ=255)
                                                               tCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGaa                                       <  1:304793/40‑1 (MQ=255)
                                                                 aCTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATg      >  1:608373/1‑71 (MQ=255)
                                                                   ttcttcGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATgg     >  1:3934213/1‑70 (MQ=255)
                                                                     cttcGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCaa  >  1:3931720/1‑71 (MQ=255)
                                                                     |                                                                      
TCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAA  >  minE/837659‑837798

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: