Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,418,334 G→A 66.2% V95I (GTC→ATC)  narP → DNA‑binding response regulator in two‑component regulatory system with NarQ or NarX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,418,3340GA66.2% 73.9 / 67.3 68V95I (GTC→ATC) narPDNA‑binding response regulator in two‑component regulatory system with NarQ or NarX
Reads supporting (aligned to +/- strand):  ref base G (10/13);  new base A (28/17);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 1.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTAT  >  minE/1418268‑1418397
                                                                  |                                                               
gCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCtt                                                             >  1:2464433/1‑71 (MQ=255)
gCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCtt                                                             >  1:2034379/1‑71 (MQ=255)
 ccTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCttt                                                            >  1:4051155/1‑71 (MQ=255)
 ccTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCttt                                                            >  1:2613106/1‑71 (MQ=255)
      cgcAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCt                                                       <  1:3514373/71‑1 (MQ=255)
        cAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGa                                                     <  1:1521351/71‑1 (MQ=255)
        cAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGa                                                     <  1:1837685/71‑1 (MQ=255)
        cAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGa                                                     <  1:3429298/71‑1 (MQ=255)
        cAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGa                                                     <  1:2906648/71‑1 (MQ=255)
             aTGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGa                                                 >  1:815351/1‑70 (MQ=255)
             aTGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGAc                                                >  1:681350/1‑71 (MQ=255)
              tGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCTATGCCTCCAGCGATGTCTTTgc                                                          <  1:1291870/60‑1 (MQ=255)
              tGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTgc                                                          <  1:612367/60‑1 (MQ=255)
              tGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTgc                                                          <  1:1373835/60‑1 (MQ=255)
              tGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATa                                                   <  1:39716/67‑1 (MQ=255)
              tGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATa                                                   <  1:1361339/67‑1 (MQ=255)
                     aCCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGgcg                                         >  1:967295/1‑70 (MQ=255)
                            aaaTTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGt                                 <  1:2101787/71‑1 (MQ=255)
                             aaTTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCt                                                       <  1:1243079/48‑1 (MQ=255)
                             aaTTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCt                                                       <  1:4054441/48‑1 (MQ=255)
                             aaTTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCt                                                       <  1:3811547/48‑1 (MQ=255)
                             aaTTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCt                                                       <  1:2861269/48‑1 (MQ=255)
                             aaTTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGtt                                <  1:4178642/71‑1 (MQ=255)
                              attattATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTa                               >  1:1818617/1‑71 (MQ=255)
                              attattATCCTGACCGTATCCGATGCCTACAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTa                               >  1:1682649/1‑71 (MQ=255)
                               ttattaTCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTAt                              <  1:1799774/71‑1 (MQ=255)
                               ttattaTCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTAt                              <  1:3612320/71‑1 (MQ=255)
                               ttattaTCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTa                               >  1:2438688/1‑70 (MQ=255)
                               ttattaTCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTa                               >  1:2667875/1‑70 (MQ=255)
                               ttattaTCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTAt                              >  1:3614735/1‑71 (MQ=255)
                               ttattaTCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTAt                              >  1:3764975/1‑71 (MQ=255)
                               ttattaTCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTAt                              >  1:3637691/1‑71 (MQ=255)
                               ttattaTCCTGACCGTATCCGATGCCTACAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTAt                              >  1:3620620/1‑71 (MQ=255)
                                tattatCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTAt                              <  1:2777949/70‑1 (MQ=255)
                                tattatCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGa                                                     >  1:711203/1‑47 (MQ=255)
                                tattatCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGa                                                     >  1:3863178/1‑47 (MQ=255)
                                tattatCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGa                                                     >  1:418569/1‑47 (MQ=255)
                                tattatCCTGACCGTATCCGATGACTCCAGCGATATCTTTGCGCTGa                                                     >  1:4007153/1‑47 (MQ=39)
                                  ttatCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGa                                                     >  1:1914448/1‑45 (MQ=255)
                                  ttatCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGa                                                     >  1:372287/1‑45 (MQ=255)
                                     tCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGa                                                     <  1:1842521/42‑1 (MQ=255)
                                      ccTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTg                        >  1:279848/1‑70 (MQ=255)
                                          aCCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGa                                                     >  1:2527629/1‑37 (MQ=255)
                                          aCCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGa                                                     >  1:2022970/1‑37 (MQ=255)
                                          aCCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGa                                                     >  1:4070007/1‑37 (MQ=255)
                                          aCCGTATCCGATGCCTCCAGCGATATCTTTTCGCTGATAGACGCAGGCGCAGACg                                   >  1:3960116/1‑55 (MQ=39)
                                          aCCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGaaa                     <  1:2384924/69‑1 (MQ=255)
                                          aCCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGa                   >  1:2327021/1‑71 (MQ=255)
                                            cGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACa                 <  1:3015172/71‑1 (MQ=255)
                                             gTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGtt                                <  1:2934232/55‑1 (MQ=255)
                                             gTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGtt                                <  1:4219478/55‑1 (MQ=255)
                                             gTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGtt                                <  1:3326493/55‑1 (MQ=255)
                                              tATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGc               <  1:1637693/71‑1 (MQ=255)
                                              tATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGtt                                <  1:2111969/54‑1 (MQ=255)
                                                tCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGa             >  1:1873940/1‑71 (MQ=255)
                                                 ccGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGa                       <  1:828771/60‑1 (MQ=255)
                                                  cGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGAcc           >  1:2427108/1‑71 (MQ=255)
                                                   gATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGaaa                     >  1:2758204/1‑60 (MQ=255)
                                                   gATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGaaa                     >  1:702211/1‑60 (MQ=255)
                                                   gATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGaaa                     >  1:3026909/1‑60 (MQ=255)
                                                   gATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGaaa                     >  1:2900662/1‑60 (MQ=255)
                                                      gCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGAccc          >  1:3816148/1‑68 (MQ=255)
                                                        cTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAg     <  1:32207/71‑1 (MQ=255)
                                                          ccAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGtt                         >  1:210840/1‑49 (MQ=255)
                                                           cAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTg                           >  1:498073/1‑46 (MQ=255)
                                                           cAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTg                           >  1:1964232/1‑46 (MQ=255)
                                                           cAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTAt  >  1:3079359/1‑71 (MQ=255)
                                                           cAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTAt  <  1:4072277/71‑1 (MQ=255)
                                                                  |                                                               
GCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTAT  >  minE/1418268‑1418397

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: