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.0% 97.8 / ‑5.9 33*749* (TAG→TAAptsPPEP‑protein phosphotransferase enzyme I, GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (19/13);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGC  >  NC_000913/2966078‑2966323
                                                                                                              |                                                                                                                                       
gCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCa                                                                                                             >  1:285191/1‑139 (MQ=255)
          tGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACg                                                                                                   <  1:313564/139‑1 (MQ=255)
            cgcgcCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCcgcg                                                                                                      >  2:656042/1‑134 (MQ=255)
            cgcgcCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCcgcg                                                                                                      <  1:656042/134‑1 (MQ=255)
                        ttCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTg                                                                                     <  1:1157735/139‑1 (MQ=255)
                                cTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGa                                                                             >  2:1062184/1‑139 (MQ=255)
                                           cGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTc                                                                  >  2:1673345/1‑139 (MQ=255)
                                                   aGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCa                                                          <  2:285191/139‑1 (MQ=255)
                                                     cacaAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAgctcgc                                                                                                                <  1:1325463/83‑1 (MQ=255)
                                                     cacaAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAgcccgc                                                                                                                >  2:1325463/1‑83 (MQ=255)
                                                      acaAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGtt                                                       >  1:270191/1‑139 (MQ=255)
                                                          ggggACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGa                                                                                                     <  1:1779412/89‑1 (MQ=255)
                                                          ggggACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGa                                                                                                     >  2:1779412/1‑89 (MQ=255)
                                                                ctggctggTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAg                                             >  2:118670/1‑139 (MQ=255)
                                                                    ctggTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                      <  1:184494/126‑1 (MQ=255)
                                                                    ctggTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                      >  2:184494/1‑126 (MQ=255)
                                                                                     ttAAAAGATATGTATATGATCCGCGTTATAACCCTCCGCGAGTCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGc                                                            >  1:1269173/1‑103 (MQ=255)
                                                                                     ttAAAAGATATGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGc                                                            <  2:1269173/103‑1 (MQ=255)
                                                                                             tatGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACtt                                                                   >  1:1056870/1‑88 (MQ=255)
                                                                                             tatGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACtt                                                                   <  2:1056870/88‑1 (MQ=255)
                                                                                              atGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCt                    >  1:1713216/1‑134 (MQ=255)
                                                                                              atGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCt                    <  2:1713216/134‑1 (MQ=255)
                                                                                               tGTATATGATCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCa              >  2:418744/1‑139 (MQ=255)
                                                                                                       aTCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTg                                                                                     <  2:1142564/60‑1 (MQ=255)
                                                                                                       aTCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTg                                                                                     >  1:1142564/1‑60 (MQ=255)
                                                                                                       aTCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCt                    <  2:992413/125‑1 (MQ=255)
                                                                                                       aTCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCt                    >  1:992413/1‑125 (MQ=255)
                                                                                                       aTCCGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCACTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAAt      >  1:1332243/1‑139 (MQ=255)
                                                                                                         ccGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                      >  1:84271/1‑89 (MQ=255)
                                                                                                         ccGCGTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTg                                                      <  2:84271/89‑1 (MQ=255)
                                                                                                          cgcgTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATaa                                                                                           >  2:1309225/1‑51 (MQ=255)
                                                                                                          cgcgTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATaa                                                                                           <  1:1309225/51‑1 (MQ=255)
                                                                                                          cgcgTTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgcg   >  1:1433381/1‑139 (MQ=255)
                                                                                                           gcgtTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTg                                                                                     <  2:1522941/56‑1 (MQ=255)
                                                                                                           gcgtTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTg                                                                                     >  1:1522941/1‑56 (MQ=255)
                                                                                                           gcgtTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgcgc  <  1:418744/139‑1 (MQ=255)
                                                                                                           gcgtTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgcgc  <  1:2022287/139‑1 (MQ=255)
                                                                                                           gcgtTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgcgc  >  2:789566/1‑139 (MQ=255)
                                                                                                              |                                                                                                                                       
GCGATTGAAATGCGCGCCGCAGGTTTCAGGCGCTCCAAAGGTGCGAATAATAGCACAAGGGGACCTGGCTGGTTGCCGGATACCGTTAAAAGATATGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGC  >  NC_000913/2966078‑2966323

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