Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 247,122 A→G 38.6% G255G (GGA→GGG 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,1220AG38.6% 42.0 / 58.5 88G255G (GGA→GGGphoRsensory histidine kinase in two‑component regulatory system with PhoB
Reads supporting (aligned to +/- strand):  ref base A (32/22);  new base G (19/15);  total (51/37)
Fisher's exact test for biased strand distribution p-value = 8.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

GGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCTCAATGAAAAG  >  minE/247063‑247188
                                                           |                                                                  
ggCGCGGTACGTGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGGTTGGTGAAg                                                           <  1:153763/69‑1 (MQ=255)
ggCGCGGTACGTGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGGTTGGTGAAg                                                           <  1:2203558/69‑1 (MQ=255)
ggCGCGGTACGTGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGGTTGGTGAAg                                                           <  1:1727722/69‑1 (MQ=255)
ggCGCGGTACGTGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGGTTGGTGAAg                                                           <  1:582741/69‑1 (MQ=255)
ggCGCGGTACGTGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGGTTGGTGAAg                                                           <  1:3176013/69‑1 (MQ=255)
ggCGCGGTACGTGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGGTTGGTGAAg                                                           <  1:218869/69‑1 (MQ=255)
ggCGCGGTACGTGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGGTTGGTGAAg                                                           <  1:1551202/69‑1 (MQ=255)
ggCGCGGTACGTGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGGTTGGTGAAg                                                           <  1:321733/69‑1 (MQ=255)
ggCGCGGTACGTGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGGTTGGTGAAg                                                           <  1:2407390/69‑1 (MQ=255)
ggCGCGGTACGTGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGGTTGGTGAAg                                                           <  1:2827922/69‑1 (MQ=255)
ggCGCGGTACGTGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGGTTGGTGAAg                                                           <  1:2814/69‑1 (MQ=255)
ggCGCGGTACGTGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGGTTGGTGAAg                                                           <  1:765335/69‑1 (MQ=255)
ggCGCGGTACGTGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGGTTGGTGAAg                                                           <  1:1358755/69‑1 (MQ=255)
ggCGCGGTACGTGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGGTTGGTGAAg                                                           <  1:2783829/69‑1 (MQ=255)
ggCGCGGTACGTGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGGTTGGTGAAg                                                           <  1:3372048/69‑1 (MQ=255)
 gcgcGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGc                                                          >  1:854417/1‑69 (MQ=255)
 gcgcGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGc                                                          >  1:70872/1‑69 (MQ=255)
 gcgcGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGc                                                          <  1:1234323/69‑1 (MQ=255)
  cgcgGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGTAAGGACTGGTGAAGCa                                                         <  1:1292887/69‑1 (MQ=255)
  cgcgGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCa                                                         >  1:2920086/1‑69 (MQ=255)
  cgcgGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCa                                                         >  1:2204727/1‑69 (MQ=255)
  cgcgGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCa                                                         <  1:3043479/69‑1 (MQ=255)
  cgcgGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCa                                                         <  1:454124/69‑1 (MQ=255)
  cgcgGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCa                                                         >  1:27221/1‑69 (MQ=255)
  cgcgGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCa                                                         <  1:2006652/69‑1 (MQ=255)
  cgcgGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCa                                                         <  1:390234/69‑1 (MQ=255)
  cgcgGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCa                                                         <  1:875319/69‑1 (MQ=255)
  cgcgGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCa                                                         >  1:1459577/1‑69 (MQ=255)
  cgcgGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCa                                                         >  1:549984/1‑69 (MQ=255)
  cgcgGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCa                                                         >  1:1603297/1‑69 (MQ=255)
  cgcgGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCa                                                         >  1:62295/1‑69 (MQ=255)
  cgcgGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCa                                                         <  1:1537113/69‑1 (MQ=255)
  cgcgGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCa                                                         >  1:723100/1‑69 (MQ=255)
   gcgGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCa                                                         <  1:741343/68‑1 (MQ=255)
                    gTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGaa                                       <  1:2856434/69‑1 (MQ=255)
                             cATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTc                                          >  1:2007464/1‑57 (MQ=255)
                             cATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGc                              >  1:91767/1‑69 (MQ=255)
                             cATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGc                              >  1:717355/1‑69 (MQ=255)
                             cATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGc                              >  1:3175498/1‑69 (MQ=255)
                             cATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGc                              >  1:937745/1‑69 (MQ=255)
                             cATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGc                              >  1:3188503/1‑69 (MQ=255)
                             cATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGc                              >  1:33359/1‑69 (MQ=255)
                             cATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGc                              >  1:1580496/1‑69 (MQ=255)
                             cATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGc                              >  1:154124/1‑69 (MQ=255)
                             cATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGc                              >  1:2037781/1‑69 (MQ=255)
                             cATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGc                              >  1:2443665/1‑69 (MQ=255)
                             cATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGc                              >  1:2130987/1‑69 (MQ=255)
                             cATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGc                              >  1:2442808/1‑69 (MQ=255)
                             cATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGc                              >  1:1540857/1‑69 (MQ=255)
                             cATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGc                              >  1:2289204/1‑69 (MQ=255)
                             cATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGc                              >  1:1124279/1‑69 (MQ=255)
                             cATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGc                              >  1:1601438/1‑69 (MQ=255)
                             cATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTCTGACGCTGTCGAAAATAGAAGc                              >  1:1590969/1‑68 (MQ=255)
                                      gCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGa                      <  1:1899049/68‑1 (MQ=255)
                                        aGACCCAGCGGATGGAATGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGc                   <  1:2864919/69‑1 (MQ=255)
                                        aGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGc                   <  1:367641/69‑1 (MQ=255)
                                        aGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGc                   <  1:1662129/69‑1 (MQ=255)
                                        aGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGc                   <  1:519741/69‑1 (MQ=255)
                                        aGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGc                   <  1:2399160/69‑1 (MQ=255)
                                        aGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGc                   <  1:2793/69‑1 (MQ=255)
                                         gACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCa                  >  1:2021788/1‑69 (MQ=255)
                                         gACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCa                  >  1:2302124/1‑69 (MQ=255)
                                           cccAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCAtt                <  1:3209255/69‑1 (MQ=255)
                                           cccAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCAtt                <  1:389637/69‑1 (MQ=255)
                                           cccAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCAtt                >  1:2906728/1‑69 (MQ=255)
                                           cccAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCAtt                <  1:2906998/69‑1 (MQ=255)
                                           cccAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCAtt                <  1:930353/69‑1 (MQ=255)
                                           cccAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCAt                 >  1:3247783/1‑68 (MQ=255)
                                              aGCGGATTGAAGGGTTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTg              >  1:2817775/1‑68 (MQ=255)
                                              aGCGGATGGAAGGGTTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGAc                     >  1:1809842/1‑61 (MQ=255)
                                              aGCGGATGGAAGGGTTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCAttt               >  1:2658399/1‑67 (MQ=255)
                                              aGCGGATGGAAGGGTTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTg              >  1:1432990/1‑68 (MQ=255)
                                              aGCGGATGGAAGGGTTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTg              >  1:1575246/1‑68 (MQ=255)
                                              aGCGGATGGAAGGGTTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTg              >  1:1075289/1‑68 (MQ=255)
                                              aGCGGATGGAAGGGTTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTg              >  1:1381065/1‑68 (MQ=255)
                                              aGCGGATGGAAGGGTTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTg              >  1:1084069/1‑68 (MQ=255)
                                              aGCGGATGGAAGGGTTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTg              >  1:2703112/1‑68 (MQ=255)
                                              aGCGGATGGAAGGGTTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTg              >  1:583497/1‑68 (MQ=255)
                                              aGCGGATGGAAGGGTTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTg              >  1:373527/1‑68 (MQ=255)
                                              aGCGGATGGAAGGGTTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTg              >  1:1942559/1‑68 (MQ=255)
                                              aGCGGATGGAAGGGTTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTg              >  1:3226410/1‑68 (MQ=255)
                                              aGCGGATGGAAGGGTTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTg              >  1:3186397/1‑68 (MQ=255)
                                              aGCGGATGGAAGGGTTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTg              >  1:2237159/1‑68 (MQ=255)
                                              aGCGGATGGAAGGGTTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTg              >  1:3031379/1‑68 (MQ=255)
                                              aGCGGATGGAAGGGTTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTg              >  1:3006595/1‑68 (MQ=255)
                                              aGCGGATGGAAGGGTTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTg              >  1:296373/1‑68 (MQ=255)
                                              aGCGGATGGAAGGGTTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTg              >  1:2741779/1‑68 (MQ=255)
                                                         ggACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCTCAATGAAAAg  <  1:2226701/69‑1 (MQ=255)
                                                           |                                                                  
GGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCTCAATGAAAAG  >  minE/247063‑247188

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: