Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,418,334 G→A 23.8% 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,3340GA23.8% 81.1 / 22.6 42V95I (GTC→ATC) narPDNA‑binding response regulator in two‑component regulatory system with NarQ or NarX
Reads supporting (aligned to +/- strand):  ref base G (15/17);  new base A (2/8);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

ATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGG  >  minE/1418266‑1418403
                                                                    |                                                                     
aTGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTc                                                                     >  1:360774/1‑71 (MQ=255)
aTGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTc                                                                     >  1:929314/1‑71 (MQ=255)
 tGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCt                                                                    <  1:114017/71‑1 (MQ=255)
   ccTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCTATGCCTCCAGCGATGTCtt                                                                   >  1:510995/1‑70 (MQ=255)
   ccTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCtt                                                                   <  1:1963660/70‑1 (MQ=255)
   ccTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCtt                                                                   >  1:1602486/1‑70 (MQ=255)
   ccTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCtt                                                                   <  1:1544162/70‑1 (MQ=255)
      tGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCttt                                                                  >  1:2142081/1‑68 (MQ=255)
        cgcAGGTATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCt                                                             <  1:998006/71‑1 (MQ=255)
        cgcAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCt                                                             <  1:1085541/71‑1 (MQ=255)
            gggATGGCGTTACCTCGCAAATTATTATCATGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTgt                                                           >  1:383472/1‑68 (MQ=255)
            gggATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGa                                                           >  1:890314/1‑69 (MQ=255)
            gggATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATa                                                         <  1:1652041/71‑1 (MQ=255)
              gATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATg                                                                       >  1:353256/1‑55 (MQ=255)
              gATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGAGATg                                                                       >  1:622624/1‑55 (MQ=255)
                     ttACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGt                                                                      <  1:386505/49‑1 (MQ=255)
                      tACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGt                                                                      <  1:1203313/48‑1 (MQ=255)
                      tACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGgc                                                <  1:292461/70‑1 (MQ=255)
                       aCCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGgcg                                               >  1:1074081/1‑70 (MQ=255)
                        ccGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCa                                             <  1:435796/71‑1 (MQ=255)
                        ccGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCa                                             <  1:1562304/71‑1 (MQ=255)
                        ccGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCa                                             <  1:505867/71‑1 (MQ=255)
                        ccGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCa                                             <  1:810108/71‑1 (MQ=255)
                        ccGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCa                                             <  1:2070515/71‑1 (MQ=255)
                                attattATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGt                                       >  1:445529/1‑69 (MQ=255)
                                  tattatCCTGACCGTATCCGATGCCTCCAGCGATGt                                                                      <  1:481336/36‑1 (MQ=255)
                                  tattatCCTGACCGTATCCGATGCCTCCAGCGATGt                                                                      <  1:2098322/36‑1 (MQ=255)
                                  tattatCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATa                                                         <  1:1603150/49‑1 (MQ=255)
                                      atCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGt                                >  1:881412/1‑70 (MQ=255)
                                         cTGACCGTATACGATGCCTCCAGCGATGTCTTTGCGCTGATAGaa                                                      >  1:1962298/1‑44 (MQ=255)
                                            aCCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGaaa                           <  1:1936281/69‑1 (MQ=255)
                                            aCCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAG‑‑GCAGACGGTTATCTGTTGAAAGAc                        <  1:1981095/70‑1 (MQ=255)
                                             ccGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGAc                        <  1:916907/71‑1 (MQ=255)
                                                    cGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGAcc                 >  1:1996174/1‑71 (MQ=255)
                                                    cGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGAcc                 >  1:870251/1‑71 (MQ=255)
                                                    cGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGAcc                 >  1:1401534/1‑71 (MQ=255)
                                                          cTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGaaa                           <  1:2049019/55‑1 (MQ=255)
                                                          cTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGaaa                           <  1:1764144/55‑1 (MQ=255)
                                                             cAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGa                         >  1:950455/1‑54 (MQ=255)
                                                             cAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGa                         >  1:1661268/1‑54 (MQ=255)
                                                             cAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCg                    <  1:1292934/59‑1 (MQ=255)
                                                               gCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTg      <  1:1943164/71‑1 (MQ=255)
                                                               gCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTg      <  1:335887/71‑1 (MQ=255)
                                                                  aTGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTg                              >  1:179943/1‑44 (MQ=255)
                                                                    gTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGaa                            >  1:356574/1‑44 (MQ=255)
                                                                    gTCTTTGCGCTGATAGAAGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTgg  >  1:1075004/1‑70 (MQ=255)
                                                                    |                                                                     
ATGCCTTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTGG  >  minE/1418266‑1418403

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: