Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,171,411 G→A E310K (GAG→AAG)  qseC → sensor histidine kinase QseC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,171,4110GA94.4% 41.6 / ‑5.9 18E310K (GAG→AAG) qseCsensor histidine kinase QseC
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (12/5);  total (12/6)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

GAACGACGCTTTACCTCCGACGCAGCTCACGAACTTCGTAGCCCGTTAACGGCGCTGAAAGTGCAAACCGAAGTTGCGCAGCTCTCTGACGATGATCCGCAGGCGCGGAAAAAAGCACTGCTCCAATTACATTCCGGGATCGATCGCGCTACTCGTCTGGTTGATCAACTGCTCACGCTATCGCGGCTGGACTCACTGGATAACCTTCAGGACGTCGCGGAGATCCCGCTTGAAGATCTCCTGCAATCGTCGGTGATGGATATTTACCACACGGCGCAGCAGGCGAAAATTGACGTGCGACTGACACTCAATGCCCACAGCATCAAACGCACCGGGCAACCGCTATTGCTAAGTTTGTTGGTGCGAAATTTGCTGGATAACGCCGTGCGCTACAGTCCACAGGGCAGCGTGGTAGACGTCACGCTGAATGCT  >  NC_000913/3171192‑3171623
                                                                                                                                                                                                                           |                                                                                                                                                                                                                    
