Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,593,986 Δ1 bp 17.4% coding (69/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,9860T.17.4% 57.8 / 8.9 23coding (69/912 nt)PP_3170LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (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 = 7.18e-01

GGAACCTTGCATTGCCTGAAATTGCACCCATCTATGGCGGTTCGATCAGCGTCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCGCTGCTGCTGCCAGCCTGGGCCTGACGCGTTCGGCGGTGGGCAAGAGCATCGCCCGCCTGGAAGCCCGCATGGCGGTGCGTCTGCTCAACCGCACTACGCGCAAGCTG  >  NC_002947/3593854‑3594091
                                                                                                                                    |                                                                                                         
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)
                                                                                                                                    |                                                                                                         
GGAACCTTGCATTGCCTGAAATTGCACCCATCTATGGCGGTTCGATCAGCGTCAGGAAAACACCATGGAACCCAACCGCTTCGGCGACATCGCCGCCTTCGTCAGCGCCGCGAAGGCCGGCAGCTTCACCGCTGCTGCTGCCAGCCTGGGCCTGACGCGTTCGGCGGTGGGCAAGAGCATCGCCCGCCTGGAAGCCCGCATGGCGGTGCGTCTGCTCAACCGCACTACGCGCAAGCTG  >  NC_002947/3593854‑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: