Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 247,251 A→G 8.3% K298K (AAA→AAG phoR → sensory histidine kinase in two‑component regulatory system with PhoB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE247,2510AG 8.3% 159.0 / 7.1 60K298K (AAA→AAGphoRsensory histidine kinase in two‑component regulatory system with PhoB
Reads supporting (aligned to +/- strand):  ref base A (27/28);  new base G (3/2);  total (30/30)
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.97e-01

TGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTG  >  minE/247182‑247318
                                                                     |                                                                   
tGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGaaaaaa                                                                     >  1:2487335/1‑70 (MQ=255)
    aaGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGAc                                                                <  1:1306295/71‑1 (MQ=255)
         tGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAGCAGACATTTa                                                           >  1:2646430/1‑71 (MQ=255)
         tGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAGCAGACATTTa                                                           >  1:874740/1‑71 (MQ=255)
         tGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAGCAGACATTTa                                                           >  1:3906309/1‑71 (MQ=255)
              tgCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCttt                                                      >  1:2036274/1‑71 (MQ=255)
                     gatgCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTgaga                                                  >  1:1159806/1‑68 (MQ=255)
                     gatgCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTgaga                                                  >  1:175076/1‑68 (MQ=255)
                       tgctgcGCGTTGTTGAGCGCTAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGagatagat                                             <  1:2566518/71‑1 (MQ=255)
                       tgctgcGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGagatagat                                             <  1:969077/71‑1 (MQ=255)
                       tgctgcGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGagatagat                                             <  1:3249163/71‑1 (MQ=255)
                       tgctgcGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGagatagat                                             <  1:2302473/71‑1 (MQ=255)
                       tgctgcGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGagatagat                                             <  1:188984/71‑1 (MQ=255)
                       tgctgcGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGagatagat                                             <  1:4132240/71‑1 (MQ=255)
                        gctgcGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAGAAACAGACATTTa                                                           >  1:916852/1‑56 (MQ=255)
                        gctgcGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTa                                                           >  1:1717463/1‑56 (MQ=255)
                        gctgcGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTa                                                           >  1:3562333/1‑56 (MQ=255)
                        gctgcGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTa                                                           >  1:3714705/1‑56 (MQ=255)
                        gctgcGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTa                                                           >  1:2370018/1‑56 (MQ=255)
                        gctgcGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGagatagat                                             <  1:4059792/70‑1 (MQ=255)
                            cgcgTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACgg                                        >  1:2372788/1‑71 (MQ=255)
                                 tgttgAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTc                                    >  1:413026/1‑70 (MQ=255)
                                  gttgAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTc                                    >  1:2379067/1‑69 (MQ=255)
                                  gttgAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTc                                    >  1:1315369/1‑69 (MQ=255)
                                  gttgAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTc                                    >  1:3845404/1‑69 (MQ=255)
                                  gttgAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTc                                    >  1:242869/1‑69 (MQ=255)
                                  gttgAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTc                                    >  1:604630/1‑69 (MQ=255)
                                  gttgAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTc                                    >  1:539600/1‑69 (MQ=255)
                                  gttgAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTc                                    >  1:4110309/1‑69 (MQ=255)
                                  gttgAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCa                                   <  1:922317/70‑1 (MQ=255)
                                   ttgAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCaa                                  <  1:486371/70‑1 (MQ=255)
                                   ttgAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCaa                                  <  1:3148149/70‑1 (MQ=255)
                                          cgAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTg                          <  1:1884124/71‑1 (MQ=255)
                                                 cAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACgc                                        >  1:1124341/1‑49 (MQ=255)
                                                 cAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGgt                               >  1:4021464/1‑59 (MQ=255)
                                                 cAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGgt                               >  1:3440168/1‑59 (MQ=255)
                                                 cAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGgt                               >  1:3649081/1‑59 (MQ=255)
                                                 cAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGc                        <  1:2992110/66‑1 (MQ=255)
                                                 cAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGc                        <  1:1121687/66‑1 (MQ=255)
                                                 cAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGc                        <  1:1001642/66‑1 (MQ=255)
                                                 cAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACg                    >  1:3979656/1‑70 (MQ=255)
                                                 cAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACg                    >  1:1338295/1‑70 (MQ=255)
                                                  aGACTCTGATTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACg                    <  1:3394844/69‑1 (MQ=255)
                                                  aGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACg                    <  1:1117152/69‑1 (MQ=255)
                                                     ctctGAGTCAGAAAAAACAGACATTTACCTTTgaga                                                  >  1:917622/1‑36 (MQ=255)
                                                     ctctGAGTCAGAAAAAACAGACATTTACCTTTgaga                                                  >  1:3907713/1‑36 (MQ=255)
                                                     ctctGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGAt               <  1:3650140/71‑1 (MQ=255)
                                                     ctctGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGAt               <  1:1815639/71‑1 (MQ=255)
                                                     ctctGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGAt               <  1:241572/71‑1 (MQ=255)
                                                     ctctGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGAt               <  1:2655247/71‑1 (MQ=255)
                                                     ctctGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGAt               <  1:2109700/71‑1 (MQ=255)
                                                      tctGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGTCAACGAAg                 <  1:260852/68‑1 (MQ=255)
                                                      tctGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAg                 <  1:3379399/68‑1 (MQ=255)
                                                      tctGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAg                 <  1:1156307/68‑1 (MQ=255)
                                                      tctGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAg                 <  1:2512006/68‑1 (MQ=255)
                                                            tCAGAAAAAGCAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTAc        <  1:1563336/71‑1 (MQ=255)
                                                            tCAGAAAAAGCAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTAc        <  1:2629461/71‑1 (MQ=255)
                                                            tCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGa                >  1:350590/1‑63 (MQ=255)
                                                              aGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTc                            >  1:608824/1‑49 (MQ=255)
                                                               gAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCa     <  1:2468961/71‑1 (MQ=255)
                                                               gAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCa     <  1:846426/71‑1 (MQ=255)
                                                                  aaaaCAGACATTTGCCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGt   >  1:4201026/1‑70 (MQ=255)
                                                                  aaaaCAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTTTGGCAACGAAGa                <  1:3327024/57‑1 (MQ=255)
                                                                  aaaaCAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGa                <  1:1926851/57‑1 (MQ=255)
                                                                  aaaaCAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGa                <  1:146657/57‑1 (MQ=255)
                                                                  aaaaCAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGa                <  1:3716046/57‑1 (MQ=255)
                                                                  aaaaCAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGa                <  1:919715/57‑1 (MQ=255)
                                                                  aaaaCAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTg  >  1:3862970/1‑71 (MQ=255)
                                                                  aaaaCAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTg  >  1:3015755/1‑71 (MQ=255)
                                                                     aCAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAg            >  1:4081446/1‑58 (MQ=255)
                                                                     |                                                                   
TGAAAAGGTTGATGTGCCGATGATGCTGCGCGTTGTTGAGCGCGAGGCTCAGACTCTGAGTCAGAAAAAACAGACATTTACCTTTGAGATAGATAACGGCCTCAAGGTGTCTGGCAACGAAGATCAGCTACGCAGTG  >  minE/247182‑247318

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: