Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,872,077 T→A 22.3% D211V (GAC→GTC)  ptsP ← phosphoenolpyruvate‑dependent regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,872,0770TA22.3% 40.6 / 5.7 22D211V (GAC→GTC) ptsPphosphoenolpyruvate‑dependent regulator
Reads supporting (aligned to +/- strand):  ref base T (11/6);  new base A (2/3);  total (13/9)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTCAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAGGCCACG  >  NC_002947/5871936‑5872219
                                                                                                                                             |                                                                                                                                              
gcTTCAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGGGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGACATCGACg                                                                                                                                        <  2:514907/149‑1 (MQ=255)
 cTTCAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGACATCGACgc                                                                                                                                       <  2:286882/150‑2 (MQ=255)
  ttCAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGt                                                                                                                                      >  1:146980/1‑150 (MQ=255)
   tCAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTc                                                                                                                                     <  2:252558/150‑1 (MQ=255)
     aaGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCtt                                                                                                                                   >  1:474689/1‑150 (MQ=255)
                                 gcACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCAAACCTTCCAAGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTGCAGGTcggccggc                                                                                                       >  1:180743/1‑150 (MQ=255)
                                            cGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATcac                                                                                            >  1:604068/1‑150 (MQ=255)
                                                       gggTGGCCAGTTTGGCCGACAGCTTGCGCCTGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGcc                                                                                 >  1:383571/1‑150 (MQ=255)
                                                       gggTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGcc                                                                                 >  1:244416/1‑150 (MQ=255)
                                                           ggCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACTACCGCACGCCCCAc                                                                             <  2:507735/150‑1 (MQ=255)
                                                                gTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCATCGCgg                                                                        >  1:278408/1‑150 (MQ=255)
                                                                                         gtcgCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATAGCCGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCAc                                               >  1:263731/1‑150 (MQ=255)
                                                                                                acacCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAgcgc                                        <  2:604068/150‑1 (MQ=255)
                                                                                                                    gAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGcc                    >  2:62186/1‑150 (MQ=255)
                                                                                                                    gAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGcc                    <  1:84834/150‑1 (MQ=255)
                                                                                                                      aCAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCt                  >  1:510419/1‑150 (MQ=255)
                                                                                                                         gCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGAc               <  1:356595/150‑1 (MQ=255)
                                                                                                                               gTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCa         >  2:334611/1‑150 (MQ=255)
                                                                                                                                   accGTCGATGACATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGCCACCAGGCGAACCCGGCACACCCACGAGGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAGGcc     >  2:10948/3‑150 (MQ=255)
                                                                                                                                     cGTCGATGACATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGaa                                            >  2:77035/1‑109 (MQ=255)
                                                                                                                                     cGTCGATGACATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGaa                                            <  1:77035/109‑1 (MQ=255)
                                                                                                                                      gTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAGGCCACg  <  1:84258/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                              
TCTTCAAGCATCATCAGATAGACGTCGAACAGCGCACGCTCTTCCGGGCGCAGCTGGGTGGCCAGTTTGGCCGACAGCTTGCGCATGTCGTCGCGCACACCTTCCAGGGCGTTCTGGAACAGCTTGAGTTCGGCGTCGATGTCATCGACGGTCTTGTCCGGGACCACTTCCAGGTCGGCCGGCGGCAACATCACCACCGCACGCCCCACCGCGGCACCAGGCGAACCCGGCACACCCACGAAGCGCGCTTCCTGAATACCTTTGCCCTGACGGCCCAGGCCACG  >  NC_002947/5871936‑5872219

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