Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,684 A→C 7.8% N381T (AAT→ACT) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,6840AC 7.8% 157.0 / 7.3 64N381T (AAT→ACT) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base C (2/3);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

TGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTG  >  minE/832616‑832752
                                                                    |                                                                    
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)
          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)
                  gcTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCg                                                   >  1:2032484/1‑70 (MQ=255)
                      ggCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACgg                                              >  1:2288422/1‑71 (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: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: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:5462661/69‑1 (MQ=255)
                                         agtGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTgg                           <  1:2090708/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)
                                                       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)
                                                            tAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGc                              >  1:1573237/1‑49 (MQ=255)
                                                            tAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGc                              >  1:4182093/1‑49 (MQ=255)
                                                            tAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGc                              >  1:5345206/1‑49 (MQ=255)
                                                            tAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATgg         >  1:4799910/1‑70 (MQ=255)
                                                               tACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCg                       <  1:1178768/53‑1 (MQ=255)
                                                               tACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATc     <  1:4637853/71‑1 (MQ=255)
                                                                aCCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCg                       <  1:2057519/52‑1 (MQ=255)
                                                                  cAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCtt   <  1:1573428/70‑1 (MQ=255)
                                                                  cAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTg  <  1:1672421/71‑1 (MQ=255)
                                                                  cAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTg  <  1:4891240/71‑1 (MQ=255)
                                                                    |                                                                    
TGCTGATCTGGTTGATGCGCTGGGCCAGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTG  >  minE/832616‑832752

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: