Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,786,965 C→T A50T (GCC→ACC)  ppsA ← phosphoenolpyruvate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,786,9650CT85.2% 74.9 / ‑5.0 27A50T (GCC→ACC) ppsAphosphoenolpyruvate synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (2/21);  minor base G (4/0);  total (6/21)
Fisher's exact test for biased strand distribution p-value = 8.55e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCGCACCGCAAAAGAGGCGTTTTCGTCATCGGCGGAAAGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTG‑TGGCGAAACCATTCGGAACGGAAACACCCATTCCGG‑AAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGAT  >  NC_000913/1786730‑1787189
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                 
gagCGCACCGCAAAAGAGGCGTTTTCGTCATCGGCGGAAAGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTG‑TGGTGAAACCATTcggaac                                                                                                                                                                                                                    <  1:224439/251‑1 (MQ=255)
          acaatgagGCGTTTCTGTCCTCGGCGGTAAGCTGTGCATAGGTTTCGCGTATGGCGTTTTCCTGCTCAGGCTGGAAGGGAGTTTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTGCGCAAGCTGATTAACATCGTCAATATCCTTTGTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCTCTTTGGTCCAGAAACTGGTTAAACGCGTTGGCGGTTG‑TGGTGAAACCATTcggttctgtagcaccc                                                                                                                                                                                                          <  2:381535/246‑14 (MQ=255)
                                 gcggAAAGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCAGCGGTTG‑TGGTGAAACCATTCGGAACGGAAACACCCATTCCGG‑AAAGATTAGTAATCAt                                                                                                                                                                                  <  1:87457/251‑1 (MQ=255)
                                     aaagctggtgctaagCTTCGCGGATGGCGTTTACCAGCACAGGCTGGAAGGGAGTGACGCTAATCCACTGGCGGATGTGCGCTCCCGCTTTCGCAAGCTTAGTAACATCGTCAATATCCGTTTTAGTCAGCAGTTCATAAATGCGCTGGTTTACTCCGCGTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTG‑TGGTGAAACCATTCGGAACGGAAACCCCCATTCCGG‑AAAGATTAGTAATCACTTCa                                                                                                                                                                              <  2:256811/237‑1 (MQ=255)
                                                   gCTTCGCGGATGGCGGTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTG‑TGGTGAAACCATTCGGAACGGAAACACCCATTCCGG‑AAAGATTAGTAATCATTTCACCCAGGGAGGCAtt                                                                                                                                                                <  1:440785/251‑1 (MQ=255)
                                                       cgcgGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTG‑TGGTGAAACCATTCGGAACGGAAACACCCATTCCGG‑AAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTg                                                                                                                                                            <  1:496104/251‑1 (MQ=255)
                                                                            cAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCCGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCCGAAACTGGTTAAACGCGTCGGCGGTTGGTGGGTAAACC‑TTCGGAACGGGAAAACCCCTTACGG‑GGAGACTAATTATCATTTCTCCCCTGGGGGGATTTTTTTCCCCAACCCTGTCTaacac                                                                                                                                        >  2:17730/1‑246 (MQ=255)
                                                                                  ggAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTG‑TGGTGAAACCATTCGGAACGGAAACACCCATTCCGG‑AAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACAtcattca                                                                                                                                   <  1:92782/249‑1 (MQ=255)
                                                                                           tgtCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCCGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAACCTGGTTAAAAGCGTCGGCGGTTG‑TGGGGAATACATTCGGGAAGGAGAAACCCCTTCCGGGACAG‑TTTGTAATCATTTCACCCCCGGAGGGACTTTTGCCCCCCACCCCGTCTCCAacattcttcccccggtgg                                                                                                                        >  2:100333/1‑239 (MQ=255)
                                                                                                            cGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACAGCGGCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGGTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTG‑TGGTGAAACCATTCGGAACGGAAACACCCATGCCGG‑AAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAg                                                                                                       <  1:218823/251‑1 (MQ=255)
                                                                                                               aTTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTG‑TGGTGAAACCATTCGGAACGGAAACACCCATTCCGG‑AAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGt                                                                                                                       <  1:195852/232‑1 (MQ=255)
                                                                                                               aTTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTG‑TGGGGAAACCATTCGGAACGGAAACACCCATTCCGG‑GAAGATTAGTAATCATTTCACCCCGGGAGGCATTTTTGCCCCCCACCCTGTCTACCTCCTTCCTTCCGAggtggtagatcggacgggcgccggg                                                                                                    >  2:195852/1‑232 (MQ=255)
                                                                                                                          cGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTG‑TGGTGAAACCATTCGGAACGGAAACACCCATTCCGG‑AAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCAt                                                                                         <  1:30335/251‑1 (MQ=255)
                                                                                                                              ttCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTG‑TGGTGAAACCATTCGGAACGGAAACACCCATTCCGG‑AAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCAttgt                                                                                      <  1:100333/250‑1 (MQ=255)
                                                                                                                                 gCAAGCTTAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTG‑TGGTGAAACCATTCGGAACGGAAACACCCATTCCGG‑AAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGa                                                                                  <  1:506274/251‑1 (MQ=255)
                                                                                                                                                  gTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTTTTTTACGCTGTTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTG‑TGGTGAAACCATTCGGAACGGAACCACCCATTCCGG‑AAAGATTAGTAATCATTGCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCtttt                                                                 <  1:297294/251‑1 (MQ=255)
                                                                                                                                                                                                ccgcTTGGGTCCAGAAACTGATTAAACGCGGCGGCGGTTG‑TGGTGAAACCAGTCGGAACGGAAACACCCATTCCGG‑AAAGATTAGTAATCATTTCCCCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCtt                   <  2:160114/251‑1 (MQ=255)
                                                                                                                                                                                                    tttGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTG‑TGGTGAAACCATTCGGAACGGAAACACCCATTCCGG‑AAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGa                                                                                                 >  1:281886/1‑169 (MQ=255)
                                                                                                                                                                                                    tttGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTG‑TGGTGAAACCATTCGGAACGGAAACACCCATTCCGG‑AAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGa                                                                                                 <  2:281886/169‑1 (MQ=255)
                                                                                                                                                                                                     ttGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTG‑TGGTGAAACCATTCGGAACGGAAACACCCATTCCGG‑AAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCa               <  1:17730/250‑1 (MQ=255)
                                                                                                                                                                                                     ttGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTG‑TGGTGAAACCATTCGGAACGGAAACACCCATTCCGG‑AAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCGTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGc                <  1:324241/249‑1 (MQ=255)
                                                                                                                                                                                                              aaaCTGGTTAAACGCGTCGGCGGTTG‑TGGTGAAACCATTCGGAACGGAAACACCCATTCCGG‑AAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCg     <  1:94729/251‑1 (MQ=255)
                                                                                                                                                                                                                aCTGGTTAAACGCGTCGGCGGTTG‑TGGGGGAACCATTAGGAAAGGAAACACCCACTCCGG‑AAAGATTAAGAATCATTTCCCCCAGGGAGGCCTTTTTTCCCCCAACCCTGTCTACATCATTCATGCCGAGATGGGTATACCCAAACACCAGCGGGGACGAGCCATTGTTGGACATAGAACAATCCTTTTGTGATAAATGAACCGTTTGTGAAAAACATTTCTTCGCCTTTATCTTTTCATATTTAAcccgt   >  1:247003/1‑248 (MQ=255)
                                                                                                                                                                                                                  gggTTAAACGCGTCGGCGGTTG‑TGGTGAAACCATTCGGAACGGAAACACCCATTCCGG‑AAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCgg                                                                                                    <  1:352632/151‑1 (MQ=255)
                                                                                                                                                                                                                   ggTTAAACGCGTCGGCGGTTG‑TGGTGAAACCATTCGGAACGGAAACACCCATTCCGG‑AAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGAt  <  1:314742/249‑1 (MQ=255)
                                                                                                                                                                                                                                   ggTTG‑TGGTGAAACCATTCGGAACGGAAACACCCATTCCGG‑AAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATTCCAAAGCACCAGCGGTGACGAGcc                                                                                           <  2:314867/144‑1 (MQ=255)
                                                                                                                                                                                                                                   ggTTG‑TGGTGAAACCATTCGGAACGGAAACACCCATTCCGG‑AAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGcc                                                                                           >  1:314867/1‑144 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                 
GAGCGCACCGCAAAAGAGGCGTTTTCGTCATCGGCGGAAAGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTG‑TGGCGAAACCATTCGGAACGGAAACACCCATTCCGG‑AAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGAT  >  NC_000913/1786730‑1787189

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