Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,170,445 T→G 8.9% L56R (CTG→CGG)  PP_3669 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,170,4450TG 8.9% 132.4 / 5.6 45L56R (CTG→CGG) PP_3669LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (25/16);  new base G (2/2);  total (27/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

ACATGAACCTGCTGGTGCTGTTCGAAGCCCTGATGATCGAACGCAACCTGACCCGCGTGGGTGAAAAGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTCGACCCGGCCAACAGCTGCAACGTGTTCCGCCTGGGC  >  NC_002947/4170307‑4170593
                                                                                                                                          |                                                                                                                                                    
aCATGAACCTGCTGGTGCTGTTCGAAGCCCTGATGATCGAACGCAACCTGACCCGCGTGGGTGAAAAGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATc                                                                                                                                                  <  6:29473/143‑1 (MQ=255)
aCATGAACCTGCTGGTGCTGTTCGAAGCCCTGATGATCGAACGCAACCTGACCCGCGTGGGTGAAAAGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATc                                                                                                                                                  >  5:29473/1‑143 (MQ=255)
  aTGAACCTGCTGGTGCTGTTCGAAGCCCTGATGATCGCACGCAACCTGACCCGCGTGGGTGAAAAGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCGGATCAGCAGCgg                                                                                                                                          <  4:38834/149‑1 (MQ=255)
           ctgGTGCTGTTCGAAGCCCTGATGATCGAACGCAACCTGACCCGCGTGGGTGAAAAGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCa                                                                                                                                  <  1:300299/148‑1 (MQ=255)
                  tGTTCGAAGCCCTGATGATCGAACGCAACCTGACCCGCGTGGGTGAAAAGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCGGATCAGCAGCgg                                                                                                                                          <  1:70290/133‑1 (MQ=255)
                  tGTTCGAAGCCCTGATGATCGAACGCAACCTGACCCGCGTGGGTGAAAAGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCGGATCAGCAGCgg                                                                                                                                          >  2:70290/1‑133 (MQ=255)
                  tGTTCGAAGCCCTGATGATCGAACGCAACCTGACCCGCGTGGGTGAAAAGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGAGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAaccc                                                                                                                          >  2:35061/1‑149 (MQ=255)
                           cccTGATGATCGAACGCAACCTGACCCGCGTGGGTGAAAAGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACgg                                                                                                                 <  5:265798/149‑1 (MQ=255)
                               gatgatCGAACGCAACCTGACCCGCGTGGGTGAAAAGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCa                                                                                                             >  1:151203/1‑149 (MQ=255)
                                           cAACCTGACCCGCGTGGGTGAAAAGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTg                                                                                                 <  2:151203/149‑1 (MQ=255)
                                            aaCCTGACCCGCGTGGGTGAAAAGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGc                                                                                                >  7:152813/1‑149 (MQ=255)
                                             aCCTGACCCGCGTGGGTGAAAAGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCt                                                                                               <  7:158333/149‑1 (MQ=255)
                                                tGACCCGCGTGGGTGAAAAGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGaa                                                                                            >  1:308618/1‑149 (MQ=255)
                                                    ccGCGTGGGTGAAAAGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTgg                                                                                        >  3:65370/1‑149 (MQ=255)
                                                             tGAAAAGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCa                                                                               >  3:105481/1‑149 (MQ=255)
                                                              gAAAAGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTCATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCAt                                                                              >  6:120241/1‑149 (MQ=255)
                                                                aaaGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATgg                                                                            <  2:180608/149‑1 (MQ=255)
                                                                 aaGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGCCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGa                                                                           >  5:98305/1‑149 (MQ=255)
                                                                  aGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGAc                                                                          >  1:203607/1‑149 (MQ=255)
                                                                   gCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACg                                                                         >  6:181854/1‑149 (MQ=255)
                                                                    cTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCTGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGCGCCCGGCCATGGACGt                                                                        >  7:164760/1‑149 (MQ=255)
                                                                          aTCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGAttt                                                                   <  2:203607/148‑1 (MQ=255)
                                                                                 aaTCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGcc                                                           <  4:105481/149‑1 (MQ=255)
                                                                                 aaTCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGcc                                                           <  1:35061/149‑1 (MQ=255)
                                                                                   tCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCAt                                                         >  8:215588/1‑149 (MQ=255)
                                                                                      aCCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCTCGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCAt                                                                              <  7:326138/125‑1 (MQ=255)
                                                                                      aCCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCaa                                                                              >  8:326138/1‑124 (MQ=255)
                                                                                      aCCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAg                                                      >  3:35804/1‑149 (MQ=255)
                                                                                         gTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCcg                                                   >  4:307951/1‑149 (MQ=255)
                                                                                             gcgcCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGTCATGGACGTGATTTCCGCTGCCATCAGCCGCGcc                                               >  4:275172/1‑149 (MQ=255)
                                                                                              cgccgcCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGccc                                              >  5:54781/1‑149 (MQ=255)
                                                                                                ccgccCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCcgc                                            <  4:65370/149‑1 (MQ=255)
                                                                                                    ccTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAg                                         <  8:5314/148‑1 (MQ=255)
                                                                                                    ccTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAg                                         >  7:5314/1‑148 (MQ=255)
                                                                                                           cgccTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTCGAccc                                 <  3:307951/149‑1 (MQ=255)
                                                                                                               tGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACCTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTCGACCCGGc                              >  5:94091/1‑148 (MQ=255)
                                                                                                               tGCGCGAGCTGTTAGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACACCACGGGCCATGCACATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCAGCTGCCATCAGCCGCGCCCGCGAGTTCGACCCGGc                              >  6:304194/1‑148 (MQ=255)
                                                                                                                        tGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTc                                      >  1:176458/1‑131 (MQ=255)
                                                                                                                        tGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTc                                      <  2:176458/131‑1 (MQ=255)
                                                                                                                            cgacgaCCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTCGACCCGGCCAACAGCTGCAACg                >  6:165956/1‑149 (MQ=255)
                                                                                                                                 ggccGCTGCGGATCAGCAGCTGCCGGGCCATGGAACCCACCCCACGGGCCATCCAGATCAATGCAGAACTGGGCCCGGCCAAGGACGTGATATCCGCTGCCATCAGCCGCGCCCGCGAGATCGACCCGGCCAACAGCTGCAACGTGtta           >  2:137497/3‑148 (MQ=255)
                                                                                                                                   cggccgctgATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTCGACCCGGCCAACAGCTGCAACGTGTTCCg         >  7:284864/6‑149 (MQ=255)
                                                                                                                                   ccGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTCGACCCGGCCAACAGCTGCAACGTGTTCCg         >  4:35263/1‑149 (MQ=255)
                                                                                                                                    cGCTGCTGATGCGCATCGGCCGGGCCATGGAACCCTTCGCACGGGCCATGCAGATCTTTGCTGAACTGTGCCCGGCCATTGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTCGACCGGGCCATCAGCTTCAACGTGTTCCGc        <  8:164760/149‑1 (MQ=255)
                                                                                                                                    cGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTCGACCCGGCCAACAGCTGCAACGTGTTCCGc        <  5:181854/149‑1 (MQ=255)
                                                                                                                                    cGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTCGACCCGGCCAACAGCTGCAACGTGTTCCGc        >  8:324472/1‑149 (MQ=255)
                                                                                                                                         ctgATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTCGACCCGGCCAACAGCTGCAACGTGTTCCGCCTggg   <  7:215588/149‑1 (MQ=255)
                                                                                                                                          tgATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTCGACCCGGCCAACAGCTGCAACGTGTTCCGCCTGGGc  >  5:173105/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
ACATGAACCTGCTGGTGCTGTTCGAAGCCCTGATGATCGAACGCAACCTGACCCGCGTGGGTGAAAAGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTCGACCCGGCCAACAGCTGCAACGTGTTCCGCCTGGGC  >  NC_002947/4170307‑4170593

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