Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,787,826 Δ1 bp coding (382/834 nt) ppsR → posphoenolpyruvate synthetase regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,787,8260G.94.4% 62.1 / ‑2.3 18coding (382/834 nt)ppsRposphoenolpyruvate synthetase regulatory protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (8/9);  total (8/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGCTTTTGCCAGGATATCGTTCAGGCGCTGGTTGCCCCGCTACAACAAGAGATGAAACTGGATCCAACGCCGATTGCTCATCGTACCCATGGCCTTAACCCTAATAATCTCAATAAATATGATGCGCGCATTGCGGCGATTGATTACACCCTCGCCCACGATGACGGCATTTCGTTGCGCAATCTGGACCAGGCTCAGGTGATCCTGCTCGGTGTTTCTCGCTGTGGTAAAACCCCCACCAGT  >  NC_000913/1787701‑1787945
                                                                                                                             |                                                                                                                       
aaGGCTTTTGCCAGGATATCGTGCAGGCGATGGTTGCACCGCTACAACAAGAGATGAAACTGGATCCAACGCCGATTGCTCATCGTACCCATGGCCTTAACCCTAATAATCTCAATAAATATGATCGCGCATTGCGGCga                                                                                                          <  1:374813/140‑1 (MQ=255)
 aGGCTTTTACCAGGATATCGTTAAGGCGCTGATTGACCCTCTACAACAGGAGATGAAACTGGATCCAACGCCGATTGCTCATCGTACCCATGGCCTTAACCCTAATAATCTCAATAAATATGATCGCGCATTGCGGCgat                                                                                                         <  2:373794/140‑1 (MQ=255)
                        aGGCGCTGGTTGCCCCGCTACAACAAGAGATGAAACTGGATCCAACGCCGATTGCTCATCGTACCCATGGCCTTAACCCTAATAATCTCAATAAATATGATCGCGCATTGCGGCGATTGATTACACCCTCGCCCACGATg                                                                                  >  1:776652/1‑140 (MQ=255)
                         ggCGCTGGTTGCCCCGCTACAACAAGAGATGAAACTGGATCCAACGCCGATTGCTCATCGTACCCATGGCCTTAACCCTAATAATCTCAATAAATATGATCGCGCATTGCGGCGATTGATTACACCCTCGCCCACGATGa                                                                                 >  1:626701/1‑140 (MQ=255)
                             cTGGTTGCCCCGCTACAACAAGAGATGAAACTGGATCCAACGCCGATTGCTCATCGTACCCATGGCCTTAACCCTAATAATCTCAATAAATATGATCGCGCATTGCGGCGATTGCTTACACCCTCGCCCACGATGACGGc                                                                             >  1:304300/1‑140 (MQ=255)
                                                      tGAAACTGGATCCAACGCCGATTGCTCATCGTACCCATGGCCTTAACCCTAATAATCTCAATAAATATGATCGCGCATTGCGGCGATTGATTACACCCTCGCCCACGATGACGGCATTTCGTTGCGCAATCTGGACCAgg                                                    <  1:525160/140‑1 (MQ=255)
                                                       gAAACTGGATCCAACGCCGATTGCTCATCGTACCCATGGCCTTAACCCTAATAATCTCAATAAATATGATCGCGCATTGCGGCGATTGATTACACCCTCGCCCACGATGACGGCATTTCGTTGCGCAATCTGGACCAGGc                                                   <  2:626701/140‑1 (MQ=255)
                                                       gAAACTGGATCCAACGCCGATTGCTCATCGTACCCATGGCCTTAACCCTAATAATCTCAATAAATATGATCGCGCATTGCGGCGATTGATTACACCCTCGCCCACGATGACGGCATTTCGTTGCGCAATCTGGACCAGGc                                                   <  2:62366/140‑1 (MQ=255)
                                                              gATCCAACGCCGATTGCTCATCGTACCCATGGCCTTAACCCTAATAATCTCAATAAATATGATCGCGCATTGCGGCGATTGATTACACCCTCGCCCACGATGACGGCATTTCGTTGCGCAATCTGGACCAGGCTCAGGTg                                            >  1:551267/1‑140 (MQ=255)
                                                               aTCCAACGCCGATTGCTCATCGTACCCATGGCCTTAACCCTAATAATCTCAATAAATATGATCGCGCATTGCGGCGATTGATTACACCCTCGCCCACGATGATGGCATTTAGTTGCGAAATCTGAACCAGGCTCAGGTGa                                           >  1:234225/1‑140 (MQ=255)
                                                                          aTTGCTCATCGTACCCATGGCCTTAACCCTAATAATCTCAATAAATATGATCGCGCATTGCGGCGATTGATTACACCCTCGCCCACGATGACGGCATTTCGTTGCGCAATCTGGACCAGGCTCAgg                                              >  1:455919/1‑126 (MQ=255)
                                                                          aTTGCTCATCGTACCCATGGCCTTAACCCTAATAATCTCAATAAATATGATCGCGCATTGCGGCGATTGATTACACCCTCGCCCACGATGACGGCATTTCGTTGCGCAATCTGGACCAGGCTCAgg                                              <  2:455919/126‑1 (MQ=255)
                                                                                     tACCCATGGCCTTAACCCTAATAATCTCAATAAATATGATCGCGCATTGCGGCGATTGATTACACCCTCGCCCACGATGACGGCATTTCGTTGCGCAATCTGGACCAGGCTCAGGTGATCCTGCTCGGTGTTTCTCGCtg                     <  2:776652/140‑1 (MQ=255)
                                                                                     tACCCATGGCCTTAACCCTAATAATCTCAATAAATATGAGCGCGCATTTCGGCGATTGATTACACCCTCGCCCACGATGACGGCATTTCGTTGCGCAATCTGGACCAGGCTCAGGTGATCCTGCTCGGTGTTTCTCGCtg                     <  2:304300/140‑1 (MQ=255)
                                                                                                  aaCCCTAATAATCTCAATAAATATGATCGCGCATTGCGGCGATTGATTACACCCTCGCCCACGATGACGGCATTTCGTTGCGCAATCTGGACCAGGCTCAGGTGATCCTGCTCGGTGTTTCTCGCTGTGGTAAAAccccc        >  1:194716/1‑140 (MQ=255)
                                                                                                    cccTAATAATCTCAATAAATATGATCGCGCATTGCGGCGATTGATTACACCCTCGCCCACGATGACGGCATTTCGTTGCGCAATCTGGACCAGGCTCAGGTGATCCTGCTCGGTGTTTCTCGCTGTGGTAAAACCCccac      <  2:551267/140‑1 (MQ=255)
                                                                                                    cccTAATAATCTCAATAAATATGATCGCGCACTGCTGCGATTGATTACTCCCTCGCCCACGATGACGGCATTTCGTTGCGCAATCTGGACCAGGCTCAGGTGATCCTGCTCGGTGTTTCTCGCTGTGGTAAAACCCccac      <  2:234225/140‑1 (MQ=255)
                                                                                                        aataatCTCAATAAATATGATCGCGCATTGCGGCGATTGATTACACCCTCGCCCACGATGACGGCATTTCGTTGCGCAATCTGGACCAGGCTCAGGTGATCCTGCTCGGTGTTTCTCGCTGTGGTAAAACCCCCACCAGt  >  1:99598/1‑140 (MQ=255)
                                                                                                                             |                                                                                                                       
AAGGCTTTTGCCAGGATATCGTTCAGGCGCTGGTTGCCCCGCTACAACAAGAGATGAAACTGGATCCAACGCCGATTGCTCATCGTACCCATGGCCTTAACCCTAATAATCTCAATAAATATGATGCGCGCATTGCGGCGATTGATTACACCCTCGCCCACGATGACGGCATTTCGTTGCGCAATCTGGACCAGGCTCAGGTGATCCTGCTCGGTGTTTCTCGCTGTGGTAAAACCCCCACCAGT  >  NC_000913/1787701‑1787945

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