Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,103 T→C 15.9% A187A (GCT→GCC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,1030TC15.9% 146.5 / 19.3 63A187A (GCT→GCCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (27/26);  new base C (3/7);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATT  >  minE/832035‑832172
                                                                    |                                                                     
gaggTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTgg                                                                     <  1:2809470/71‑1 (MQ=255)
   gTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCTGACCGTGTTGCTGGt                                                                    <  1:3855664/69‑1 (MQ=255)
     gaacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTgg                                                                     >  1:788178/1‑66 (MQ=255)
     gaacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGttt                                                                  <  1:351880/69‑1 (MQ=255)
             tGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAAtt                                                        >  1:1041434/1‑71 (MQ=255)
                tCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGttt                                                                  <  1:2612761/58‑1 (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:636593/69‑1 (MQ=255)
                       cgcTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGt                                              <  1:3671646/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)
                                             aaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGcc                                                            >  1:1054923/1‑35 (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:2791558/1‑46 (MQ=255)
                                                   ggCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGa                                           >  1:3263008/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:1062666/42‑1 (MQ=255)
                                                         gACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATg                                         <  1:2907780/42‑1 (MQ=255)
                                                         gACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATg                                         <  1:954096/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:1338445/1‑70 (MQ=255)
                                                                ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt      >  1:3640714/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)
                                                                   cTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGAtt  <  1:276185/71‑1 (MQ=255)
                                                                   cTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGAt   <  1:3205353/70‑1 (MQ=255)
                                                                    |                                                                     
GAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATT  >  minE/832035‑832172

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: