Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 625,068 A→G 9.3% intergenic (‑159/+24) PP_0537 ← / ← ppa transcriptional regulator/inorganic pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947625,0680AG 9.3% 117.1 / 2.9 43intergenic (‑159/+24)PP_0537/ppatranscriptional regulator/inorganic pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (20/19);  new base G (2/2);  total (22/21)
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.03e-01

AGCGGCGGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTGCTCGA  >  NC_002947/624922‑625216
                                                                                                                                                  |                                                                                                                                                     
aGCGGCGGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAgcg                                                                                                                                                    <  1:230293‑M2/150‑4 (MQ=255)
 gcggcgGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAgcgg                                                                                                                                                   >  1:209222‑M2/1‑146 (MQ=255)
    gcggGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAgcggggg                                                                                                                                                <  2:537189‑M2/150‑8 (MQ=255)
      gggTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAgcggggttt                                                                                                                                              >  1:791633‑M2/1‑141 (MQ=255)
             cAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCa                                                                                                                                       >  2:597522/1‑150 (MQ=255)
             cAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCa                                                                                                                                       >  1:683400/1‑150 (MQ=255)
              aGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAgcgg                                                                                                                                                   >  2:563578‑M2/1‑133 (MQ=255)
              aGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAgcgg                                                                                                                                                   <  1:563578‑M2/137‑5 (MQ=255)
              aGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCa                                                                                                                                       >  2:153777/1‑149 (MQ=255)
              aGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCa                                                                                                                                       <  1:153777/149‑1 (MQ=255)
              aGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCAAGCAGTGAGCGCGTTCCATTTAGGCATGAAAAACCCCGCAAAAGCGGGGTTTCCCTgaaa                                                                                                                                      >  2:552841/1‑147 (MQ=255)
                      aaTGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATc                                                                                                                              >  2:7584/1‑150 (MQ=255)
                       aTGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCa                                                                                                                             <  1:417119/150‑1 (MQ=255)
                               ggTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGt                                                                                                                    <  2:473811/150‑1 (MQ=255)
                                 tgtTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAg                                                                                                                  <  1:7584/150‑1 (MQ=255)
                                      aCGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGc                                                                                                             <  2:683400/150‑1 (MQ=255)
                                               gTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGCGACCGAttt                                                                                                    <  1:552841/150‑1 (MQ=255)
                                                tACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGgcggggtttttcatgcctgaa                                                                                                                                  >  1:853737‑M2/1‑99 (MQ=255)
                                                tACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGgcggggtttttcatgcctgaa                                                                                                                                  <  2:853737‑M2/120‑22 (MQ=255)
                                                tACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGCGACCGATTTg                                                                                                   <  2:992228/150‑1 (MQ=255)
                                                           gCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAgcggg                                                                                                                                                  <  2:659071‑M2/93‑6 (MQ=255)
                                                           gCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAgcggg                                                                                                                                                  >  1:659071‑M2/1‑88 (MQ=255)
                                                                 tGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGggc                                                                                  <  2:355483/150‑1 (MQ=255)
                                                                            ggTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGgcgc                                                                       <  2:439492/150‑1 (MQ=255)
                                                                               gtAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGAGTTCCATTCAGTCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCtt                                                                    >  2:567650/1‑150 (MQ=255)
                                                                                   gACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGgcggggtttttcat                                                                                                                                         >  1:1080477‑M2/1‑64 (MQ=255)
                                                                                   gACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGgcggggtttttcat                                                                                                                                         <  2:1080477‑M2/78‑15 (MQ=255)
                                                                                   gACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCa                                                                <  1:901210/150‑1 (MQ=255)
                                                                                          tAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCg                                                         >  2:88398/1‑150 (MQ=255)
                                                                                          tAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTAGGTA‑CGGCCACCTTTATGGTGCGACCGATTTGGTGCTAGCGCGGCGGTCGGCGTCGTCTCCTTCCCAGCCTTCg                                                         >  1:23265/1‑150 (MQ=255)
                                                                                            ggATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGAt                                                       >  2:3038/1‑150 (MQ=255)
                                                                                             gATTTGTTCCCAGCAGTGCGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGCAAGCAGCTACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATc                                                      >  2:494916/1‑150 (MQ=255)
                                                                                             gATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATc                                                      <  2:230293/150‑1 (MQ=255)
                                                                                                               aGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTTGCCAGGCTTG‑CGTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTtc                                    >  2:1138351/1‑150 (MQ=255)
                                                                                                                    gTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCgag                               >  2:1207166/1‑150 (MQ=255)
                                                                                                                          ttCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGGAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCttt                         <  2:23265/150‑1 (MQ=255)
                                                                                                                               gCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCttt                         >  2:1143985/1‑145 (MQ=255)
                                                                                                                               gCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCttt                         <  1:1143985/145‑1 (MQ=255)
                                                                                                                                    aaaaaCCCCGCCAAAgcggggtttctcatgcctgaatggaacgcgctcactgctgggaacaaatcctacacagtcctacaccagcatcctacaaacagcctcttttcgtacaaatcccgtcaacacccggttcatttccttactggcaac                <  2:452002‑M2/150‑136 (MQ=255)
                                                                                                                                        aCCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAaaga           >  2:1056900/1‑150 (MQ=255)
                                                                                                                                          cccGCCAAAGCGGGGTTTCCCTGCAGTNGAAATCAGCC‑CTNGTAGGCNGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAaagaag         >  1:750673/1‑150 (MQ=255)
                                                                                                                                          cccGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAaagaag         >  1:1277175/1‑150 (MQ=255)
                                                                                                                                              ccAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTGCt     <  1:1238925/150‑1 (MQ=255)
                                                                                                                                                 aaGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTGCTCGa  <  1:537189/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AGCGGCGGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTGCTCGA  >  NC_002947/624922‑625216

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