Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,170,450 C→A 8.5% R58S (CGC→AGC)  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,4500CA 8.5% 142.4 / 3.9 47R58S (CGC→AGC) PP_3669LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (27/16);  new base A (2/2);  total (29/18)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

ATGAACCTGCTGGTGCTGTTCGAAGCCCTGATGATCGAACGCAACCTGACCCGCGTGGGTGAAAAGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTCGACCCGGCCAACAGCTGCAACGTGTTCCGCCTGGGCCTT  >  NC_002947/4170309‑4170596
                                                                                                                                             |                                                                                                                                                  
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                                            >  7:5314/1‑148 (MQ=255)
                                                                                                  ccTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAg                                            <  8:5314/148‑1 (MQ=255)
                                                                                                         cgccTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTCGAccc                                    <  3:307951/149‑1 (MQ=255)
                                                                                                             tGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACCTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTCGACCCGGc                                 >  5:94091/1‑148 (MQ=255)
                                                                                                             tGCGCGAGCTGTTAGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACACCACGGGCCATGCACATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCAGCTGCCATCAGCCGCGCCCGCGAGTTCGACCCGGc                                 >  6:304194/1‑148 (MQ=255)
                                                                                                                      tGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTc                                         <  2:176458/131‑1 (MQ=255)
                                                                                                                      tGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTc                                         >  1:176458/1‑131 (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)
                                                                                                                                           tCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTCGACCCGGCCAACAGCTGCAACGTGTTCCGCCTGGGCCtt  >  2:22467/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
ATGAACCTGCTGGTGCTGTTCGAAGCCCTGATGATCGAACGCAACCTGACCCGCGTGGGTGAAAAGCTGTTTATCACCCAATCCACCGTCAGCGCCGCCCTCGCCCGCCTGCGCGAGCTGTTCGACGACCCGCTGCTGATCCGCAGCGGCCGGGCCATGGAACCCACCCCACGGGCCATGCAGATCTTTGCTGAACTGGGCCCGGCCATGGACGTGATTTCCGCTGCCATCAGCCGCGCCCGCGAGTTCGACCCGGCCAACAGCTGCAACGTGTTCCGCCTGGGCCTT  >  NC_002947/4170309‑4170596

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