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,9650CT90.0% 47.5 / ‑6.0 20A50T (GCC→ACC) ppsAphosphoenolpyruvate synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (8/10);  minor base G (2/0);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCACCGCAAAAGAGGCGTTTTCGTCATCGGCGGAAAGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGA  >  NC_000913/1786731‑1787203
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                               
agCGCACCGCAAAAGAGGCGTTTTCGTCATCGGCGGAAAGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGGGGTGAAACCATTcgggacg                                                                                                                                                                                                                                 >  1:180080/1‑251 (MQ=255)
                                    aaaGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGTGAAAACCCTCGGAACGGAAACAACCCATCCGGGAAGATTAGTAATCATTTcc                                                                                                                                                                                            >  2:90736/1‑250 (MQ=255)
                                          gtgCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGTGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAgg                                                                                                                                                                                      <  2:248994/251‑1 (MQ=255)
                                                              gCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGTGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGGAATCATTTCACCCCGGGAGGCATTTTTTCCCCCAAc                                                                                                                                                                  >  1:187280/1‑251 (MQ=255)
                                                                 ttttCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTTGTCCAGAAACTGGTTAAACGCGTCGGCGGTTTTGGGGAAAACAATTGGCACGGAAACACCCCATCCGGAAAGATTAGTAATCATTTTACCCAAGGAGGGAtttttttcccccacccc                                                                                                                                                               >  2:300425/1‑239 (MQ=255)
                                                                      cAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGTGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTa                                                                                                                                                          <  1:90736/251‑1 (MQ=255)
                                                                       aGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGTGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTAc                                                                                                                                                         <  2:180080/251‑1 (MQ=255)
                                                                           cAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAAAATCGACAATATCCGTTTTATCCAACAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGTGAAACCATTCGGAAAGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCCCCCAGGGAGGCCTTTTTGCCCCCAACCCTGTCCAcatcc                                                                                                                                                     >  1:188396/1‑250 (MQ=255)
                                                                                        aGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGTGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGtt                                                                                                                                        <  1:105360/251‑1 (MQ=255)
                                                                                             cGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGTGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTta                                                                                                                                   <  1:156952/251‑1 (MQ=255)
                                                                                                                     cgcCCGCGTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTATGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGTGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGatc                                                                                                          <  2:227593/251‑3 (MQ=255)
                                                                                                                                          tAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGTGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACa                                                                                       >  1:57286/1‑250 (MQ=255)
                                                                                                                                             cATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGTGAAAACATTCGGAACGGAAACACCCCTTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATcc                                                                                   >  2:269474/1‑251 (MQ=255)
                                                                                                                                                                      cagTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGTGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTgag                                                          <  2:187280/251‑1 (MQ=255)
                                                                                                                                                                                  gcgcTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGGGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGGAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATAACAAAGCACCAGCGGTGAAGAGACATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGGTTTAGAAACACCTTTCt                                              >  2:294824/1‑251 (MQ=255)
                                                                                                                                                                                                            gAAACTGGTTAAACGCGTCGGCGGTTGTGGTGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAAcc                    <  1:260491/251‑1 (MQ=255)
                                                                                                                                                                                                              aaCTGGTTAAACGCGTCGGCGGTTGTGGTTAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCgg                  <  1:306822/251‑1 (MQ=255)
                                                                                                                                                                                                                         cgcgTCGGCGGTTGTGGTGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAAc       <  1:67748/251‑1 (MQ=255)
                                                                                                                                                                                                                               ggcggTTGTGGTGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGCGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGa  >  1:1758/1‑250 (MQ=255)
                                                                                                                                                                                                                               ggcggATGTGGTGAAACCATTCGGAACGGAAACACCCCTTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGa  >  2:180581/1‑250 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                               
AGCGCACCGCAAAAGAGGCGTTTTCGTCATCGGCGGAAAGCTGTGCATAGGCTTCGCGGATGGCGTTTTCCAGCTCAGGCTGGAAGGGAGTGTCGATAATCCACTGGCGGATTTGCGCGCCCGCTTTCGCAAGCTGAGTAACATCGTCAATATCCGTTTTATCCAGCAGTTCATAAATGCGCTGGTTTACGCCGCTTTGGTCCAGAAACTGGTTAAACGCGTCGGCGGTTGTGGCGAAACCATTCGGAACGGAAACACCCATTCCGGAAAGATTAGTAATCATTTCACCCAGGGAGGCATTTTTGCCCCCAACCCTGTCTACATCATTCATGCCGAGTTGGTTATACCAAAGCACCAGCGGTGACGAGCCATTGTTGGACATCGAACAATCCTTTTGTGATAAATGAACGGTTTGAGAAACACATTTCTGCGCATTTATCTTTGCATATTTAACCGGATGAAAAAAACGGTGA  >  NC_000913/1786731‑1787203

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