Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,674 T→G 8.2% F378V (TTT→GTT)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,6740TG 8.2% 151.1 / 7.4 61F378V (TTT→GTT) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (32/24);  new base G (2/3);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAA  >  minE/832607‑832741
                                                                   |                                                                   
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)
           ctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGc                                                       >  1:475166/1‑71 (MQ=255)
           ctgATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGc                                                       >  1:5028076/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:4205509/67‑1 (MQ=255)
                   gTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGaaa                                                   <  1:3994189/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:2697215/1‑57 (MQ=255)
                                      aaaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTgcgc                                          >  1:781686/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:3755832/1‑71 (MQ=255)
                                       aaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTg                           >  1:4676397/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:1353414/54‑1 (MQ=255)
                                                ggAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACgg                                   <  1:222466/54‑1 (MQ=255)
                                                ggAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACgg                                   <  1:2235080/54‑1 (MQ=255)
                                                ggAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCt                  <  1:3559558/71‑1 (MQ=255)
                                                  agtGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTgg                <  1:4568169/69‑1 (MQ=255)
                                                  agtGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTgg                <  1:2090708/69‑1 (MQ=255)
                                                  agtGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTgg                <  1:5462661/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)
                                                                gCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCaa  <  1:3167222/71‑1 (MQ=255)
                                                                gCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCa   >  1:2257604/1‑70 (MQ=255)
                                                                gCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCa   >  1:3800650/1‑70 (MQ=255)
                                                                 ccTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTc                                <  1:888479/40‑1 (MQ=255)
                                                                   |                                                                   
GCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAA  >  minE/832607‑832741

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: