Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,418,334 G→A 41.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,3340GA41.8% 48.3 / 97.3 79V95I (GTC→ATC) narPDNA‑binding response regulator in two‑component regulatory system with NarQ or NarX
Reads supporting (aligned to +/- strand):  ref base G (16/30);  new base A (20/13);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 3.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

TTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTG  >  minE/1418271‑1418402
                                                               |                                                                    
ttGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTc                                                                    <  1:2851162/66‑1 (MQ=255)
ttGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTc                                                                    <  1:4069200/66‑1 (MQ=255)
ttGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTc                                                                    <  1:1464453/66‑1 (MQ=255)
 tGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTc                                                                    <  1:532225/65‑1 (MQ=255)
  gcgcAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCTATGTCTTTgc                                                               <  1:2251669/69‑1 (MQ=255)
  gcgcAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTgcgc                                                             >  1:3780646/1‑71 (MQ=255)
  gcgcAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTgc                                                               <  1:2952815/69‑1 (MQ=255)
  gcgcAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTgc                                                               <  1:3120540/69‑1 (MQ=255)
  gcgcAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTgc                                                               <  1:479050/69‑1 (MQ=255)
  gcgcAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTgc                                                               <  1:1816591/69‑1 (MQ=255)
       gggATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATa                                                        <  1:441667/71‑1 (MQ=255)
       gggATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCACTGATa                                                        <  1:2801501/71‑1 (MQ=255)
         gATGGCGTTACCGCGCAAATTATTATCGTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGa                                                      >  1:1182445/1‑71 (MQ=255)
         gATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGa                                                      >  1:4114358/1‑71 (MQ=255)
         gATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGa                                                      >  1:4059716/1‑71 (MQ=255)
         gATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGa                                                      >  1:1224311/1‑71 (MQ=255)
          aTGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTg                                                                >  1:1589783/1‑60 (MQ=255)
          aTGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGAc                                                     >  1:413044/1‑71 (MQ=255)
          aTGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGAc                                                     >  1:2836799/1‑71 (MQ=255)
          aTGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGAc                                                     >  1:2839513/1‑71 (MQ=255)
          aTGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGAc                                                     >  1:794061/1‑71 (MQ=255)
          aTGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGAc                                                     >  1:3843157/1‑71 (MQ=255)
          aTGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGAc                                                     >  1:442333/1‑71 (MQ=255)
          aTGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGAc                                                     >  1:3921525/1‑71 (MQ=255)
             gCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTg                                                           <  1:744252/62‑1 (MQ=255)
             gCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGc                                                   <  1:2566895/70‑1 (MQ=255)
              cGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAg                                                 >  1:2033123/1‑71 (MQ=255)
                ttACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCttt                                                                 <  1:1122304/53‑1 (MQ=255)
                ttACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGAt                                                         <  1:2854284/61‑1 (MQ=255)
                  aCCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGgcg                                              <  1:1283544/70‑1 (MQ=255)
                  aCCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGgcg                                              >  1:3461494/1‑70 (MQ=255)
                  aCCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGgcg                                              >  1:445220/1‑70 (MQ=255)
                  aCCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGgcg                                              >  1:4122240/1‑70 (MQ=255)
                  aCCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGgcg                                              <  1:1531111/70‑1 (MQ=255)
                  aCCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGAAgg                                                >  1:4019608/1‑68 (MQ=255)
                  aCCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATAACTTTGCGCTGATAGACGCAGgcg                                              <  1:1063822/70‑1 (MQ=255)
                  aCCGAGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGgcg                                              >  1:3621797/1‑70 (MQ=255)
                            ttattaTCCTGACCGTATCCGATGCCTCCAGCGATATc                                                                    <  1:781016/38‑1 (MQ=255)
                            ttattaTCCTGACCGTATCCGATGCCTCCAGCGATATc                                                                    <  1:3938932/38‑1 (MQ=255)
                            ttattaTCCTGACCGTATCCGATGCCTCCAGCGATATc                                                                    <  1:1900947/38‑1 (MQ=255)
                            ttattaTCCTGACCGTATCCGATGCCTCCAGCGATATc                                                                    <  1:1523108/38‑1 (MQ=255)
                            ttattaTCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCt                                                            <  1:381346/46‑1 (MQ=255)
                            ttattaTCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCt                                                            <  1:2216913/46‑1 (MQ=255)
                            ttattaTCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCt                                                            <  1:4095690/46‑1 (MQ=255)
                              attatCCTGAACGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATc                                  <  1:106161/70‑1 (MQ=255)
                               ttatCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTg                                <  1:2941352/71‑1 (MQ=255)
                                 atCCTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGgcg                                              <  1:2131441/55‑1 (MQ=255)
                                    cTGACCGTATCCGATGCCTCCAGCGATATCTTTGCGCTGATAGAAGCAGGCGCAGACGGTTATCTGTTGa                            <  1:1323713/70‑1 (MQ=255)
                                       aCCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGa                        <  1:3374849/71‑1 (MQ=255)
                                       aCCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGa                        <  1:829329/71‑1 (MQ=255)
                                       aCCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGa                        <  1:2754096/71‑1 (MQ=255)
                                           tATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGaa                           <  1:2600374/64‑1 (MQ=255)
                                             tCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGtt                                     <  1:55274/52‑1 (MQ=255)
                                               cGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGc                    >  1:1929908/1‑67 (MQ=255)
                                               cGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGc                    >  1:3561514/1‑67 (MQ=255)
                                               cGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGc                    >  1:1478961/1‑67 (MQ=255)
                                               cGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGc                    >  1:2765934/1‑67 (MQ=255)
                                               cGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGAcc                >  1:1107099/1‑71 (MQ=255)
                                               cGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGAcc                >  1:2824864/1‑71 (MQ=255)
                                               cGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGAcc                >  1:2399797/1‑71 (MQ=255)
                                               cGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGAcc                >  1:1537142/1‑71 (MQ=255)
                                               cGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGAcc                >  1:2645725/1‑71 (MQ=255)
                                               cGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGAcc                >  1:2670738/1‑71 (MQ=255)
                                               cGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGAcc                >  1:1386831/1‑71 (MQ=255)
                                               cGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGAcc                >  1:2764047/1‑71 (MQ=255)
                                               cGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGAcc                >  1:1277533/1‑71 (MQ=255)
                                               cGATGCCTCCAGCGATATCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGAcc                >  1:1254502/1‑71 (MQ=255)
                                                gATGCCTCCAGCGATGTCTTTGCGCTGATAGGCGCAGGCGCAGACGGTTATCTGTTGaaa                          <  1:3686685/60‑1 (MQ=255)
                                                gATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGaaa                          <  1:3748865/60‑1 (MQ=255)
                                                  tGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGa                        >  1:3755850/1‑60 (MQ=255)
                                                  tGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGa                        >  1:3475957/1‑60 (MQ=255)
                                                     cTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAg          <  1:3418895/71‑1 (MQ=255)
                                                      tCCAGCGATGTCTTTGCGCTTATAGACGCAGGCGCAGACGGTTATCTGTTGAAAg                         <  1:155829/55‑1 (MQ=255)
                                                      tCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAg                         <  1:954941/55‑1 (MQ=255)
                                                      tCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAg                         <  1:1776858/55‑1 (MQ=255)
                                                      tCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAg                         <  1:2830863/55‑1 (MQ=255)
                                                      tCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACa                      <  1:3716664/58‑1 (MQ=255)
                                                        cAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTAt       >  1:3515924/1‑71 (MQ=255)
                                                             aTGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTg  <  1:2160144/71‑1 (MQ=255)
                                                               |                                                                    
TTGCGCAGGGATGGCGTTACCGCGCAAATTATTATCCTGACCGTATCCGATGCCTCCAGCGATGTCTTTGCGCTGATAGACGCAGGCGCAGACGGTTATCTGTTGAAAGACAGCGACCCGGAAGTATTGCTG  >  minE/1418271‑1418402

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