Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,466 C→G 24.9% T395S (ACC→AGC) ‡ 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,4660CG24.9% 82.1 / 23.9 48T395S (ACC→AGC) ‡narHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base C (19/17);  new base G (7/5);  total (26/22)
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

GTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGC  >  minE/836397‑836535
                                                                     |                                                                     
gTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAc                                                                       >  1:1452578/1‑70 (MQ=255)
 tAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAccg                                                                     >  1:1734793/1‑71 (MQ=255)
      aCGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGt                                                                <  1:1073971/71‑1 (MQ=255)
       cGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGCCCGCCgg                                                                 >  1:2609738/1‑69 (MQ=255)
       cGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGt                                                                >  1:3186378/1‑70 (MQ=255)
       cGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGt                                                                >  1:2610465/1‑70 (MQ=255)
       cGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGt                                                                >  1:641943/1‑70 (MQ=255)
         gCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGAt                                                             <  1:2401520/71‑1 (MQ=255)
            tgcTGCCCGCCATCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATAcc                                                          >  1:3664896/3‑71 (MQ=255)
              cTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCa                                                         <  1:1266547/70‑1 (MQ=255)
                 ccAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCCGTGATACCAAAc                                                      <  1:2405206/70‑1 (MQ=255)
                  cAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCg                                                    <  1:426982/71‑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)
                        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:2141125/70‑1 (MQ=255)
                                       ttCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACg                                <  1:3651336/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:721740/1‑71 (MQ=255)
                                              aCAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTgg                        >  1:3192528/1‑71 (MQ=255)
                                              aCAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTgg                        >  1:3529363/1‑71 (MQ=255)
                                               cAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGc                       >  1:1492931/1‑71 (MQ=255)
                                               cAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGc                       <  1:1249856/71‑1 (MQ=255)
                                                aGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCg                      <  1:1119820/71‑1 (MQ=255)
                                                  tATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCa                                        <  1:1286605/51‑1 (MQ=255)
                                                  tATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACGGATGATGGCGAt                    >  1:263630/1‑71 (MQ=255)
                                                  tATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTGAAACGGATGATGGCGAt                    >  1:140326/1‑71 (MQ=255)
                                                   aTCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGt                                                  <  1:1456546/40‑1 (MQ=255)
                                                    tCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGc                           <  1:190613/62‑1 (MQ=255)
                                                      tGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCa                                        >  1:535422/1‑47 (MQ=255)
                                                               tgctgACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACaaa       >  1:3777084/1‑71 (MQ=255)
                                                                 ctgACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACg     >  1:1953207/1‑71 (MQ=255)
                                                                  tgACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAAc      >  1:1255969/1‑69 (MQ=255)
                                                                   gACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACg                                <  1:3094015/42‑1 (MQ=255)
                                                                    aCCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTg   <  1:862594/70‑1 (MQ=255)
                                                                    aCCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGc  <  1:351318/71‑1 (MQ=255)
                                                                     ccgccgGTGATACCAAACCGGTACTGCGCGCACTGAAACGTa                              >  1:2808009/1‑42 (MQ=255)
                                                                     |                                                                     
GTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGC  >  minE/836397‑836535

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: