Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,115 A→G 18.7% P191P (CCA→CCG narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,1150AG18.7% 170.3 / 33.2 80P191P (CCA→CCGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (37/28);  new base G (4/11);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 4.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

TGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCT  >  minE/832048‑832185
                                                                   |                                                                      
tGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAAtt                                                                     >  1:1041434/1‑71 (MQ=255)
    tgccgccGCTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCgg                                                                 >  1:935703/2‑71 (MQ=255)
    cgccgcATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCg                                                                  >  1:1298201/1‑70 (MQ=255)
    cgccgcATCTAACGTTTACACCATCAAAAACTACGGCCCCGACCGTGTTGCTGGTTTCTCGCCAATTCCgg                                                                 >  1:3104682/1‑71 (MQ=255)
     gccgccGCTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCgg                                                                 <  1:3031776/70‑1 (MQ=255)
        gcATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAAt                                                             <  1:2572548/71‑1 (MQ=255)
        gcATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAAt                                                             <  1:2560064/71‑1 (MQ=255)
         cATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATg                                                            >  1:3678561/1‑71 (MQ=255)
          cgcTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGt                                                           <  1:3671646/69‑1 (MQ=255)
          cgcTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGt                                                           <  1:636593/69‑1 (MQ=255)
            cTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTc                                                          >  1:3169106/1‑70 (MQ=255)
               aCGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCg                                                         <  1:3300731/68‑1 (MQ=255)
                   ttACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGtt                                                   >  1:1988771/1‑70 (MQ=255)
                    tACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGttt                                                  <  1:866329/70‑1 (MQ=255)
                    tACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTc                                                 >  1:3519133/1‑71 (MQ=255)
                      gacCATCAAAACTTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGTCGATGGttt                                                  <  1:2066270/67‑1 (MQ=255)
                      gacCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCtt                                               <  1:1990199/70‑1 (MQ=255)
                      cacCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCtt                                               <  1:1292348/71‑1 (MQ=255)
                       acCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCt                                                <  1:3504019/69‑1 (MQ=255)
                       acCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTa                                              >  1:3564092/1‑71 (MQ=255)
                        ccATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTAc                                             <  1:3482110/71‑1 (MQ=255)
                         cATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCTTACg                                            <  1:2068268/71‑1 (MQ=255)
                             aaaacCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCTTACGccg                                         <  1:757389/70‑3 (MQ=255)
                               aaaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCaa                                                              >  1:2730605/1‑47 (MQ=255)
                               aaaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCaa                                                              >  1:2802201/1‑47 (MQ=255)
                               aaaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCgg                                      <  1:3312136/71‑1 (MQ=255)
                               aaaCTACGGCCCGGACCGTGTTGCTGGTTTATCGCCAATTCCGGCaa                                                              >  1:942336/1‑47 (MQ=255)
                                aaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCggg                                     <  1:3059410/71‑1 (MQ=255)
                                aaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCgg                                      >  1:371405/1‑70 (MQ=255)
                                  cTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCtt                                               <  1:2516922/59‑1 (MQ=255)
                                    aCGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTa                                              >  1:204929/1‑58 (MQ=255)
                                     cGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGTCGATGGtt                                                   >  1:1641533/1‑52 (MQ=255)
                                     cGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCa                                 <  1:1010444/70‑1 (MQ=255)
                                      ggCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGa                                                        >  1:3263008/1‑46 (MQ=255)
                                      ggCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGa                                                        >  1:2791558/1‑46 (MQ=255)
                                       gCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGtt                                                   <  1:2548387/50‑1 (MQ=255)
                                         ccGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTa                            >  1:436276/1‑71 (MQ=255)
                                           ggACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtc                          >  1:3637623/1‑71 (MQ=255)
                                           ggACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtc                          >  1:616363/1‑71 (MQ=255)
                                            gACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTc                                                          <  1:2185862/38‑1 (MQ=255)
                                            gACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATgg                                                     >  1:1314347/1‑43 (MQ=255)
                                            gACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATg                                                      <  1:954096/42‑1 (MQ=255)
                                            gACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATg                                                      <  1:1062666/42‑1 (MQ=255)
                                            gACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATg                                                      <  1:2907780/42‑1 (MQ=255)
                                            gACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtct                         <  1:1695461/71‑1 (MQ=255)
                                            gACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtct                         <  1:1066680/71‑1 (MQ=255)
                                             aCCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtctc                        <  1:14805/71‑1 (MQ=255)
                                               cGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtctctc                      <  1:1882145/71‑1 (MQ=255)
                                                  gtTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAt                           >  1:2504553/1‑63 (MQ=255)
                                                   ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                   >  1:1370585/1‑70 (MQ=255)
                                                   ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                   >  1:3640714/1‑70 (MQ=255)
                                                   ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                   >  1:1338445/1‑70 (MQ=255)
                                                   ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCAAGCTATCTCTCGCt                   >  1:3011493/1‑70 (MQ=255)
                                                    tGCCGGTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTgc                 >  1:1516964/1‑70 (MQ=255)
                                                     gCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGAt                >  1:3524794/1‑71 (MQ=255)
                                                      ccggTTTCTCGCCGATCCCGGCTATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCtt               <  1:1692044/69‑1 (MQ=255)
                                                      cTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGAtt               <  1:276185/71‑1 (MQ=255)
                                                      cTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGAt                <  1:3205353/70‑1 (MQ=255)
                                                       cggTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTg              <  1:1763539/70‑1 (MQ=255)
                                                         gTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGc            <  1:1921210/71‑1 (MQ=255)
                                                         gTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                 >  1:688888/1‑66 (MQ=255)
                                                         gTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                 >  1:1134028/1‑66 (MQ=255)
                                                         gTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGc            <  1:3392997/71‑1 (MQ=255)
                                                          tttCTCGCCAATTCCGGCAATGTCGATGGTTTCTTa                                              >  1:118393/1‑36 (MQ=255)
                                                           ttCTCGCCGATCCCGGCGATGTCGATGGTTTCTTa                                              >  1:2212534/1‑35 (MQ=25)
                                                            tctcGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGATATCTCTCGCTGATTGGCgg          >  1:503240/1‑70 (MQ=255)
                                                             ctcGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTa        >  1:34490/1‑71 (MQ=255)
                                                                gCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCAcct     <  1:2778830/71‑3 (MQ=255)
                                                                gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTg                                   <  1:609973/41‑1 (MQ=255)
                                                                gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTa                            >  1:3001070/1‑48 (MQ=255)
                                                                gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                 >  1:754515/1‑59 (MQ=255)
                                                                gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTAc       >  1:3864520/1‑69 (MQ=255)
                                                                gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTAc       >  1:1179512/1‑69 (MQ=255)
                                                                gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACtt     <  1:604534/71‑1 (MQ=255)
                                                                gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACtt     <  1:3831561/71‑1 (MQ=255)
                                                                gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACtt     <  1:741763/71‑1 (MQ=255)
                                                                gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACtt     >  1:2991987/1‑71 (MQ=255)
                                                                 ccAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTg    >  1:2483212/1‑71 (MQ=255)
                                                                  cAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGc   >  1:3388541/1‑71 (MQ=255)
                                                                  cAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGc   <  1:3323066/71‑1 (MQ=255)
                                                                   aaTTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTgg             >  1:3641737/1‑60 (MQ=255)
                                                                   aaTTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCt  >  1:1457850/1‑71 (MQ=255)
                                                                   |                                                                      
TGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCT  >  minE/832048‑832185

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: