Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,446 G→A 22.0% Q388Q (CAG→CAA 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,4460GA22.0% 98.3 / 23.1 50Q388Q (CAG→CAAnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (21/18);  new base A (6/5);  total (27/23)
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

AGACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCG  >  minE/836380‑836515
                                                                  |                                                                     
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)
                         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)
                               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: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:3529363/1‑71 (MQ=255)
                                                               aCAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTgg    >  1:3192528/1‑71 (MQ=255)
                                                               aCAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTgg    >  1:721740/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)
                                                                  |                                                                     
AGACGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCG  >  minE/836380‑836515

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: