Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,497 A→G 7.4% A405A (GCA→GCG 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,4970AG 7.4% 184.7 / 4.7 67A405A (GCA→GCGnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base A (32/30);  new base G (3/2);  total (35/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 = 1.00e+00

GCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATACCCGT  >  minE/836431‑836566
                                                                  |                                                                     
gCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTg                                                                    >  1:68702/1‑70 (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: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)
                  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)
                                    cgccgGTGATACCAAACCGGTACTGCGCGCACTGAAAc                                                                >  1:1782418/1‑38 (MQ=255)
                                     gccgGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGa                              <  1:3328494/71‑1 (MQ=255)
                                      ccgGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGa                              >  1:1654898/1‑70 (MQ=255)
                                        ggTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAAc                           >  1:2637478/1‑71 (MQ=255)
                                        ggTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAAc                           >  1:3136597/1‑71 (MQ=255)
                                            aTACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGtt                       <  1:2773572/71‑1 (MQ=255)
                                            aTACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGtt                       <  1:3892482/71‑1 (MQ=255)
                                            aTACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGtt                       <  1:2858762/71‑1 (MQ=255)
                                            aTACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGtt                       <  1:2069211/71‑1 (MQ=255)
                                            aTACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGtt                       <  1:663707/71‑1 (MQ=255)
                                            aTACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGtt                       <  1:3222740/71‑1 (MQ=255)
                                             tACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGt                                   <  1:2990046/58‑1 (MQ=255)
                                              aCCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGAt                                                   >  1:1998264/1‑41 (MQ=255)
                                              aCCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGAt                                                   >  1:675893/1‑41 (MQ=255)
                                              aCCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCg                         <  1:2118596/67‑1 (MQ=255)
                                               ccAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACg                                    >  1:2734011/1‑55 (MQ=255)
                                               ccAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACg                                    >  1:2390871/1‑55 (MQ=255)
                                               ccAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGa                     <  1:3570913/70‑1 (MQ=255)
                                               ccAAACCGGTACTGCGCGCACTGAAACCTATGCTGGc                                                      <  1:1992818/37‑1 (MQ=255)
                                                cAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACg                   >  1:1338379/1‑71 (MQ=255)
                                                cAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACg                   >  1:1592283/1‑71 (MQ=255)
                                                cAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACg                   >  1:2032387/1‑71 (MQ=255)
                                                 aaaCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACa                                        <  1:3699867/49‑1 (MQ=255)
                                                 aaaCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACgg                  >  1:1705896/1‑71 (MQ=255)
                                                  aaCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGt                                               >  1:617399/1‑41 (MQ=255)
                                                  aaCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACgg                  >  1:3178845/1‑70 (MQ=255)
                                                  aaCCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACgg                  >  1:837261/1‑70 (MQ=255)
                                                    ccGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTaa               <  1:185515/71‑1 (MQ=255)
                                                    ccGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTa                <  1:2467750/70‑1 (MQ=255)
                                                      ggTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTa                >  1:3566989/1‑68 (MQ=255)
                                                      ggTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAg             <  1:3332000/71‑1 (MQ=255)
                                                         aCTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACgg                  >  1:74235/1‑63 (MQ=255)
                                                          cTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTg          <  1:2489591/70‑1 (MQ=255)
                                                           tGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGa         <  1:937292/70‑1 (MQ=255)
                                                           tGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGAAAAGt            >  1:2052142/1‑67 (MQ=255)
                                                            gcgcgcACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATa       >  1:1856843/1‑71 (MQ=255)
                                                            gcgcgcACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATa       >  1:1713103/1‑71 (MQ=255)
                                                               cgcACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATAccc    <  1:844213/71‑1 (MQ=255)
                                                                 cACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAAcc                          <  1:2583551/47‑1 (MQ=255)
                                                                 cACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGtt           <  1:2926113/62‑1 (MQ=255)
                                                                 cACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATACCCGt  <  1:654493/71‑1 (MQ=255)
                                                                 cACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATACCCGt  <  1:3522257/71‑1 (MQ=255)
                                                                  |                                                                     
GCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGCGTCATTACAAACGTGCTGAAACCGTTGACGGTAAAGTTGATACCCGT  >  minE/836431‑836566

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: