Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,141,133 C→T 28.7% E363E (GAG→GAA)  ptsA ← putative PTS multiphosphoryl transfer protein PtsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,141,1330CT28.7% 12.4 / 6.0 14E363E (GAG→GAA) ptsAputative PTS multiphosphoryl transfer protein PtsA
Reads supporting (aligned to +/- strand):  ref base C (4/6);  new base T (2/2);  total (6/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 = 4.01e-01

TTGGCGGATTCCAGCGCCTGACAAAAAATGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGT  >  NC_000913/4140912‑4141362
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                     
ttGGCGGAGTCCAGCTCCTAACAAAAAATGTTGTACAACTCGCTTTCGGCCGGTGCGCTGGTGCTATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCAACATTGCCGAATGCGGCCTGCGCTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGc                                                                                                                                                                                  <  1:123785/275‑1 (MQ=255)
                                                               cGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGGCGTGCGCGATTGTAGCGGCAATTTCACTGCGGATACCGTCCGCGGTACGGGCTTGTTGGGTCAGCCAGACACGGTGTTGCTCCCGCAg                                                                                                                                                                                                                               >  1:221223/1‑167 (MQ=255)
                                                                        taGAGCATTTCAGTGCTGAACAAACCAACGCCTTCCGCACCATTGACGAATGCGGCCTGCGCTTGCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCGTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGTGTACTGTACGCGTGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAg                                                                                                          <  1:6363/275‑1 (MQ=255)
                                                                                ttCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACGGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCCGGGCGTCCTGTGCGCGGGCTTATTGCTGata                                                                                                                                                                                                 >  1:17997/1‑180 (MQ=255)
                                                                                    gggCGGAACAAACCACCGCCTTCCGCACCATTGCCGAATGCGGCCGGCGCTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGAtaataa                                                                                                                                                                                             <  1:154621/178‑1 (MQ=255)
                                                                                                    cGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGTTCACGCCGGGCGTCCTGTACGCGGGcttc                                                                                                                                                                                                          >  1:23930/1‑151 (MQ=255)
                                                                                                    cGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGTTCACGCAGGGCGTCCTGTACGCGGGcttc                                                                                                                                                                                                          >  1:23915/1‑151 (MQ=255)
                                                                                                        ttCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGTTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGc                                                                                                                                                                                                      <  2:23915/151‑1 (MQ=255)
                                                                                                        ttCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGTTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGc                                                                                                                                                                                                      <  2:23930/151‑1 (MQ=255)
                                                                                                                                                  aGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGCGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCCACCACAATCGCCCCGGCGTTGCCGTCGATATAAATTGTTTGTTGCTGCCACGGGGTAAGGGCCTCAATATCCacac                                                                                         >  1:9585/1‑218 (MQ=255)
                                                                                                                                                             cGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACg                                                                                                                                                                                                                                  >  1:202027/1‑70 (MQ=255)
                                                                                                                                                                 aaTTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAgg                                                                                                                                                                                                                              <  2:202027/70‑1 (MQ=255)
                                                                                                                                                                                                              ccAGACACGCTGTTGCTCACGCAGGGCGTCCTTTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCg                                                                                <  1:165230/167‑1 (MQ=255)
                                                                                                                                                                                                                        tGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGt  <  1:185786/235‑1 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                     
TTGGCGGATTCCAGCGCCTGACAAAAAATGTTGTACAACTCGCTTTCGCCCGGTGCGCTGGTGCGATCCATATAGAGCATTTCAGTGCGGAACAAACCAACGCCTTCCGCACCATTGCCGAATGCGGCCTGCGCTTCCACGGAGTGAGCGATGTTAGCGGCAATTTCAATGCGGATACCGTCAGCGGTACGGGCTTGTTGGGTCAGCCAGACACGCTGTTGCTCACGCAGGGCGTCCTGTACGCGGGCTTCTTGCTGATAATAACGAGCTACGGCTTCCCCTGGCTCAACCACAATCGCCCCGGCGTTGCCGTCGATATAAATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGT  >  NC_000913/4140912‑4141362

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