Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,966,188 C→T *749* (TAG→TAA ptsP ← PEP‑protein phosphotransferase enzyme I, GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,966,1880CT97.4% 120.2 / ‑3.4 39coding (2247/2247 nt)ptsPPEP‑protein phosphotransferase enzyme I, GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (21/17);  minor base . (1/0);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

TCTTGTTATTTCGCTGAAAGAGAACAGCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGC  >  NC_000913/2966052‑2966323
                                                                                                                                        |                                                                                                                                       
tCTTGTTATTTCGCTGAAAGAGAACAGCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTta                                                                                                                                       <  1:1125107/139‑1 (MQ=255)
            gCTGAAAGAGAACAGCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCcgcg                                                                                                                           >  2:1008751/1‑139 (MQ=255)
                  agagAACAGCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAg                                                                                                                     <  2:745873/139‑1 (MQ=255)
                           cGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTATAACCCTCCGCGAATCAGCCCGCCCATg                                                                                                           >  1:1187280/1‑139 (MQ=255)
                                 aaaTGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCcgcg                                                                                                      <  2:16500/139‑1 (MQ=255)
                                           cgcAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATg                                                                                                           >  1:1020434/1‑124 (MQ=255)
                                           cgcAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATg                                                                                                           <  2:1020434/124‑1 (MQ=255)
                                                          cTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGa                                                                             >  1:98406/1‑139 (MQ=255)
                                                                 ggTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGaa                                                                      <  1:1008751/139‑1 (MQ=255)
                                                                     cGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTc                                                                  >  2:1287562/1‑139 (MQ=255)
                                                                                 caAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGc                                                                                       <  1:525691/106‑1 (MQ=255)
                                                                                 caAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGc                                                                                       >  2:525691/1‑106 (MQ=255)
                                                                                    ggggACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTgcgc                                                   >  1:11373/1‑139 (MQ=255)
                                                                                    ggggACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                      >  1:1287638/1‑136 (MQ=255)
                                                                                    ggggACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                      <  2:1287638/136‑1 (MQ=255)
                                                                                    ggggACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCACGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTgcgc                                                   >  1:32811/1‑139 (MQ=255)
                                                                                      ggACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCtt                                                 <  1:1287562/139‑1 (MQ=255)
                                                                                       gACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTc                                                >  2:532818/1‑139 (MQ=255)
                                                                                              ctggTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTa                                         >  2:744248/1‑139 (MQ=255)
                                                                                                ggTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACg                                       <  2:863902/139‑1 (MQ=255)
                                                                                                            ccGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGAtt                           >  2:384911/1‑139 (MQ=255)
                                                                                                               ttAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGc                                                                                       >  1:209054/1‑76 (MQ=255)
                                                                                                               ttAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGc                                                                                       <  2:209054/76‑1 (MQ=255)
                                                                                                                   aaGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCt                    <  2:11373/139‑1 (MQ=255)
                                                                                                                       tatGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACtt                                                                   >  1:831324/1‑88 (MQ=255)
                                                                                                                       tatGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACtt                                                                   <  2:831324/88‑1 (MQ=255)
                                                                                                                         tGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCa              >  1:291486/1‑139 (MQ=255)
                                                                                                                              atGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGc                                                                                                  >  2:391588/1‑50 (MQ=255)
                                                                                                                              atGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGc                                                                                                  <  1:391588/50‑1 (MQ=255)
                                                                                                                              atGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATa                                                                                            >  2:874462/1‑56 (MQ=255)
                                                                                                                              atGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATa                                                                                            <  1:874462/56‑1 (MQ=255)
                                                                                                                               tGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                      >  1:1031463/1‑93 (MQ=255)
                                                                                                                               tGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                      <  2:1031463/93‑1 (MQ=255)
                                                                                                                               tGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCa        >  1:672187/1‑139 (MQ=255)
                                                                                                                                gATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCaa       >  2:360418/1‑139 (MQ=255)
                                                                                                                                   ccGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATaaa                                                                                          >  2:214978/1‑53 (MQ=255)
                                                                                                                                   ccGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATaaa                                                                                          <  1:214978/53‑1 (MQ=255)
                                                                                                                                    cgcgTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCtt                                                 >  1:768830/1‑93 (MQ=255)
                                                                                                                                    cgcgTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCtt                                                 <  2:768830/93‑1 (MQ=255)
                                                                                                                                     gcgtTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGa                                                                                                     <  2:676282/40‑1 (MQ=255)
                                                                                                                                     gcgtTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGc                                                                                                  <  1:983329/43‑1 (MQ=255)
                                                                                                                                     gcgtTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGc                                                                                                  >  2:983329/1‑43 (MQ=255)
                                                                                                                                     gcgtTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGGTGCGCAAGATTTTCTGCTTCGGCATAATCAATgcgc  >  1:824918/1‑139 (MQ=255)
                                                                                                                                     gcgtTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgcgc  >  2:735741/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                       
TCTTGTTATTTCGCTGAAAGAGAACAGCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGC  >  NC_000913/2966052‑2966323

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