Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,967,707 T→A D243V (GAT→GTT)  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,967,7070TA97.1% 232.3 / ‑1.2 69D243V (GAT→GTT) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (39/28);  total (40/29)
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.46e-01

CCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTC  >  NC_000913/2967612‑2967803
                                                                                               |                                                                                                
ccGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCAACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGaa                                                                                              >  2:1374019/1‑100 (MQ=255)
 cGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGaaa                                                                                             >  2:261942/1‑100 (MQ=255)
  gTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGaaaa                                                                                            >  2:1573763/1‑100 (MQ=255)
       tACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAgcc                                                                                       >  1:1286982/1‑100 (MQ=255)
       tACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAgc                                                                                        >  2:667606/1‑99 (MQ=255)
        aCCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAgccg                                                                                      >  2:518773/1‑100 (MQ=255)
         ccGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAgccgc                                                                                     <  1:1654843/100‑1 (MQ=255)
         ccGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAgccgc                                                                                     >  2:713714/1‑100 (MQ=255)
         ccGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGACCGAAAGCAGGTGCGAGTAAAGAACGAAAATacac                                                                                       >  1:1463009/1‑95 (MQ=255)
          cGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAgccgcc                                                                                    >  2:591762/1‑100 (MQ=255)
                     ccACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTtg                                                                         <  1:627978/100‑1 (MQ=255)
                     ccACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTtg                                                                         >  2:1686616/1‑100 (MQ=255)
                      cACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTtgt                                                                        <  1:775386/100‑1 (MQ=255)
                      cACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTtgt                                                                        >  2:355897/1‑100 (MQ=255)
                       aCCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTtgtg                                                                       <  1:531198/100‑1 (MQ=255)
                            gCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACcg                                                                  >  2:553578/1‑100 (MQ=255)
                             ccTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACcgg                                                                 <  1:268630/100‑1 (MQ=255)
                              cTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACcggc                                                                >  2:1652020/1‑100 (MQ=255)
                                ttATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACcggcgg                                                              <  2:562565/100‑1 (MQ=255)
                                ttATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACcggcgg                                                              >  1:782179/1‑100 (MQ=255)
                                 tATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCggcggc                                                             <  1:348393/100‑1 (MQ=255)
                                  aTCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGACAGCAGGTGCGAGTAACGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACcgg                                                                 >  1:1706877/1‑95 (MQ=255)
                                  aTCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCa                                                            >  1:311952/1‑100 (MQ=255)
                                      aCCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGc                                                        >  2:837003/1‑100 (MQ=255)
                                      aCCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGACGGCAAAGc                                                        >  2:743982/1‑100 (MQ=255)
                                       ccTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCg                                                       >  2:1095754/1‑100 (MQ=255)
                                       ccTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCg                                                       <  2:202138/100‑1 (MQ=255)
                                       ccTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCg                                                       <  1:1480469/100‑1 (MQ=255)
                                        cTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGt                                                      >  1:1414058/1‑100 (MQ=255)
                                        cTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGt                                                      <  2:1189645/100‑1 (MQ=255)
                                          cGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGttt                                                    <  1:1014610/100‑1 (MQ=255)
                                          cGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCAGCGGCAAAGCGttt                                                    >  2:1467004/1‑100 (MQ=255)
                                            gCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGc                                                  <  2:1015309/100‑1 (MQ=255)
                                              aaaCAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTg                                                <  1:1283852/100‑1 (MQ=255)
                                               aaCAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGt                                               >  2:924142/1‑100 (MQ=255)
                                                aCAATTCGCGACGCAGCCGGTTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTa                                              <  2:1010779/100‑1 (MQ=255)
                                                aCAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTa                                              >  1:1519139/1‑100 (MQ=255)
                                                 cAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAg                                             >  1:492070/1‑100 (MQ=255)
                                                  aaTTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAgc                                            >  1:1113138/1‑100 (MQ=255)
                                                  aaTTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAgc                                            >  1:1348469/1‑100 (MQ=255)
                                                   aTTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAgcg                                           >  1:748462/1‑100 (MQ=255)
                                                    ttCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAgcgg                                          <  2:1436784/100‑1 (MQ=255)
                                                     tCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTCgcggc                                         >  1:1605573/1‑100 (MQ=255)
                                                        cgACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGaa                                      <  2:392936/100‑1 (MQ=255)
                                                        cgACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGaa                                      >  2:319430/1‑100 (MQ=255)
                                                           cGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGAGTTTGCTGTAGCGGCGAAAc                                    >  1:1656718/1‑99 (MQ=255)
                                                            gCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTc                                  <  1:850798/100‑1 (MQ=255)
                                                             cAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGaaa                                     >  1:1740062/1‑96 (MQ=255)
                                                                ccGGGTGTCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGttt                              <  1:324423/100‑1 (MQ=255)
                                                                   ggTGTCCGAAAGCAGGTGCGAGTAAAGAACGCAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGcc                           <  1:82140/100‑1 (MQ=255)
                                                                    gtgtCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCg                          <  2:1739859/100‑1 (MQ=255)
                                                                     tgtCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGc                         >  2:1354121/1‑100 (MQ=255)
                                                                      gtCCGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGct                        <  2:336747/100‑1 (MQ=255)
                                                                        ccGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGcttc                      >  1:314572/1‑100 (MQ=255)
                                                                         cGAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGcttct                     >  1:154217/1‑100 (MQ=255)
                                                                          gAAAGCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGcttctt                    >  2:1511736/1‑100 (MQ=255)
                                                                             aGCATGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCa                 <  1:1722046/100‑1 (MQ=255)
                                                                              gCAGGTGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAg                <  1:1652009/100‑1 (MQ=255)
                                                                                   tGCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCAAAACTCGTTTGCCGCTTCTTCCAGCGccc           >  2:39611/1‑100 (MQ=255)
                                                                                    gCGAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGcccc          <  1:1275943/100‑1 (MQ=255)
                                                                                      gAGTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCgg        <  1:261941/100‑1 (MQ=255)
                                                                                        gTAAATAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTGGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTc      <  1:1672763/100‑1 (MQ=255)
                                                                                        gTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTc      <  1:836998/100‑1 (MQ=255)
                                                                                        gTAAAGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTc      >  2:643927/1‑100 (MQ=255)
                                                                                          aaaGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAg    >  2:1063000/1‑100 (MQ=255)
                                                                                          aaaGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTACCGCTTCTTCCAGCGCCCCGGTCAg    >  2:140570/1‑100 (MQ=255)
                                                                                           aaGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGt   <  1:1245964/100‑1 (MQ=255)
                                                                                           aaGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGt   <  2:770133/100‑1 (MQ=255)
                                                                                            aGAACGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGAGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTc  >  2:103563/1‑100 (MQ=255)
                                                                                               |                                                                                                
CCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTC  >  NC_000913/2967612‑2967803

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