Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,458,070 C→G 11.8% intergenic (+142/‑63) PP_2150 → / → sthA thiamine biosynthesis lipoprotein/pyridine nucleotide transhydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,458,0700CG11.8% 130.3 / 4.6 48intergenic (+142/‑63)PP_2150/sthAthiamine biosynthesis lipoprotein/pyridine nucleotide transhydrogenase
Reads supporting (aligned to +/- strand):  ref base C (22/20);  new base G (2/4);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGTGAAAGCACAAGCAGGGATGATTGCCGCCAGGGAATGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCA  >  NC_002947/2457922‑2458213
                                                                                                                                                     |                                                                                                                                               
ccAGTGAAAGCACAAGCAGGGATGATTGCCGCCAGGGAATGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGccc                                                                                                                                                 >  4:94581/1‑149 (MQ=255)
   gTGAAAGCACAAGCAGGGATGATTGCCGCCAGGGAATGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCcgg                                                                                                                                              >  4:59052/1‑149 (MQ=255)
    tGAAAGCACAAGCAGGGATGATTGCCGCCAGGGAATGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCcggc                                                                                                                                             <  1:57554/149‑1 (MQ=255)
    tGAAAGCACAAGCAGGGATGATTGCCGCCAGGGAATGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCcggc                                                                                                                                             >  1:77421/1‑149 (MQ=255)
         gCACAAGCAGGGATGATTGCCGCCAGGGAATGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGCCGCCGGGGCGCCAt                                                                                                                                        <  2:254389/149‑1 (MQ=255)
           acaAGCAGGGATGATTGCCGCCAGGGAATGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGa                                                                                                                                      >  7:79092/1‑149 (MQ=255)
            caAGCAGGGATGATTGCCGCCAGGGAATGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGCCGCCGGGGCGCCATGag                                                                                                                                     <  8:190147/149‑2 (MQ=255)
             aaGCAGGGATGATTGCCGCCAGGGAATGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATc                                                                                                                                    >  8:138300/1‑149 (MQ=255)
                                cAGGGAATGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGt                                                                                                                 <  3:59052/149‑1 (MQ=255)
                                     aaTGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCtt                                                                                                            <  8:79092/149‑1 (MQ=255)
                                           tGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGc                                                                                                      <  3:94581/149‑1 (MQ=255)
                                                     gTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCGCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTg                                                                                              >  5:261853/1‑147 (MQ=255)
                                                      tAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCcggcg                                                                                                                                            <  8:61477/100‑1 (MQ=255)
                                                      tAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCcggcg                                                                                                                                            >  7:61477/1‑100 (MQ=255)
                                                       aGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCgcga                                                                                                                         <  8:186689/118‑1 (MQ=255)
                                                        gTGCGGGCAAAATTGGCCTACGACGAGACAAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGa                                                                                         >  7:216931/1‑149 (MQ=255)
                                                              gCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTc                                                                                                >  5:228445/1‑136 (MQ=255)
                                                              gCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTc                                                                                                <  6:228445/136‑1 (MQ=255)
                                                              gCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGCCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCGTGCGGGCCAGTTCTGAAGGAGTACg                                                                                    >  3:75576/1‑148 (MQ=255)
                                                               cAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCa                                                                                  >  7:15147/1‑149 (MQ=255)
                                                                aaaaTTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCAt                                                                                 <  5:122871/149‑1 (MQ=255)
                                                                    ttGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGt                                                                                       >  3:138809/1‑139 (MQ=255)
                                                                    ttGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGt                                                                                       <  4:138809/139‑1 (MQ=255)
                                                                       gcagacgacgAGACCAAGTGTTAATGTGCGCGGAGTTCAGGGTTGATTAGTCTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGt                                                                           <  8:216931/144‑1 (MQ=255)
                                                                       gCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTc                                                                          <  2:77421/149‑1 (MQ=255)
                                                                         cTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTGCTGAAGGAGTACGCATGGCTGTCTa                                                                        >  3:101471/1‑149 (MQ=255)
                                                                          tACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTa                                                                        >  1:238094/1‑149 (MQ=255)
                                                                                         ggTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAgtg                                                        >  2:100761/1‑149 (MQ=255)
                                                                                          gTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAgtgg                                                       >  4:276811/1‑149 (MQ=255)
                                                                                          gTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCATGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAgtgg                                                       >  5:209910/1‑149 (MQ=255)
                                                                                           ttAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGCCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAgtggt                                                      <  4:75576/149‑1 (MQ=255)
                                                                                                             aGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTccggcccggc                                    >  3:176442/1‑149 (MQ=255)
                                                                                                                   gaatAGACTGCAACC‑GAATTTTCTCATGGCGCCCCGGCTGCATGATCGAGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGTATGAGTACGCATGGCTGTCTACAACTACGACGTAGTTGTTCTGGGTTCCGGCCCTGCCGGTg                               <  6:295222/145‑1 (MQ=255)
                                                                                                                    aTTAGACTGCACCC‑GAATTTTCTCATGGCGCCGCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGaa                             <  6:261853/149‑1 (MQ=255)
                                                                                                                     ttAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAg                            <  1:100761/149‑1 (MQ=255)
                                                                                                                     ttAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAg                            <  8:15147/149‑1 (MQ=255)
                                                                                                                           tGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCgg                      >  2:217978/1‑149 (MQ=255)
                                                                                                                                   gAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAgtggtg                                                     >  8:267052/1‑111 (MQ=255)
                                                                                                                                   gAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAgtggtg                                                     <  7:267052/111‑1 (MQ=255)
                                                                                                                                   gAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGaa               <  1:217978/149‑1 (MQ=255)
                                                                                                                                         tctcATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGaa                             <  1:98037/129‑1 (MQ=255)
                                                                                                                                         tctcATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGaa                             >  2:98037/1‑129 (MQ=255)
                                                                                                                                          atcATGCCGCCGGGGCGCCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTa                                                           <  5:135031/97‑1 (MQ=255)
                                                                                                                                          atcATGCCGCCGGGGCGCCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCAAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTa                                                           >  6:135031/2‑98 (MQ=255)
                                                                                                                                           tcATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACgccgcc        >  5:107656/1‑148 (MQ=255)
                                                                                                                                             aTGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCa       >  1:153493/1‑147 (MQ=255)
                                                                                                                                              tGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAg                            <  6:104868/125‑1 (MQ=255)
                                                                                                                                               ggCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAg    >  2:49808/1‑148 (MQ=255)
                                                                                                                                                gcgcCCCGGCGGCATGATCATGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCa  <  6:209910/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CCAGTGAAAGCACAAGCAGGGATGATTGCCGCCAGGGAATGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCC‑GAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCA  >  NC_002947/2457922‑2458213

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