Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,997,686 C→A 9.5% H124Q (CAC→CAA)  PP_3526 → AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,997,6860CA 9.5% 115.6 / 4.1 42H124Q (CAC→CAA) PP_3526AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (19/19);  new base A (2/2);  total (21/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 = 6.06e-01

TCAGCGACCTGGGGCTGGAACTGGTGGCCACCGAGCTGTCCGCTGCGGCGCTGAAGGAGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTGGACGACTACGGTTGCAGCATCCACCCCGAGCAACGTGCCAGCCTGTCGGAGCGCAGCCCGCAAACCCGCATCACCCCGGC  >  NC_002947/3997541‑3997826
                                                                                                                                                 |                                                                                                                                            
tCAGCGACCTGGGGCTGGAACTGGTGGCCACCGAGCTGTCCGCTGCGGCGCTGAAGGAGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACgg                                                                                                                                            >  8:254504/1‑148 (MQ=255)
tCAGCGACCTGGGGCTGGAACTGGTGGCCACCGAGCTGTCCGCTGCGGCGCTGAAGGAGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACgg                                                                                                                                            >  5:303674/1‑148 (MQ=255)
tCAGCGACCTGGGGCTGGAACTGGTGGCCACCGAGCTGTCCGCTGCGGCGCTGAAGGAGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCAAGGc                                                                                                                                           >  5:92973/1‑149 (MQ=255)
   gCGACCTGGGGCTGGAACTGGTGGCCACCGAGCTGTCCGCTGCGGCGCTGAAGTAGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCTCCGCCCACGTCATg                                                                                                                                        <  3:214882/149‑1 (MQ=255)
   gCGACCTGGGGCTGGAACTGGTGGCCACCGAGCTGTCCGCTGCGGCGCTGAAGGAGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATg                                                                                                                                        <  6:223884/149‑1 (MQ=255)
             gCTGGAACTGGTGGCCACCGAGCTGTCCGCTGCGGCGCTGAAGGAGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCAAGGCCATGCCGTGggcg                                                                                                                              <  7:260902/149‑1 (MQ=255)
                 gAACTGGTGGCCACCGAGCTGTCCGCTGCGGCGCTGAAGGAGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCAAGGCCATGCCGTGggcggc                                                                                                                            <  6:92973/147‑1 (MQ=255)
                    cTGGTGGTGACCGAGCTGTCCGCTGCGGCGCTGAAGGAGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCAGTCAAGTACCCGGAACTCGATCGCCTGCTATACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTtgtg                                                                                                                       <  7:221705/149‑1 (MQ=255)
                       gtggCCACCGAGCTGTCCGCTGCGGCGCTGAAGGAGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGaa                                                                                                                    <  7:254504/149‑1 (MQ=255)
                            cACCGAGCTGTCCGCTGCGGCGCTGAAGGAGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGgcg                                                                                                               >  8:104201/1‑149 (MQ=255)
                                 aGCTGTCCGCTGCGGCGCAGAAGGAGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCggg                                                                                                          >  4:23268/1‑149 (MQ=255)
                                      tCCGCTGCGGCGCTGAAGGAGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCt                                                                                                     >  2:37855/1‑149 (MQ=255)
                                        cGCTGCGGCGCTGAAGGAGCTGGACCTTCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTg                                                                                                   <  4:272703/149‑1 (MQ=255)
                                        cGCTGCGGCGCTGAAGGAGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCtt                                                                                                    <  1:37855/148‑1 (MQ=255)
                                                 gcTGAAGGAGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCg                                                                                          >  7:87048/1‑149 (MQ=255)
                                                     aaGGAGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGt                                                                                      >  3:64988/1‑149 (MQ=255)
                                                     aaGGAGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTAACCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGt                                                                                      >  5:152406/1‑149 (MQ=255)
                                                      aGGAGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTg                                                                                     <  6:303674/149‑1 (MQ=255)
                                                         aGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTg                                                                                  <  7:104201/149‑1 (MQ=255)
                                                         aGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTg                                                                                  <  5:40869/149‑1 (MQ=255)
                                                         aGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTg                                                                                  <  5:45549/149‑1 (MQ=255)
                                                           cTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTGga                                                                                <  8:87048/149‑1 (MQ=255)
                                                             ggACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTGgac                                                                               >  2:159795/1‑148 (MQ=255)
                                                              gACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTGgac                                                                               <  4:64988/147‑1 (MQ=255)
                                                                    ctgGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGc                                                                                    <  7:245516/136‑1 (MQ=255)
                                                                    ctgGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGc                                                                                    >  8:245516/1‑136 (MQ=255)
                                                                        tggtCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGCACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTGGACGACTACGGTTgc                                                                   <  6:152406/149‑1 (MQ=255)
                                                                              gcggcgGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCt                                                                                                            <  6:20965/102‑1 (MQ=255)
                                                                              gcggcgGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCt                                                                                                            >  5:20965/1‑102 (MQ=255)
                                                                                                   tCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTGga                                                                                <  7:98109/109‑1 (MQ=255)
                                                                                                   tCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTGga                                                                                >  8:98109/1‑109 (MQ=255)
                                                                                                     aaGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTGGACGACTACgg                                                                       <  6:47042/116‑1 (MQ=255)
                                                                                                     aaGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTGGACGACTACgg                                                                       >  5:47042/1‑116 (MQ=255)
                                                                                                          cccGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTGGACGACTACGGTTGCAGCATCCACCCCGAGCAACGTGCCAGCCTGTCg                                  >  5:225013/1‑148 (MQ=255)
                                                                                                            cGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGTTGCTGGACGACTACGGTTGCAGCATCCACCCCGAGCAACGTGCCAGCCTGTCGGAg                               >  2:10409/1‑149 (MQ=255)
                                                                                                                    gATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTGGACGACTACgg                                                                       >  6:238450/1‑101 (MQ=255)
                                                                                                                    gATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTGGACGACTACgg                                                                       <  5:238450/101‑1 (MQ=255)
                                                                                                                    gATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTGGACGACTACGGTTGCAGCATCCACCCCGAGCAACGTGCCAGCCTGTCGGAGCGCAGccc                       >  7:247747/1‑149 (MQ=255)
                                                                                                                               aaaCGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTGGACGACTACGGTTGCAGCATCCACCCCGAGCAACGTGCCAGCCTGTCGGAGCGCAGCCCGCAAACCCGCa            >  2:170774/1‑149 (MQ=255)
                                                                                                                                aaCGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTGGACGACTACGGTTGCAGCATCCACCCCGAGCAACGTGCCAGCCTGTCGGAGCGCAGCCCGCAAACCCGCAt           >  3:124659/1‑149 (MQ=255)
                                                                                                                                        cgccgcCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTGGACGACTACGGTTGCAGCATCCACCCCGAGCAACGTGCCAGCCTGTCGGAGCGCAGCCCGCAAACCCGCATCACCCCg    <  3:259873/148‑1 (MQ=255)
                                                                                                                                           cgccCAAGGCCATGCCGTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTGGACGACTACGGTTGCAGCATCCACCCCGAGCAACGTGCCAGCCTGTCGGAGCGCAGCCCACAAACCCGCATCACCCCGGc  >  8:48040/1‑147 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
TCAGCGACCTGGGGCTGGAACTGGTGGCCACCGAGCTGTCCGCTGCGGCGCTGAAGGAGCTGGACCTGCTGGTGGTCTGCGGCGGCTTGCGCACACCACTCAAGTACCCGGAACTCGATCGCCTGCTAAACGATTGCGCCGCCCACGGCATGGCCTTGGGCGGCTTGTGGAACGGCGCCTGGTTCCTTGGCCGTGCCGGGGTGCTGGACGACTACGGTTGCAGCATCCACCCCGAGCAACGTGCCAGCCTGTCGGAGCGCAGCCCGCAAACCCGCATCACCCCGGC  >  NC_002947/3997541‑3997826

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