Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,284 A→G 19.3% A334A (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,2840AG19.3% 116.5 / 21.0 57A334A (GCA→GCGnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base A (22/24);  new base G (4/7);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 7.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

GGCGATTAAAGACGGTATTCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGT  >  minE/836215‑836349
                                                                     |                                                                 
ggCGATTAAAGACGGTATTCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCaa                                                                  <  1:1635607/71‑1 (MQ=255)
    aTTAAAGACGGTATTCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCaatgg                                                               <  1:5605430/70‑1 (MQ=255)
      tAAAGACGGTATTCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCaatggaat                                                            >  1:4002175/1‑71 (MQ=255)
       aaaGACGGTATTCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCaatggaatg                                                           >  1:1518194/1‑71 (MQ=255)
       aaaGACGGTATTCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCaatggaatg                                                           >  1:1494951/1‑71 (MQ=255)
         aGACGGTATTCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCaat                                                                 >  1:1561358/1‑63 (MQ=255)
         aGACGGTATTCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCaat                                                                 >  1:2592401/1‑63 (MQ=255)
         aGACGGTATTCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCaat                                                                 >  1:4236583/1‑63 (MQ=255)
         aGACGGTATTCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAatggaatgga                                                         <  1:5009456/71‑1 (MQ=255)
             ggTATTCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACt                                                     <  1:1338139/71‑1 (MQ=255)
             ggTATTCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACt                                                     <  1:1109023/71‑1 (MQ=255)
                aTTCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGaaa                                                       >  1:2772955/1‑66 (MQ=255)
                aTTCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGaaa                                                       >  1:4799966/1‑66 (MQ=255)
                aTTCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGaaa                                                       >  1:1062226/1‑66 (MQ=255)
                aTTCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTgg                                                   >  1:4791987/1‑70 (MQ=255)
                aTTCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTg                                                    <  1:62726/69‑1 (MQ=255)
                aTTCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTg                                                    <  1:233653/69‑1 (MQ=255)
                 ttCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGTaaa                                                       <  1:2730733/65‑1 (MQ=255)
                 ttCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGgcg                                                 <  1:3288644/71‑1 (MQ=255)
                      cTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGcgctgcc                                            >  1:526092/1‑71 (MQ=255)
                         aGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCaa                                                                  >  1:5470080/1‑46 (MQ=255)
                         aGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGcgctgccgct                                         >  1:3569823/1‑71 (MQ=255)
                         aGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGcgctgccgc                                          >  1:5205604/1‑70 (MQ=255)
                             tgaTTGACGCTGCCCAGCGTTCGCCAGTCTACAAAATGGCGATGGACTGGAAACTGGCACTACCGTTGCa                                      <  1:882684/68‑1 (MQ=255)
                               aTTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCaa                                                                  >  1:4375477/1‑40 (MQ=255)
                               aTTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTgg                                                   >  1:5401080/1‑55 (MQ=255)
                               aTTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATcc                                   >  1:4703232/1‑71 (MQ=255)
                                ttGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCaatggaatgg                                                          >  1:5291052/1‑47 (MQ=255)
                                 tGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCaatgg                                                               <  1:5098649/41‑1 (MQ=255)
                                 tGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGgcgc                                                <  1:2782649/56‑1 (MQ=255)
                                    aGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCa                                      <  1:4263137/63‑1 (MQ=255)
                                    aGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCa                                      <  1:5246044/63‑1 (MQ=255)
                                       cGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATc                           <  1:5504348/71‑1 (MQ=255)
                                         cacaGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGaa                               >  1:1356357/1‑65 (MQ=255)
                                          acaGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGca                        <  1:1276949/71‑1 (MQ=255)
                                           cagcagTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGcac                       <  1:4454388/71‑1 (MQ=255)
                                           cagcagTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGcac                       <  1:1163351/71‑1 (MQ=255)
                                              cagTCGACGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACt                    >  1:1075077/1‑71 (MQ=255)
                                               gttCGCCAGTCTACAAAATGGCGATGGACTGGAAACTGGCACTACCGTTGCACCCTGAATACCGcacc                      >  1:3885699/3‑67 (MQ=255)
                                                   gCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGcc                 >  1:62041/1‑69 (MQ=255)
                                                   gCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGa               <  1:1814838/71‑1 (MQ=255)
                                                   gCCAGTCTACAATATGGCGATGGACTGGAAACTGGCACTACCGTTGCACCCTGAATACCGCACCCTGCCGa               <  1:4152390/71‑1 (MQ=255)
                                                      ggTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATgg            >  1:1938163/1‑71 (MQ=255)
                                                        tttATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCg                          <  1:1365978/55‑1 (MQ=255)
                                                        tttATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCg                          <  1:462214/55‑1 (MQ=255)
                                                        tttATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCg                          <  1:1955855/55‑1 (MQ=255)
                                                        tttATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCacac                     >  1:5177414/1‑60 (MQ=255)
                                                         ttATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCt         >  1:3313728/1‑71 (MQ=255)
                                                          tataAAATGGCAATGGAATGGAAACTGGcgctgccgctgc                                       >  1:5595983/1‑40 (MQ=255)
                                                          tataAAATGGCAATGGAATGGAAACTGGcgctgccgctgc                                       >  1:2803711/1‑40 (MQ=255)
                                                          tataAAATGGCAATGGAATGGAAACTGGcgctgccgctgc                                       >  1:4928706/1‑40 (MQ=255)
                                                          tataAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGaa                               >  1:5086200/1‑48 (MQ=255)
                                                          tacaaaaTGGCGATGGACTGGAAACTGGCACTACCGTTGCACCCTGAATACCGCACCCTGCCGATGGTCTg        >  1:5603427/4‑71 (MQ=255)
                                                          tacaaaaTGGCGATGGACTGGAAACTGGCACTACCGTTGCACCCTGAATACCGCACCCTGCCGATGGTCTg        >  1:2633138/4‑71 (MQ=255)
                                                          tacaaaaTGGCGATGGACTGGAAACTGGCACTACCGTTGCACCCTGAATACCGCACCCTGCCGATGGTCTg        >  1:3339406/4‑71 (MQ=255)
                                                            taAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTgg       >  1:2601604/1‑70 (MQ=255)
                                                            taAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGt      <  1:2329371/71‑1 (MQ=255)
                                                            taAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGt      <  1:4416471/71‑1 (MQ=255)
                                                             aaaaTGGCGATGGACTGGAAACTGGCACTACCGTTGCAccctg                                 <  1:991576/43‑3 (MQ=12)
                                                             aaaaTGGCGATGGACTGGAAACTGGCACTACCGTTGCAccctg                                 <  1:4885262/43‑3 (MQ=12)
                                                             aaaaTGGCGATGGACTGGAAACTGGCACTACCGTTGCAccctg                                 <  1:3595060/43‑3 (MQ=12)
                                                             aaaaTGGCGATGGACTGGAAACTGGCACTACCGTTGCAccctg                                 <  1:3262412/43‑3 (MQ=12)
                                                              aaaTGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCg                >  1:3084986/1‑59 (MQ=255)
                                                              aaaTGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTAc    >  1:5058895/1‑71 (MQ=255)
                                                                aTGGCGATGGACTGGAAACTGGCACTACCGTTGCACCCTGAATACCGCACCCTGCCGATGGTCTGGTACGt  <  1:1080605/71‑1 (MQ=255)
                                                                    cAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAat                              <  1:4996949/39‑1 (MQ=255)
                                                                    cAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAat                              <  1:2842852/39‑1 (MQ=255)
                                                                    cAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAat                              <  1:1162407/39‑1 (MQ=255)
                                                                    cAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAat                              <  1:709536/39‑1 (MQ=255)
                                                                     aatggaatggAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGa               <  1:5294670/53‑1 (MQ=255)
                                                                     aatggaatggAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGa               <  1:3671006/53‑1 (MQ=255)
                                                                     aatggaatggAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGa               <  1:2309694/53‑1 (MQ=255)
                                                                     aatggaatggAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGa               <  1:1083093/53‑1 (MQ=255)
                                                                     |                                                                 
GGCGATTAAAGACGGTATTCCGCTGAGCGTTATTGAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGT  >  minE/836215‑836349

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