Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,118 T→C 18.5% I192I (ATT→ATC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,1180TC18.5% 172.8 / 32.8 81I192I (ATT→ATCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (38/28);  new base C (4/11);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 4.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

TGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTT  >  minE/832048‑832186
                                                                      |                                                                    
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:2802201/1‑47 (MQ=255)
                               aaaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCaa                                                               >  1:2730605/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:616363/1‑71 (MQ=255)
                                           ggACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtc                           >  1:3637623/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:1066680/71‑1 (MQ=255)
                                            gACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtct                          <  1:1695461/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:3640714/1‑70 (MQ=255)
                                                   ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                    >  1:1338445/1‑70 (MQ=255)
                                                   ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                    >  1:1370585/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:1134028/1‑66 (MQ=255)
                                                         gTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                  >  1:688888/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:1179512/1‑69 (MQ=255)
                                                                gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTAc        >  1:3864520/1‑69 (MQ=255)
                                                                gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACtt      <  1:3831561/71‑1 (MQ=255)
                                                                gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACtt      <  1:604534/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:3323066/71‑1 (MQ=255)
                                                                  cAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGc    >  1:3388541/1‑71 (MQ=255)
                                                                   aaTTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTgg              >  1:3641737/1‑60 (MQ=255)
                                                                   aaTTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCt   >  1:1457850/1‑71 (MQ=255)
                                                                     ttCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCt   >  1:3323626/1‑69 (MQ=255)
                                                                     ttCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGAGGTACTTGCtt  >  1:3808029/1‑70 (MQ=255)
                                                                      tCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAt                            >  1:1182559/1‑43 (MQ=255)
                                                                      tCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCtt  >  1:2817410/1‑69 (MQ=255)
                                                                      |                                                                    
TGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTT  >  minE/832048‑832186

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: