Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 625,065 C→T 9.6% intergenic (‑156/+27) 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,0650CT 9.6% 105.4 / 3.6 42intergenic (‑156/+27)PP_0537/ppatranscriptional regulator/inorganic pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (19/18);  new base T (2/2);  total (22/20)
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.90e-01

AGCGGCGGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTGCT  >  NC_002947/624922‑625213
                                                                                                                                               |                                                                                                                                                     
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                                                                                                                                    <  1:153777/149‑1 (MQ=255)
              aGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCa                                                                                                                                    >  2:153777/1‑149 (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                                                                                                                                               >  1:659071‑M2/1‑88 (MQ=255)
                                                           gCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAgcggg                                                                                                                                               <  2:659071‑M2/93‑6 (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                                                                                                                                      <  2:1080477‑M2/78‑15 (MQ=255)
                                                                                   gACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCTTTGgcggggtttttcat                                                                                                                                      >  1:1080477‑M2/1‑64 (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)
                                                                                                                                               |                                                                                                                                                     
AGCGGCGGGTTGCCAGTAAGGAAATGAACCGGGTGTTGACGGGATTTGTACGAAAAGAGGCTGTTTGTAGGATGCTGGTGTAGGACTGTGTAGGATTTGTTCCCAGCAGTGAGCGCGTTCCATTCAGGCATGAAAAACCCCGCCAAAGCGGGGTTTCCCTGCAGTCGAAATCAGCC‑CTTGTAGGCAGCGACCGATTTGGTGATCGCGGCACGGGCGGCGTCGGCGCCTTCCCAGCCTTCGATCTTGACCCACTTGCCCTTCTCGAGATCTTTGTAGTTGGCAAAGAAGTGCT  >  NC_002947/624922‑625213

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: ‑