Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,593,985 Δ1 bp 17.4% coding (68/912 nt) PP_3170 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,593,9850C.17.4% 65.7 / 9.7 23coding (68/912 nt)PP_3170LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (13/6);  new base . (2/2);  total (15/8)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

ATCAGGAACCTTGCATTGCCTGAAATTGCACCCATCTATGGCGGTTCGATCAGCGTCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCGCTGCTGCTGCCAGCCTGGGCCTGACGCGTTCGGCGGTGGGCAAGAGCATCGCCCGCCTGGAAGCCCGCATGGCGGTGCGTCTGCTCAACCGCACTACGCGCAAGCTG  >  NC_002947/3593850‑3594091
                                                                                                                                       |                                                                                                          
aTCAGGAACCTTGCATTGCCTGAAATTGCACCCATCTATGGCGGTTCGATCAGCGTCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCgc                                                                                                            <  1:231380/136‑1 (MQ=255)
    ggAACCTTGCATTGCCTGAAATTGCACCCTTCTATGGCGGTTCGATCAGCGTCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCgctgct                                                                                                        >  1:18035/1‑136 (MQ=255)
    ggAACCTTGCATTGCCTGAAATTGCACCCATCTATGGCGGTTCGATCAGCGTCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCgctgct                                                                                                        >  1:2088/1‑136 (MQ=255)
    ggAACCTTGCATTGCCTGAAATTGCACCCATCTATGGCGGTTCGATCAGCGTCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCgctgct                                                                                                        >  1:307885/1‑136 (MQ=255)
     gAACCTTGCATTGCCTGAAATTGCACCCATCTATGGCGGTTCGATCAGCGTCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCgctgctg                                                                                                       >  1:324305/1‑136 (MQ=255)
             cATTGCCTGAAATTGCACCCATCTATGGCGGTTCGATCAGCGTCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCGCTGCTGCTGCCAGc                                                                                               >  1:203381/1‑136 (MQ=255)
             cATTGCCTGAAATTGCACCCATCTATGGCGGTTCGATCAGCGTCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCGCTGCTGCTGCCAGc                                                                                               >  1:434817/1‑136 (MQ=255)
             cATTGCCTGAAATTGCACCCATCTATGGCGGTTCGATCAGCGTCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCGCTGCTGCTGCCAGc                                                                                               >  1:432101/1‑136 (MQ=255)
                 gCCTGAAATTGCACCCATCTATGGCGGTTCGATCAGCGTCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCGCTGCTGCTGCCAGCCTgg                                                                                           >  1:309427/1‑136 (MQ=255)
                               ccATCTATGGCGGTTCGATCAGCGTCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCGCTGCTGCTGCCAGCCTGGGCCTGACGCGTTcg                                                                             >  2:564785/1‑136 (MQ=255)
                               ccATCTATGGCGGTTCGATCAGCGTCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCGCTGCTGCTGCCAGCCTGGGCCTGACGCGTTcg                                                                             >  2:491665/1‑136 (MQ=255)
                               ccATCTATGGCGGTTCGATCAGCGTCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCGCTGCTGCTGCCAGCCTGGGCCTGACGCGTTcg                                                                             <  1:288646/136‑1 (MQ=255)
                                    tATGGCGGTTCGATCAGCGTCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACC‑‑‑GCTGCTGCCAGCCTGGGCCTGACGCGTTCGGCGGTggg                                                                     >  1:395748/1‑136 (MQ=255)
                                    tATGGCGGTTCGATCAGCGTCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACC‑‑‑GCTGCTGCCAGCCTGGGCCTGACGCGTTCGGCGGTggg                                                                     >  1:398731/1‑136 (MQ=255)
                                    tATGGCGGTTCGATCAGCGTCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGAGCCGCGAAGGCCGGCAGCTTCACAGCTGCTGCTGCCAGCCTGGGCCTGACGCGTTCGGCGGt                                                                        >  1:382808/1‑136 (MQ=255)
                                     aTGGCGGTTCGATCAGCGTCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCGCTGCTGCTGCCAGCCTGGGCCTGACGCGTTCGGCGGTg                                                                       <  1:286439/136‑1 (MQ=255)
                                                       tCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACC‑‑‑GCTGCTGCCAGCCTGGGCCTGACGCGTTCGGCGGTGGGCAAGAGCATCGCCCGCCTg                                                  <  2:398731/136‑1 (MQ=255)
                                                       tCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACC‑‑‑GCTGCTGCCAGCCTGGGCCTGACGCGTTCGGCGGTGGGCAAGAGCATCGCCCGCCTg                                                  <  2:395748/136‑1 (MQ=255)
                                                            aaaaCACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCGCTGCTGCTGCCAGCCTGGGCCTGACGCGTTCGGCGGTGGGCAAGAGCATCGCCCGCCTGGa                                                >  2:425688/1‑136 (MQ=255)
                                                              aaCACCAAGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGGCAGCGCCGCGAAGGCCGGCAGCTTCACCGCTGCTGCTGCCAGCCTGGGCCTGACGCGTTCGGCGGTGGCCAAGAGCATCGCCCGCCTGGAAg                                              <  2:128047/136‑1 (MQ=255)
                                                                  ccATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCGCTGCTGCTGCCAGCCTGGGCCTGACGCGTTCGGCGGTGGGCAAGAGCATCGCCCGCCTGGAAGCCCg                                          >  1:434410/1‑136 (MQ=255)
                                                                  ccATGGAACCCAACCGCATCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGACGGCAGCTTCACCGCTGCTGCTGCCAGCCTGGGCCTGACGCGTTCGGCGGTGGCCAAGAGCATCGCCCCCCTGGAACCCCg                                          >  2:359234/1‑136 (MQ=255)
                                                                        aaCCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGACCGGCAGCTTCACCGCTGCTGCTGCCAGCCTGGGCCTGACGCGTTCGGCGGTGGGCAAGAGCATCGCCCGCCTGGAAGCCCGCATGGc                                    >  2:521451/1‑136 (MQ=255)
                                                                                    tCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCGCTGCTGCTGCCAGCCTGGGCCTGACGCGTTCGGCGGTGGGCAAGAGCATCGCCCGCCTGGAAGCCCGCATGGCGGTGCGTCTGCt                        <  2:307885/136‑1 (MQ=255)
                                                                                     cGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCGCTGCTGCTGCCAGCCTGGGCCTGACGCGTTCGGCGGTGGGCAAGAGCATCGCCCGCCTGGAAGCCCGCATGGCGGTGCGTCTGCTc                       >  1:363534/1‑136 (MQ=255)
                                                                                          aCATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCGCTGCTGCTGCCAGCCTGGGCCTGACGCGTTCGGCGGTGGGCAAGAGCATCGCCCGCCTGGAAGCCCGCATGGCGGTGCGTCTGCTCAACCg                  >  2:487382/1‑136 (MQ=255)
                                                                                                ccgccTTCGTCAGCGCCGCGAAGCCCGGCAGCTTCACCGCTGCTGCTGCCAGCCTGGGCCTGACGCGTTCGGCGGTGGGCAAGAGCATCGCCCGCCTGGAAGCCCGCATGGCGGTGCGTCTGCTCAACCGCACTAc            <  1:521451/136‑1 (MQ=255)
                                                                                                          cAGCGCCGCGAAGGCCGGCAGCTTCACCGCTGCTGCTGCCAGCCTGGGCCTGACGCGTTCGGCGGTGGGCAAGAGCATCGCCCGCCTGGAAGCCCGCATGGCGGTGCGTCTGCTCAACCGCACTACGCGCAAGCTg  <  1:425688/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                          
ATCAGGAACCTTGCATTGCCTGAAATTGCACCCATCTATGGCGGTTCGATCAGCGTCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCGCTGCTGCTGCCAGCCTGGGCCTGACGCGTTCGGCGGTGGGCAAGAGCATCGCCCGCCTGGAAGCCCGCATGGCGGTGCGTCTGCTCAACCGCACTACGCGCAAGCTG  >  NC_002947/3593850‑3594091

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