Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,670 C→A 7.9% V376V (GTC→GTA narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,6700CA 7.9% 194.6 / 3.1 63V376V (GTC→GTAnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (31/26);  new base A (2/3);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

GCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGC  >  minE/832607‑832730
                                                               |                                                            
gctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCttt                                                        <  1:2230023/70‑1 (MQ=255)
 ctgcGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTaa                                                      >  1:3805568/1‑71 (MQ=255)
   gcgcgcTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATa                                                    >  1:4610027/1‑71 (MQ=255)
    cgcgcTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATAc                                                   >  1:335176/1‑71 (MQ=255)
    cgcgcTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATAc                                                   >  1:3926840/1‑71 (MQ=255)
     gcgcTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATAcc                                                  <  1:3178802/71‑1 (MQ=255)
         tgctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCaa                                                >  1:4906603/1‑69 (MQ=255)
         tgctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATg                                              >  1:390499/1‑71 (MQ=255)
           ctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGcc                                                           >  1:182777/1‑56 (MQ=255)
           ctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGc                                            >  1:5028076/1‑71 (MQ=255)
           ctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGc                                            >  1:475166/1‑71 (MQ=255)
                  ggTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATgg                                     >  1:1609186/1‑71 (MQ=255)
                   gTTGATGCGCTGGGCCAGTAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAAt                                       <  1:2994825/68‑1 (MQ=255)
                   gTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTTGCCtt                                                         <  1:4291987/50‑1 (MQ=255)
                   gTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCtt                                                         <  1:3653463/50‑1 (MQ=255)
                   gTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCtt                                                         <  1:3756786/50‑1 (MQ=255)
                   gTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATg                                              >  1:4598818/1‑61 (MQ=255)
                   gTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGaaa                                        <  1:3994189/67‑1 (MQ=255)
                   gTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGaaa                                        <  1:4205509/67‑1 (MQ=255)
                   gTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAAt                                       <  1:3972065/68‑1 (MQ=255)
                   gTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAAt                                       <  1:2212994/68‑1 (MQ=255)
                   gTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATgg                                     <  1:4221993/70‑1 (MQ=255)
                   gTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATg                                      <  1:2192330/69‑1 (MQ=255)
                      gATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCg                                           >  1:3111252/1‑61 (MQ=255)
                       aTGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTgcg                                >  1:1871332/1‑71 (MQ=255)
                          cgcTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATa                                                    >  1:2382110/1‑48 (MQ=255)
                           gcTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCg                             >  1:2032484/1‑70 (MQ=255)
                               ggCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACgg                        >  1:2288422/1‑71 (MQ=255)
                                gCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCtt                                                         <  1:3616737/37‑1 (MQ=255)
                                   aGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCg                             <  1:1852836/62‑1 (MQ=255)
                                   aGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCt                    <  1:486419/71‑1 (MQ=255)
                                   aGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACAAATGGCGAAATGGTTGCGCCGAACGGTTc                     >  1:1381824/1‑70 (MQ=255)
                                    ggAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGt                                    >  1:797062/1‑54 (MQ=255)
                                    ggAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTgc                                 <  1:612997/57‑1 (MQ=255)
                                     gAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGt                                    >  1:287084/1‑53 (MQ=255)
                                      aaaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTgcgc                               >  1:781686/1‑57 (MQ=255)
                                      aaaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTgcgc                               >  1:2697215/1‑57 (MQ=255)
                                      aaaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTAtt                 >  1:497840/1‑71 (MQ=255)
                                       aaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTAt                  >  1:5661294/1‑69 (MQ=255)
                                       aaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTAt                  >  1:5679711/1‑69 (MQ=255)
                                       aaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTg                >  1:4676397/1‑71 (MQ=255)
                                       aaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTg                >  1:3755832/1‑71 (MQ=255)
                                        aaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTgg               >  1:309927/1‑71 (MQ=255)
                                        aaCAATCCGCAGTGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGAtcg                >  1:1829555/1‑68 (MQ=255)
                                        aaCAATCCGCAGTGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGAtc                 >  1:1287352/1‑68 (MQ=255)
                                             tCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGt                                    <  1:690011/45‑1 (MQ=255)
                                                ggAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACgg                        <  1:2235080/54‑1 (MQ=255)
                                                ggAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACgg                        <  1:222466/54‑1 (MQ=255)
                                                ggAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACgg                        <  1:1353414/54‑1 (MQ=255)
                                                ggAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCt       <  1:3559558/71‑1 (MQ=255)
                                                  agtGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTgg     <  1:2090708/69‑1 (MQ=255)
                                                  agtGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTgg     <  1:5462661/69‑1 (MQ=255)
                                                  agtGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTgg     <  1:4568169/69‑1 (MQ=255)
                                                  aaTGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTgg     <  1:1443698/71‑1 (MQ=255)
                                                   aTGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCt                    <  1:2674713/55‑1 (MQ=255)
                                                    tGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAAc                          <  1:1467036/48‑1 (MQ=255)
                                                    tGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTAtt                 <  1:3854500/57‑1 (MQ=255)
                                                     ggAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTc                     <  1:5639540/52‑1 (MQ=255)
                                                      gAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTgg               >  1:4543940/1‑57 (MQ=255)
                                                      gAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGc  >  1:3160069/1‑70 (MQ=255)
                                                        aaaCTGTCGCCTTTAATACCAATGGCGAAATGGtt                                   >  1:2026243/1‑35 (MQ=255)
                                                         aaCTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTggg    <  1:4571973/65‑1 (MQ=255)
                                                          aCTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGc        >  1:2670162/1‑60 (MQ=255)
                                                               |                                                            
GCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGC  >  minE/832607‑832730

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