Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,168,802 G→A E31K (GAA→AAA)  yegS → lipid kinase YegS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,168,8020GA92.9% 26.7 / ‑3.4 14E31K (GAA→AAA) yegSlipid kinase YegS
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (6/7);  total (7/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.05e-01

GGTAACCACTGAATCCTCCATACTTGAACTGACTTTTTATCCTCCGACTTTCATCCTGTTCTGACTCCACCTTTTGTTTTCTGCTCTACACTATCTACAGACCAATCATAAAGGCACATACGATCATGGCAGAATTTCCCGCCAGCTTACTGATTCTTAATGGCAAAAGTACTGACAATCTACCCTTGCGCGAAGCAATTATGCTGTTGCGTGAGGAAGGAATGACGATCCATGTGCGGGTCACCTGGGAGAAAGGCGATGCCGCACGATATGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAATTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGT  >  NC_000913/2168587‑2169019
                                                                                                                                                                                                                       |                                                                                                                                                                                                                         
ggtAACCACTGAATCCTCCATACTTGAACTGACTTTTTATCCTCCGACTTTCATCCTGTTCTGACTCCACCTTTTGTTTTCTGCTCTACACTATCTACAGACCAATCATAAAGGCACATACGATCATGGCAGAATTTCCCGCCAGCTTACTGATTCTTAATGGCAAAAGTACTGACAATCTACCCTTGCGCGAAGCAATTATGCTGTTGCGTGAGAAAGGAATg                                                                                                                                                                                                                   <  1:80781/224‑1 (MQ=255)
        cTGAATCCTCCATACTTGAACTGACTTTTTATACAACGACTTTCATCCTGTTCTGACTCCACCTTTTGTTTTATGCTCTAAACTATCTACAGACCAATCATAAAGGCACATACGATCATGGCAGAATTTCCCGCCAGCTTACTGATTCTTAATGGCAAAAGTACTGACAATCTACCCTTGCGCGAAGCAATTATGCTGTTGCGTGAGAAAGGAATGACGATCCATGTGCGGGTCACCTGGGAGAAAGGCg                                                                                                                                                                                 >  1:146917/1‑250 (MQ=255)
               ctccATACTTGAACTGACTTTTTATCCTCCGACTTTCATCCTGTTCTGACTCCACCTTTTGTTTTCTGCTCTACACTATCTACAGCCCAATCATAAAGGCACATACGATCATGGCAGAATTTCCCGCCAGCTTACTGATTCTTAATGGCAAAAGTACTGACAATCTACCCTTGCGCGAAGCAATTATGCTGTTGCGTGAGAAAGGAATGACGATCCATGTGCGGGTCACCTGGGAGAAAGGCGATGCCGc                                                                                                                                                                          >  1:155117/1‑250 (MQ=255)
                                                              gACTCCACCTTTTGTTTTCTGCTCTACACTATCTACAGACCAATCATAACGGCACATACGATCATGGCAGAATTTCCCTCCAGCTTACTGATTCTTAATGGCAAACGTACTGACACTCTACCCTTGCGCGAAACAATTATGCTGTTACGTGAGAAAGGAATGACGATCCATGTGCGGGTCACCTGGGAGACAGGCGATGCCGCACGATCTGTAGAGGAGGCACGGAAGTTTGGCGTCGCAACGGCGATTGc                                                                                                                          >  1:131271/1‑251 (MQ=255)
                                                                                                      cAATCATAAAGGCACATACGATCATGGCAGAATTTCCCGCCAGCTTACTGATTCTTAATGGCAAAAGTACTGACAATCTACCCTTGCGCGAAGCAATTATGCTGTTGCGTGAGAAAGGAATGACGATCCATGTGCTGGTCACCTGGGAGAAAGGCGATGCCGCACGATATGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCg                                                                                  <  1:119070/251‑1 (MQ=255)
                                                                                                                   ngcnnctatngtgtctgctgttCCCGCCAGCTAANGGATTCTTAATGGCAAAAGTACAGACAATCTNCNCGTGCGCGAAGCAATTATGCTGTGGCGTGAGAAAGGAATGACGATCCANGTGCGGGTCANCTGNNAGAANGNCNATGCCGCCCGATATGTAGAGGATGCCCGGNAGTTTGGCGTCGCAACGGTGAGTNNCGNNNGTNNCGATg                                                                                                            <  2:131271/192‑1 (MQ=255)
                                                                                                                   acaTACGATCATGGCAGAATTTCCCGCCAGCTTACTGATTCTTAATGGCAAAAGTACTGACAATCTACCCTTGCGCGAAGCAATTATGCTGTTGCGTGAGAAAGGAATGACGATCCATGTGCGGGTCACCTGGnag                                                                                                                                                                                        >  2:108160/1‑136 (MQ=255)
                                                                                                                              tGGCAGAATTTCCCGCCAACTTACTGATTCTGAATGGCAAAACGACTGACAATCTTCCCTTGCGCGAAGCATTGATTCTGTGGCGTGAGAAAGGAATGACGTGCCATGTTCTGTTCACCTGGGAGAAAGGCGATGCCGCACTATATGTAGAGGAGGCCCGGAAGttt                                                                                                                                              <  1:169620/167‑1 (MQ=255)
                                                                                                                                  aGAATTTCCCGCCAGCTTACTGATTCTTAATGGCAAAAGTACTGACAATCTACCCTTGCGCGAAGCAATTATGCTGTTGCGTGAGAAAGGAATGACGATCCATGTGCGGGTCACCTGGGAGAAAGGCGATGCCGCACGATATGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATAc                                                      <  1:1952/251‑1 (MQ=255)
                                                                                                                                                   tACTGATTCTTAATGGCAAAAGTACTCACAATCTACCCTTGCGCGAAGCAATTATGCCGTAGCGTGAGAAAGGAATGACGATCCATCTGCGGGTCACATGGGAGAACGGCGATGACGCACGATATGTAGAGGAGGCCCGGAAGTTTGGCCTCGCAACGGTGATTGCCGGTGGTGGCGATAGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAATTTTg                                     >  1:77540/1‑251 (MQ=255)
                                                                                                                                                    aCTGATTCTTAATGGCAAAAGTACTGGCAATCTAACCGTGCGCGAAGCAAGGATGCTGTTGCGGGAGAAAGGAATGACGACCCATGTGCGGGTCACCTGGGAGAAAGGCGATGCCGCACGATAGGTAGAGGAGGCGCGGAAGTTTGGTGTCGTATCGGTGATTGCCGGTGGTGGCTATGGCAGCATTAATGAAGTTTCTAGGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCTCTGGGAATTTTg                                     <  1:130824/250‑1 (MQ=255)
                                                                                                                                                             tAATGGCAAAAGTACTGACAAGCTACCCTTGCGCGAAGCAATTATGCTGTTGCGTGAGAAAGGAATGACGATCCATGTGCGGGTCACCTGGGAGAAAGGCGATGGCGCACGAGATGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATGGCCGGTGGTGGCGATGGCACCAGTAATGAAGTGTCTACGGCGTTGATTCAGTTTGATGGGGATGACATACCCGCGCTGGGAATTTTGCCATTAg                              <  1:65555/248‑1 (MQ=255)
                                                                                                                                                                   cAAAAGTACTGACAATCTACCCTTGCGCGAAGCAATTATGCTGTTGCGTGAGGAAGGAATGACGATCCATGTGCGGGTCACCTGGGAGAAAGGCGATGCCGCACGATATGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGACTTTTGCCATTAGGAACCGcc                      >  1:3122/1‑250 (MQ=255)
                                                                                                                                                                                       ccTTGCGCGAAGCAATTATGCTGTTGCGTGAGAAAGGAATGACGATCCATGTGCGGGTCACCTGGGAGAAAGGCGATGCCGCACGATATGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAATTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAgtgt  >  1:84750/1‑250 (MQ=255)
                                                                                                                                                                                                                       |                                                                                                                                                                                                                         
GGTAACCACTGAATCCTCCATACTTGAACTGACTTTTTATCCTCCGACTTTCATCCTGTTCTGACTCCACCTTTTGTTTTCTGCTCTACACTATCTACAGACCAATCATAAAGGCACATACGATCATGGCAGAATTTCCCGCCAGCTTACTGATTCTTAATGGCAAAAGTACTGACAATCTACCCTTGCGCGAAGCAATTATGCTGTTGCGTGAGGAAGGAATGACGATCCATGTGCGGGTCACCTGGGAGAAAGGCGATGCCGCACGATATGTAGAGGAGGCCCGGAAGTTTGGCGTCGCAACGGTGATTGCCGGTGGTGGCGATGGCACCATTAATGAAGTTTCTACGGCGTTGATTCAGTGTGAGGGGGATGACATACCCGCGCTGGGAATTTTGCCATTAGGAACCGCCAATGATTTTGCCACCAGTGT  >  NC_000913/2168587‑2169019

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