Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,443 A→G 6.4% V387V (GTA→GTG narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE836,4430AG 6.4% 177.8 / 3.7 62V387V (GTA→GTGnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base A (28/30);  new base G (2/2);  total (30/32)
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.46e-01

TGCAGCAGACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTG  >  minE/836374‑836512
                                                                     |                                                                     
tGCAGCAGACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTa                                                                       >  1:534009/1‑70 (MQ=255)
tGCAGCAGACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTa                                                                       >  1:2159994/1‑70 (MQ=255)
  cagcagACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTa                                                                       <  1:993127/68‑1 (MQ=255)
  cagcagACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACa                                                                     <  1:4062474/70‑1 (MQ=255)
  cagcagACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACa                                                                     <  1:2608575/70‑1 (MQ=255)
  cagcagACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAg                                                                    <  1:460520/71‑1 (MQ=255)
  cagcagACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAg                                                                    >  1:3973951/1‑71 (MQ=255)
   agcagACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGt                                                                   <  1:3938332/71‑1 (MQ=255)
   agcagACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGt                                                                   <  1:3523021/71‑1 (MQ=255)
     cagACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTa                                                                  <  1:1149689/70‑1 (MQ=255)
         cgcgGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTgg                                                             >  1:2122209/1‑71 (MQ=255)
          gcgGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGc                                                            >  1:3698168/1‑71 (MQ=255)
            gggTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCg                                                           >  1:3883790/1‑70 (MQ=255)
            gggTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGa                                                          <  1:2927223/71‑1 (MQ=255)
                      ggTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATctgctg                                                  <  1:1047054/69‑1 (MQ=255)
                      ggTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATctgctg                                                  <  1:1042115/69‑1 (MQ=255)
                        tAGCAACGGCATTCTGCCAGACGTAGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATc                                                       >  1:2414971/1‑62 (MQ=255)
                          gCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATctg                                                     <  1:3048575/62‑1 (MQ=255)
                            aaCGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAAt                                                        >  1:3933056/1‑57 (MQ=255)
                            aaCGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAAt                                                        >  1:3515285/1‑57 (MQ=255)
                              cgtcATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTg                                        <  1:3695159/68‑1 (MQ=255)
                               ggCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGa                                       <  1:2825019/71‑1 (MQ=255)
                               ggCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGa                                       <  1:3049400/71‑1 (MQ=255)
                                 cATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATa                                     >  1:324108/1‑71 (MQ=255)
                                   ttCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGt                                                                   <  1:27138/39‑1 (MQ=255)
                                     cTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTg                                        >  1:2542316/1‑64 (MQ=255)
                                     cTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCa                                  >  1:1357481/1‑70 (MQ=255)
                                      tGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATAc                                    >  1:1358649/1‑67 (MQ=255)
                                      tGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCaaa                                >  1:1501839/1‑71 (MQ=255)
                                        ccAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAAc                               <  1:763241/70‑1 (MQ=255)
                                          aGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCt                                                               >  1:1184506/1‑36 (MQ=255)
                                          aGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGt                                         >  1:2182173/1‑58 (MQ=255)
                                          aGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCgg                            <  1:3690589/71‑1 (MQ=255)
                                            aCGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAAtcttct                                                   <  1:1505454/46‑1 (MQ=255)
                                               tCGAAAGCCTGTGTATTCCGGTGCAATATCTCGCCaa                                                         <  1:759919/37‑1 (MQ=16)
                                               tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCaa                                                         <  1:468029/37‑1 (MQ=21)
                                               tCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGt                                         <  1:900946/53‑1 (MQ=255)
                                                 gAAAGCCTGCGTATTCCGGTGCAATATCTCGCAAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTgcg                     >  1:305060/1‑71 (MQ=255)
                                                 gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGaa                                                         >  1:1727636/1‑35 (MQ=255)
                                                 gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATctgctg                                                  >  1:2621622/1‑42 (MQ=255)
                                                 gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATctg                                                     <  1:573089/39‑1 (MQ=255)
                                                 gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAccgcc                                            <  1:599095/48‑1 (MQ=255)
                                                 gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAccgcc                                            <  1:594163/48‑1 (MQ=255)
                                                 gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATAc                                    >  1:2178356/1‑56 (MQ=255)
                                                 gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgcg                     >  1:1760449/1‑71 (MQ=255)
                                                   aaGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAccgcc                                            >  1:3962629/1‑46 (MQ=255)
                                                   aaGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgcgcg                   >  1:70665/1‑71 (MQ=255)
                                                    aGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgcgcgc                  <  1:2966477/71‑1 (MQ=255)
                                                    aGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgcgcg                   >  1:3622679/1‑70 (MQ=255)
                                                    aGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGATGATACCAAACCGGTACTgcgcg                   >  1:3739535/1‑70 (MQ=255)
                                                      ccTGCGTATTCCGGTACAGTATCTGGCGAATCTGGTGACCGCCGGTGATACCAAACCGGTACTGCGCGCAc                <  1:60223/71‑1 (MQ=255)
                                                       cTGCGTATTACGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACt               <  1:2159068/71‑1 (MQ=255)
                                                         gCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTg              >  1:1951945/1‑70 (MQ=255)
                                                          cGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCaaa                                >  1:3933348/1‑51 (MQ=255)
                                                          cGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCgg                            <  1:1218914/55‑1 (MQ=255)
                                                          cGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCgg                            <  1:3096074/55‑1 (MQ=255)
                                                              ttCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgcgc                    >  1:840668/1‑59 (MQ=255)
                                                              ttCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACg         >  1:578007/1‑70 (MQ=255)
                                                              ttCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACg         >  1:1156982/1‑70 (MQ=255)
                                                                ccGGTACAGTATCTGGCGCATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTa       >  1:1809069/1‑70 (MQ=255)
                                                                ccGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATAc                                    >  1:175459/1‑41 (MQ=255)
                                                                  ggTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGc    <  1:412767/71‑1 (MQ=255)
                                                                  ggTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGc    <  1:36380/71‑1 (MQ=255)
                                                                  ggTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGc    <  1:172290/71‑1 (MQ=255)
                                                                    tACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTg  <  1:632821/71‑1 (MQ=255)
                                                                     |                                                                     
TGCAGCAGACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTG  >  minE/836374‑836512

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: