Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,443 A→G 23.9% 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,4430AG23.9% 80.3 / 23.9 46V387V (GTA→GTGnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base A (19/16);  new base G (6/5);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCAGACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGG  >  minE/836376‑836513
                                                                   |                                                                      
cagcagACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTa                                                                        <  1:3711561/68‑1 (MQ=255)
cagcagACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTa                                                                        <  1:1729593/68‑1 (MQ=255)
    agACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATc                                                                 <  1:3849189/71‑1 (MQ=255)
       cgcgGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTgg                                                              >  1:644461/1‑71 (MQ=255)
            gTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAAt                                                         <  1:2810159/71‑1 (MQ=255)
            gTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAAt                                                         >  1:1101669/1‑71 (MQ=255)
                 cTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATctgct                                                    <  1:976080/71‑1 (MQ=255)
                     gTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAc                                                 >  1:1452578/1‑70 (MQ=255)
                      tAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAccg                                               >  1:1734793/1‑71 (MQ=255)
                          aaCGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTa                                                                   <  1:608002/47‑1 (MQ=255)
                           aCGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGt                                          <  1:1073971/71‑1 (MQ=255)
                            cGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGCCCGCCgg                                           >  1:2609738/1‑69 (MQ=255)
                            cGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGt                                          >  1:641943/1‑70 (MQ=255)
                            cGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGt                                          >  1:3186378/1‑70 (MQ=255)
                            cGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGt                                          >  1:2610465/1‑70 (MQ=255)
                             ggCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACa                                                                      >  1:2822896/1‑41 (MQ=255)
                             ggCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAAt                                                         >  1:1453048/1‑54 (MQ=255)
                              gCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGAt                                       <  1:2401520/71‑1 (MQ=255)
                                 ttCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTAt                                                                  <  1:2401522/41‑1 (MQ=255)
                                 tgcTGCCCGCCATCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATAcc                                    >  1:3664896/3‑71 (MQ=255)
                                   cTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACa                                                                      >  1:2082983/1‑35 (MQ=255)
                                   cTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCa                                   <  1:1266547/70‑1 (MQ=255)
                                      ccAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCCGTGATACCAAAc                                <  1:2405206/70‑1 (MQ=255)
                                       cAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCg                              <  1:426982/71‑1 (MQ=255)
                                        aGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATct                                                       <  1:1088336/45‑1 (MQ=255)
                                        aGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCgg                             <  1:2976921/71‑1 (MQ=255)
                                         gACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAccgc                                              <  1:2530743/53‑1 (MQ=255)
                                         gACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAccg                                               >  1:544649/1‑52 (MQ=255)
                                            gTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACt                         <  1:3147642/71‑1 (MQ=255)
                                            atCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTa                           >  1:2848576/2‑69 (MQ=255)
                                             tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGa                                                  >  1:896972/1‑45 (MQ=16)
                                             tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACt                         <  1:1755072/70‑1 (MQ=255)
                                              cGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgc                       >  1:842499/1‑71 (MQ=255)
                                               gAAAGCCTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTgc                       <  1:2414884/70‑1 (MQ=255)
                                                  aGCCTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGa                                        >  1:236314/1‑50 (MQ=255)
                                                  aGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgcgcg                    >  1:2001257/1‑70 (MQ=255)
                                                    ccTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTgcgggcgc                 <  1:711000/71‑3 (MQ=255)
                                                       gCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTg               >  1:68702/1‑70 (MQ=255)
                                                       gCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCgg                             >  1:2983601/1‑56 (MQ=255)
                                                         gTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTgcggg                    >  1:23711/1‑61 (MQ=255)
                                                           aTTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACg          >  1:721657/1‑71 (MQ=255)
                                                            ttCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACg          <  1:3651336/70‑1 (MQ=255)
                                                            ttCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACg          <  1:2141125/70‑1 (MQ=255)
                                                            ttCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACg          <  1:2835664/70‑1 (MQ=255)
                                                             tCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTa        >  1:518274/1‑71 (MQ=255)
                                                              ccGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTAt       >  1:3076874/1‑71 (MQ=255)
                                                                ggTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGc     >  1:123649/1‑71 (MQ=255)
                                                                 gTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCt    <  1:2599427/71‑1 (MQ=255)
                                                                   aCAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTgg  >  1:3529363/1‑71 (MQ=255)
                                                                   aCAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTgg  >  1:3192528/1‑71 (MQ=255)
                                                                   aCAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTgg  >  1:721740/1‑71 (MQ=255)
                                                                   |                                                                      
CAGCAGACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGG  >  minE/836376‑836513

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