ggaacacctcttACCTCCGACGCAGCTCACGAACTTCGTAGCCCGTTAACGGCGCTGAAAGTGCAAACCGAAGTTGCGCAGCTCTCTGACGATGATCCGCAGGCGCGGAAAAAAGCACTTCTCCAATTACATTCCGGGATCGATCGCGCTACTCGTCTGTTTGATCAACGGCTCACGCTATCGCCGCTTGACTCACTGGATAACCTTCAGGACGTTGCGGAGATCCCGCTTGAAGATCTCCTGCAAtcgtc                                                                                                                                                                                       <  1:223698/241‑1 (MQ=255)
                                aCTTCGTAGCCCGTTAACGGCGCTGAAAGTGCAAACCGAAGTTGCGCAACTCTCTGACGATGATCCGCAGGCGCGGAAAAAACCACTGCTCCACTTACATTCCGGGCTCGATCGCGCTACTCGTCTGGATGTTCCACTGCTCCCGCTATCGCGGCTGGACTCACTGGATAACCTTCAGGACGTCGCGAAGCTCCCGCTTGAAGATCTCCTGCAATCGTCGGTGATGGATATTTACCACACGGCGCAGCAgg                                                                                                                                                       >  1:121847/1‑251 (MQ=255)
                                   tCGTAGCCCGTTAACGGCGCTGAAAGTGCAAACCGAAGTTGCGCAGCTCTCTGACGATGATCCGCAGGCGCGGAAAAAAGCACTGCTCCAATTACATTCCGGGATCGATCGCGCTACTCGTCTGGTTGATCAACTGCTCACGCTATCGCGGCTGGACTCACTGGATAACCTTCAGGACGTCGCGAAGATCCCGCTTGAAGATCTCCTGCAATCGTCGGTGATGGATATTTACCACACGGCGCAGCAGGc                                                                                                                                                      <  1:205463/249‑1 (MQ=255)
                                           cGTTAACGGCGCTGAAAGTGCAAACCGAAGTTGCGCAGCTCTCTGACGATGATCCGCAGGCGCGGAAAAAAGCACTGCTCCAATTACATTCCGGGATCGATCGCGCTACTCGTCTGGTTGATCAACTGCTCACGCTATCGCGGCTGGACTCACTGGATAACCTTCAGGACGTCGCGAAGATCCCGCTTGAAGATCTCCTGCAATCGTCGGTGATGGATATTTACCACACGGCGCAGCAGGCGAAAATTGAc                                                                                                                                            >  1:86896/1‑251 (MQ=255)
                                                   ncgcTGAAAGTGCAAACCGAAGTTGCGCAGCTCTCTGACGATGATCCGCAGGCGCGGAAAAAAGCACTGCTCCAATTACATTCCGGGATCGATCGCGCTACTCGTCTGGTTGCTCAACTGCTCACGCTATCGCGGCTGGACTCACTGGATAACCTTCAGGACGTCGCGAAGATCCCGCTTGAAGATCTCCTGCAATCGTCGGTGATGGATATTTACCACACGGCGCAGCAGGCGAAAATTGACGTGCGACt                                                                                                                                    >  1:28309/2‑251 (MQ=255)
                                                     gcTGAAAGTGCAAACCGAAGTTGCGCAGCTCTCTGACGATGCTCCGCAGGCGCGGAAAAAAGCACTGCTCCAATTACATTCCGGGATCGATCGCGCTACTCGTCTGGTTGATCAACTGCTCACGCTATCGCGGCTGGACTCACTGGATAACCTTCAGGACGTCGCGAAGATCCCGCTTGAAGATCTCCTGCAATCGTCGGTGATGGATATTTACCACACGGCGCAGCAGGCGAAAATTGACGTGCGACTg                                                                                                                                   >  1:86349/1‑250 (MQ=255)
                                                               cAAACCGAAGTTGCGCAGCTCTCTGACGATGCTCCGCAGGCACGGAAAAAAGCACTGCTCCAATTACATTCAGGGATCGATCGCGCTCCTCGTCTGGTTGATCAACTGCTCACGCTATCGCGGCTGGACTCACTGGATAACCTTCAGGACGTCGCGAAGATCCCGCTTGAAGATCTCCTGCCATCGTCGGTGATGGATATTTACCACACGgcgc                                                                                                                                                             >  1:102636/1‑214 (MQ=255)
                                                                 aaCCGCAGTTGCGCAGCTCTCTGACGATGATCCGCAGGCGCGGAAAAACGCACTGCTCCAATTACATTCCGGGATCGCTCGCGCTACTCGTCTGGTTGATCAACTGCTCACGCTATCGCGTCTGGACTCACTGGATACCCTTCAGGACGTCGCGAAGATCCCGCTTGAAGATCTCCTGCAATCGTCGGTGATGGATATTTACCACACGGCGCAGCAGGCGAAACTTGACGTGCGACTGACACTCAATGccc                                                                                                                      >  1:137200/1‑251 (MQ=255)
                                                                           gcgcAGCTCTCTGACGCTGATCCGCAGGCGCGGAAAAAAGCACTGCTCCAATTACATTCCCGGATCGATCGCGCTACTCGTCTGGTTGATCAACTGCTCACGCTATCGCGTCTGGACTCTCTGTATAACCCTCAGGACGTCGCGAATATTCCGCTTGAATCTCTCCTGCCATCGTCGGTTATGGATATTTACCACTCGCCGCCGCAGGCGAAAACTGATGTACGACTGACTCTCAATGCCCCCAGCATCa                                                                                                             <  1:5255/250‑1 (MQ=255)
                                                                            agcAGCTCTCTGACGATGATCCGCAGGCGCGGAACAAAGCACTGCTCCAATTACATTCCGGGATCGATCGCGCTACTCGTCTGGTTGATCAACTGCTCACGCTATCGCGGCTGGACTCACTGGATAACCTTCAGGACGTCGCGAAGATCCCGCTTGAAGATCTCCTGCAATCGTCGGTGATGGATATTTACCACACGGCGCAGCAGGCGAAAATTGACGAGCGACTGACACTCAATGCCCACAGCATCaa                                                                                                            >  1:71487/2‑250 (MQ=255)
                                                                             ncAGCTCTCTGACGATGATCCGCAGGCGCGGAAAAAAGCACTGCTCCAATTACATTCCGGGATCGATCGCGCTACTCGTCTGGTTGATCAACTGCTCACGCTATCGTGGCTGGACTCACTGGATAACCTTCAGGACGTCGCGAAGATCCCGCTTGAAGATCTCCTGCAATCGTCGGTGATGGATATTTACCACACGGCGCAGCAGGCGAAAATTGACGTGCGACTGACACTCAATGCCCACAGCATCAAAc                                                                                                          >  1:48237/2‑251 (MQ=255)
                                                                                                      gcgcGGAAAAAAGCACTGCTCCAATTACATTCCGGGATCGATCTCGCTACTCGTCTGGTTGATCAACTGCTCACGCTAGCGCGGCTGGACTCACTGGATAACCTTCAGGACGTCGCGAAGATCCCGCTTGAAGATCTCCTGCAATCTTCGGTGATGGAGATTTACCACACGGCTCAGCAGGCTAAAATTGTCGTGCGACTGACACTCAATGCCCACAGCATCAAACGCACCGGGCAACCGCTATTGCTAAg                                                                                 <  1:132884/251‑1 (MQ=255)
                                                                                                                                              atcgCGCTACTCGTCTGGTTGATAAACTGCTCACGCTATCGCGGCTGGACTCACTGGATAACCTTCAGGACGCCGCGAAGATCCCGCTTGACGATCTACTGCAATCGTCGGTGATGGATATATACCACACGGCGCAGCAGGCGAAAATTGACGTGCGACTGACACTCAATGCCCACAGCATCAAACGCACCGGGCAACCGCTATTGCTAAGTTTGTTGGTGCGAAATTTGCTGGATAACGCCGTgcgc                                            >  1:79244/1‑248 (MQ=255)
                                                                                                                                                         cGTCTGGTTGATCAACTGCTCACGCTATCGCGGCTGGACTCACTGGATAACCTTCAGGACGTCGCGAAGATCCCGCTTGAAGATCTCCTGCAATCGTCGGTGATGGATATTTACCACACGGCGCAGCAGGCGAAAATTGACGTGCGACTGACACTCAATGCCCACAGCATCAAATGCACCGGGCAACCGCTATTGCTAAGTTTGTTGGTGCGAAATTTGCTGGATAACGCCGTGCGCTACAGTCCACAggn                              <  1:47305/251‑2 (MQ=255)
                                                                                                                                                              ggTTGATCAACTGCTCACGCTATCGCGGCTGGACTCACTGGATAACCTTCAGGACGTCGCGAAGATCCCGCTTGAAGATCTCCTGCAATCGTCGGTGATGGATATTTACCACACGGCGCAGCAGGCGAAAATTNACGNNNNNCtnnnnnn                                                                                                                              >  2:159324/1‑144 (MQ=255)
                                                                                                                                                              ggTTGATCAACTGCTCACGCTATCGCGGCTGGACTCACTGGATAACCTTCAGGACGTCGCGAAGATCCCGCTTGAAGATCTCCTGCAATCGTCGGTGATGGATATTTACCACACGGCGCAGCAGGCGAAAATTGACGTGCGACTGACACTCAATGCCCACAGCATCAAACGCACCGGGCAACCGCTATTGCTAAGTTTGTTGGTGCGAAAtttn                                                              <  1:159324/214‑2 (MQ=255)
                                                                                                                                                                                                 cACTGGATAACCTTCAGGACGTCGCGAAGATCCCGCTTGAAGATCTCCTGCAATCGTCGGCGCTGGATATTTACCCCAAGGCGAAGCAGGCGCAAATTGACGTGCGACTGACACTCAATGCCCACAGCATCACACGCACCGGGCAACCGCTATTGCTAAGTTTGTTGGTGCGAAATTTGCTGGATAACGCCGTGCGCTACAGTCCACAGGGCAGCGTGGTAGACGTCACGCTGAATGCt  >  1:178910/1‑239 (MQ=255)
                                                                                                                                                                                                           ccTTCAGGACGTCGCGAAGATCCCGCTTGAAGATCTCCTGCAATCGTCGGTGATGGATATTTACCACACGGCGCAGCAGGc                                                                                                                                                      >  1:41149/1‑81 (MQ=255)
                                                                                                                                                                                                                           |                                                                                                                                                                                                                    
GAACGACGCTTTACCTCCGACGCAGCTCACGAACTTCGTAGCCCGTTAACGGCGCTGAAAGTGCAAACCGAAGTTGCGCAGCTCTCTGACGATGATCCGCAGGCGCGGAAAAAAGCACTGCTCCAATTACATTCCGGGATCGATCGCGCTACTCGTCTGGTTGATCAACTGCTCACGCTATCGCGGCTGGACTCACTGGATAACCTTCAGGACGTCGCGGAGATCCCGCTTGAAGATCTCCTGCAATCGTCGGTGATGGATATTTACCACACGGCGCAGCAGGCGAAAATTGACGTGCGACTGACACTCAATGCCCACAGCATCAAACGCACCGGGCAACCGCTATTGCTAAGTTTGTTGGTGCGAAATTTGCTGGATAACGCCGTGCGCTACAGTCCACAGGGCAGCGTGGTAGACGTCACGCTGAATGCT  >  NC_000913/3171192‑3171623

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: