Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,968,268 (C)5→6 coding (167/2247 nt) ptsP ← phosphoenolpyruvate‑protein phosphotransferase PtsP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,968,2631.C97.7% 148.9 / ‑3.2 44coding (172/2247 nt)ptsPphosphoenolpyruvate‑protein phosphotransferase PtsP
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (21/22);  total (22/22)
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

TCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTG  >  NC_000913/2968124‑2968391
                                                                                                                                             |                                                                                                                                
tcttTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCATGGTTTTTTCAg                                                                                                                                   <  2:398466/140‑1 (MQ=255)
     tACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGccccc                                                                                                                             <  2:221027/140‑1 (MQ=255)
           ggggATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTc                                                                                                                        >  2:31730/1‑140 (MQ=255)
            gggATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCg                                                                                                                       >  1:508355/1‑140 (MQ=255)
                tGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCAt                                                                                                                   <  1:31730/140‑1 (MQ=255)
                  tATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCa                                                                                                                 <  1:541277/140‑1 (MQ=255)
                    tttGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAgg                                                                                                               >  2:32531/1‑140 (MQ=255)
                      tGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGgta                                                                                                             >  2:664011/1‑140 (MQ=255)
                        aaGCTGGGGTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGgtagt                                                                                                           >  1:261719/1‑140 (MQ=255)
                        aaGCTGGGGTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGccccccgggt                                                                                                                        <  1:168012/126‑6 (MQ=255)
                        aaGCTGGGGTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCcgggt                                                                                                                        >  2:168012/1‑124 (MQ=255)
                        aaGCTGGGGTGCTTTTGCGCATCTGCAAGGATTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGgtagt                                                                                                           <  2:508355/140‑1 (MQ=255)
                         aGCTGGGGTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGgtagta                                                                                                          >  1:809248/1‑140 (MQ=255)
                              gggTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCAt                                                                                                                   <  2:475499/126‑1 (MQ=255)
                              gggTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCAt                                                                                                                   >  1:475499/1‑126 (MQ=255)
                              gggTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTaacaac                                                                                                     <  2:809248/140‑1 (MQ=255)
                              gggTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTaacaac                                                                                                     >  2:236350/1‑140 (MQ=255)
                               ggTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACg                                                                                                    <  1:668164/140‑1 (MQ=255)
                                   cTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGa                                                                                                <  2:721992/140‑1 (MQ=255)
                                        gcgcATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATg                                                                                           <  1:664011/140‑1 (MQ=255)
                                        gcgcATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATg                                                                                           <  2:261719/140‑1 (MQ=255)
                                         cgcATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACg                                                                                                    <  2:183759/130‑1 (MQ=255)
                                         cgcATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACg                                                                                                    >  1:183759/1‑130 (MQ=255)
                                            aTCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCg                                                                                       <  1:150417/140‑1 (MQ=255)
                                                        tATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGAcc                                                                           >  2:617011/1‑140 (MQ=255)
                                                        tATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGAcc                                                                           >  1:30326/1‑140 (MQ=255)
                                                             gTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACa                                                                      <  2:744114/140‑1 (MQ=255)
                                                             gTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACa                                                                      <  1:604155/140‑1 (MQ=255)
                                                                ccGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGAc                                                                   <  1:748840/140‑1 (MQ=255)
                                                                         tGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATc                                                          <  2:691412/140‑1 (MQ=255)
                                                                            cAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCAt                                                       <  1:166174/140‑1 (MQ=255)
                                                                                 aGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCaa                                                  >  1:103043/1‑140 (MQ=255)
                                                                                 aGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCaa                                                  >  1:444142/1‑140 (MQ=255)
                                                                                                  atcaAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAAcc                                 <  1:32306/140‑1 (MQ=255)
                                                                                                    caAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCa                                >  2:799808/1‑140 (MQ=255)
                                                                                                                      tGCGACCGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCAtt             >  1:43310/1‑140 (MQ=255)
                                                                                                                        cGACCGCGTGGTTTTTTCAGCCCCCCGGGTCACCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCa           >  1:262660/1‑140 (MQ=255)
                                                                                                                           ccGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGcc                  <  2:540089/130‑1 (MQ=255)
                                                                                                                           ccGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGcc                  >  1:540089/1‑130 (MQ=255)
                                                                                                                           ccGCGTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGc        >  1:374179/1‑140 (MQ=255)
                                                                                                                            cgcgTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAcgcg       <  1:474332/140‑1 (MQ=255)
                                                                                                                            cgcgTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCg       <  2:374179/140‑1 (MQ=255)
                                                                                                                             gcgTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCgt      >  1:556195/1‑140 (MQ=255)
                                                                                                                              cgTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCgtg     >  2:124548/1‑140 (MQ=255)
                                                                                                                              cgTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCgtg     >  2:737047/1‑140 (MQ=255)
                                                                                                                              cgTGGTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCgtg     <  2:30326/140‑1 (MQ=255)
                                                                                                                                 ggTTTTTTCAGCCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCgtggtg  <  1:124548/140‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                
TCTTTTACGGAGGGGATGTATTTGAAGCTGGGGTGCTTTTGCGCATCTGCAAGGTTTATCGGTTCCGCCAGCCTGCCAACCAGGCCGACGATCCCTTCATCAAACGCGAGCGTTACAGTGCGACCGCGTGGTTTTTTCAGCCCCCGGGTCGCCATCAGGTAGTAACAACGTCGATCATGATCGGCCAGGTAGACCGAACAGACCTCGGTATCCATCGCAAGACAGATGTCGGTAACCAGAATATTTAACGCCTCATTCAGGCGTGGTG  >  NC_000913/2968124‑2968391

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: