Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,041,937 T→G 11.4% intergenic (+22/+93) appA → / ← etk periplasmic phosphoanhydride phosphatase/multiple inositol‑polyphosphate phosphatase/protein‑tyrosine kinase Etk

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,041,9370TG11.4% 87.2 / 3.1 43intergenic (+22/+93)appA/etkperiplasmic phosphoanhydride phosphatase/multiple inositol‑polyphosphate phosphatase/protein‑tyrosine kinase Etk
Reads supporting (aligned to +/- strand):  ref base T (20/16);  new base G (3/2);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01

GTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAA  >  NC_000913/1041805‑1042075
                                                                                                                                    |                                                                                                                                          
gTGAGACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGGTAACTGAATg                                                                                                                                  <  2:739945/143‑1 (MQ=255)
gTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTGATATACCTGAGTg                                                                                                                                  >  2:615980/1‑143 (MQ=255)
gTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGGTACCTGAATg                                                                                                                                  >  2:786403/1‑143 (MQ=255)
gTGAAACTGACCCTAGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGCTGGCAGGTCTTAAGAGATTCTTGTATGGAGCACGCACATCGGCGCGCAGTTTGTCATGCATAAAAAGGAGCATTCCGTGACCTGACTg                                                                                                                                  >  2:475244/1‑143 (MQ=255)
   aaaCTGGCCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGctc                                                                                                                               <  2:688546/143‑1 (MQ=255)
        gACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGGTAACTGAATGctcttttt                                                                                                                          <  1:426301/143‑5 (MQ=255)
                       tgAAGAGCGAAATTGGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACgaa                                                                                                           <  2:212122/143‑1 (MQ=255)
                       tgAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACgaa                                                                                                           <  2:600877/143‑1 (MQ=255)
                       tgAAGAGCGAAATGCGCAGGGCATGGGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACgaa                                                                                                           <  1:86201/143‑1 (MQ=255)
                             gCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGGTACCTGAATGCTCTGAGGCTGATGACAAACGCACAACTg                                                                                                     >  2:642808/1‑143 (MQ=255)
                              cGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCGAATCGTGACTGACGCCCGCATACCGGCGTGCAGCTTGTAATGCATACGCAACAGCATTCCGACACCTGAATGCTCTGAGGCTGATCACCACCGAAGACCTGt                                                                                                    >  2:514944/1‑143 (MQ=255)
                              cGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGCATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAACCGAAGAACTGt                                                                                                    >  2:359962/1‑143 (MQ=255)
                              cGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGt                                                                                                    >  2:652792/1‑143 (MQ=255)
                              cGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAAGCGTGAATGACGCACGCATACCGGCGTTCAGTTTGTCATGCATAAAACAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATAACAAACGAAGACCTGt                                                                                                    >  2:9251/1‑143 (MQ=255)
                                   tGCGCCGGGCATGTGTCTGTTGGCAGGTCTTAAGCAAACTGTGAATGACGCTGGTGTACCGGCGTGCAGGTTGTACTGCATAAAAACGAGCATTCAGTTACCTGAACTGGCTGAGGCTGATGACAAACGAAGAACTGTCTAAt                                                                                               <  2:597830/143‑1 (MQ=255)
                                       cAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGt                                                                                           >  1:630318/1‑143 (MQ=255)
                                        aGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTa                                                                                          >  2:277524/1‑143 (MQ=255)
                                           gCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGCCAAACGAAGAACTGTCTAATGCGTAGAc                                                                                       >  1:710250/1‑143 (MQ=255)
                                           gCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGAc                                                                                       >  2:63559/1‑143 (MQ=255)
                                              tgtgTTCGTTGGCAGGTTTTACGCAAAGCGTGAATGAAGCACGCATAGCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCgg                                                                                    <  2:81284/143‑1 (MQ=255)
                                                  ttCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGaaaa                                                                                >  1:633886/1‑143 (MQ=255)
                                                   tCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAg                                                                               >  1:125130/1‑143 (MQ=255)
                                                                ttACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAGTGCGTAGACCGGAAAAGGCGTTCAcgccgc                                                                  >  1:296291/1‑143 (MQ=255)
                                                                      aaaTCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCgg                                                            >  1:354608/1‑143 (MQ=255)
                                                                       aaTCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAAGGCATAAAAAAGAGCACTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGc                                                           <  1:665913/143‑1 (MQ=255)
                                                                              aatgaaGCACGCATACCGGCGTGCAGTTTGTAAGGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTc                                                    <  2:397407/143‑1 (MQ=255)
                                                                                 gaaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACttcccgt                                                 >  1:248538/1‑138 (MQ=255)
                                                                                  aaGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGtt                                                >  2:95339/1‑143 (MQ=255)
                                                                                   aGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGttt                                               <  2:134758/143‑1 (MQ=255)
                                                                                            aCCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCAATTGCTGATGCGGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTtc                                      >  1:179382/1‑143 (MQ=255)
                                                                                                 cGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGa                                 <  1:765586/143‑1 (MQ=255)
                                                                                                 cGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGa                                 <  1:219949/143‑1 (MQ=255)
                                                                                                    gCAGTTTGCCAGGCATAAAGAACGGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAATGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTa                              <  1:71321/143‑1 (MQ=255)
                                                                                                            tAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAAc                      >  1:692142/1‑143 (MQ=255)
                                                                                                               tGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGt                   <  1:289883/143‑1 (MQ=255)
                                                                                                                    gaagaaGAGAATTCAGTGACCTGAATGCTCTGAGTATGATGACAAACTAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAg              <  2:29536/139‑1 (MQ=255)
                                                                                                                     aaaaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTCCAGTTTTACTCTTTCTCGGAGGAACCATAACCGTAATAgg             >  2:286287/1‑142 (MQ=255)
                                                                                                                       aaaGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTta           >  1:476620/1‑143 (MQ=255)
                                                                                                                             cATTCAGGTAACTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCg     >  1:247925/1‑143 (MQ=255)
                                                                                                                              atacAGTCACCGGAAGGCTCTGAGGCGGATGACAAACGAAGAACGGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGt    <  1:775369/140‑1 (MQ=255)
                                                                                                                              aTTCAGTTCCCTGAATGCGCTGAGGCTGATGACAAATGAAGAGGTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGt    <  2:163796/143‑1 (MQ=255)
                                                                                                                              aTTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTGCTCGGAGTAACTATAACCGTAATAGTTATAGCCGt    >  2:47625/1‑143 (MQ=255)
                                                                                                                                tCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGCAGAACTGTCTAATGCGTAGACCGGCAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAAATATAACCGTAATAGTTCTAGCAGtca  >  2:620617/1‑141 (MQ=255)
                                                                                                                                    |                                                                                                                                          
GTGAAACTGACCCTGGCAGGATGTGAAGAGCGAAATGCGCAGGGCATGTGTTCGTTGGCAGGTTTTACGCAAATCGTGAATGAAGCACGCATACCGGCGTGCAGTTTGTAATGCATAAAAAAGAGCATTCAGTTACCTGAATGCTCTGAGGCTGATGACAAACGAAGAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAA  >  NC_000913/1041805‑1042075

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