Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,805,830 C→T E41K (GAA→AAA)  waaP ← lipopolysaccharide core heptose (I) kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,805,8300CT93.3% 35.2 / NA 15E41K (GAA→AAA) waaPlipopolysaccharide core heptose (I) kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (7/7);  minor base G (0/1);  total (7/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GGTAATAATAAATGATGCGCGAGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGTTTTTTGACCTCTGCAAAAGCATCTTTACCACGCCAAAGTGTGGCAAGCGGCTCTTTTAATTCAACCATCTAAACCACCTGTAATGATATCCGCGGCTTTTTCTGGCAAGC  >  NC_000913/3805612‑3805992
                                                                                                                                                                                                                          |                                                                                                                                                                  
ggTAAGAATAAATGATGCGCGAGTGAATTGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCATGCGATGAATAGCGTTCCTCTCGCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATGGTGTATAATCTCTTGTAATGTCGTCCCCTTGCGCCATTTAAGAAAAAAGCGTTTCCCGGACAGTGTAAAGCGTAATGTACGACGAg                                                                                                                                                <  1:429740/239‑1 (MQ=255)
                   ngAGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTTAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTc                                                                                                                   >  1:100259/2‑249 (MQ=255)
                                                                     aTCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTTAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGTTTTTTGACCTCTGCAAAAGCATCTTTACCACGCCAAAGTGTGGCAAGc                                                                <  1:79355/250‑1 (MQ=255)
                                                                                       cAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTTAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGTTTTTTGACCTCTGCAAAAGCATCTTTACCACGCCAAAGTGTGGCAAGCGGCTCTTTTAATTCAAc                                               <  1:432886/249‑1 (MQ=255)
                                                                                                         ccACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGCGTAGATTTTTTCTAATCTCTTTTCATGTGGTACCCTTGTGCCATTTAAGAAAATAGCTTTTCCAGGACAGTTTAAAGCGCCATGTACGACGAGTCTCCAGTTCACGAACGACCTCGCCGTTCAGTTTTTTGACCTCTGCAACAGCATCTTTACCACGACAACGTGTGGCAAGCGGCTCTTTTAATTAAACCATCTAAACCCCCTGTAATg                           >  1:248733/1‑251 (MQ=255)
                                                                                                         ccACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTGCCCGGACAGTTTAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACGTCGCCGTTCAGTTTTTGGACCTCTGCAAAAGCATCTTTACCGCGCCAAAGTGTGGCAAGCGGCTCTTTTAATTCAACCATCTAAACCACCTGTAATg                           <  1:281208/251‑1 (MQ=255)
                                                                                                          cACTCTCTGTCTGCGCCCAAAACGGGCTTCCGCAATGAGAGTATATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCTGACAGTTTAAAGCGTAATGTACGACTAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGTTTTTTGACCTCTGCAAAAGCa                                                                                           <  1:374204/186‑1 (MQ=255)
                                                                                                          cACTCNNTGNNNGCNNCCAAAACGGGCATCCGCAATGAGNGTAGATTTTTTATAATCTCTTTTAATGTCNTNCNCTTGNNCCANTTAAGAAAATNGCTTTTCCCGGACAGTTTAAAGCGTAATGTACGACGAGTCTCCAGTTCNCNAAAGACCTCGCCGTTCAGTTTTTTGACCTCTNCAAAAGCa                                                                                           >  2:374204/1‑186 (MQ=255)
                                                                                                                 tgtcttCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTGTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTGAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGTCGTTCAGTTTTTTGACCTCTGCAAAAGCATCTTTAcn                                                                                   <  1:90035/187‑2 (MQ=255)
                                                                                                                      gcgcCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCATTGTGCCCTTTAAGAAAATAGATTTTCCCGGACAGTTTAAAGCGTAATGTACGACGAGTCTCCAGTTCCCGAAAGACCTCGACGTTCAGTTTTTTGACCTCTGCAAAAGCATCTTTACCAAGCCAAAGTGTGGCAAGCGGCTCTTTTATTTCAACCATCTAAACCACCTGTAATGATATCCGCGGCtt              >  1:344331/1‑251 (MQ=255)
                                                                                                                                   gCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTTAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGTTTTTTGACCTCTGCAAAAGCATCTTTACCACGCCACAGTGTGGCAAGCGGCTCTTTTAAtt                                                   >  1:447455/1‑201 (MQ=255)
                                                                                                                                    cATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTTAAAGCg                                                                                                                                                              >  1:418411/1‑93 (MQ=255)
                                                                                                                                                                            gtcgtcCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTTAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGTTTTTTGACCTCTGCAAAAGCATCTTTACCACGCCAAAGTGTGGCAAGCGGCTCTTTTAATTCAACCATCTAAACCACCTGTAATGATATCCGCGGCTTTTTCTGGCAAGc  <  1:82800/209‑1 (MQ=255)
                                                                                                                                                                                   ccTTGTGCCATTTACGAAAATAGGTTTTCCCGGACAGTTTAAAGCGTAATGTACGACGAGTCTCCAGTTCACGACAGACCTCGCCGTTCAGTTTTTTGACCTCTGCAAAAGCATCTTTACCACGCCaa                                                                            >  1:178438/1‑128 (MQ=255)
                                                                                                                                                                                                                     cAGTTTAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGTTTTTTGACCTCTGCAAAAGCATCTTTACCACGCCAAAGTGTGGCAAGCGGCTCTTTTAATTCAACCATCTAAACCACCTGTAATGATATCCGCGGCTTTTTCTgg       <  1:364796/163‑1 (MQ=255)
                                                                                                                                                                                                                     cAGTTNNAANNNTANNGTACGACGAGTCTCCAGTTCACGNAAGACCTCGCCGTTCAGTTTTTTGACCTCNGNANAAGCNNCTTNACCACGCCAANGTGTGGCAAGCGGCTCTTTTAATTCAACCATCTAAACCACCTGTAATGNTNTCCGCGGCTTTTTCTgg       >  2:364796/1‑163 (MQ=255)
                                                                                                                                                                                                                          |                                                                                                                                                                  
GGTAATAATAAATGATGCGCGAGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGTTTTTTGACCTCTGCAAAAGCATCTTTACCACGCCAAAGTGTGGCAAGCGGCTCTTTTAATTCAACCATCTAAACCACCTGTAATGATATCCGCGGCTTTTTCTGGCAAGC  >  NC_000913/3805612‑3805992

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