Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,896,399 T→C 10.0% intergenic (+32/‑176) PP_4304 → / → sbp‑I VIC family cation transporter/sulfate ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,896,3990TC10.0% 101.0 / 3.0 42intergenic (+32/‑176)PP_4304/sbp‑IVIC family cation transporter/sulfate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (18/18);  new base C (2/2);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

CCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCC  >  NC_002947/4896255‑4896546
                                                                                                                                                 |                                                                                                                                                       
ccATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGc                                                                                                                                                     >  7:186547/1‑149 (MQ=255)
    gCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGGTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTcc                                                                                                                                                  >  8:253164/1‑148 (MQ=255)
             acgttgcgtccttgttcttgaagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgcGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGcgcg                                                                                                                                        <  3:176819‑M1/39‑1 (MQ=255)
                    gtccttgttcttgaagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTg                                                                                                                                 <  7:117620‑M1/46‑1 (MQ=255)
                       caTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCTCTGTGGGAATGTGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTAGCGGGTAATTACCCTACAACAGGCGCTGGCCGTgcgc                                                                                                                              >  5:228533/1‑149 (MQ=255)
                       caTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTgcgc                                                                                                                              >  7:108197/1‑149 (MQ=255)
                        aTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGcgcg                                                                                                                             <  8:108197/149‑1 (MQ=255)
                              gCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTcc‑‑‑c                                                                                                                     <  8:264064/148‑1 (MQ=255)
                                 aagggtgttggggcagattgccagcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCg                                                                                                                  >  8:117620‑M1/91‑148 (MQ=255)
                                    ccTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAgg                                                                                                              <  4:145464/149‑1 (MQ=255)
                                         cacaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGCCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCgg                                                                                                         <  7:253164/149‑1 (MQ=255)
                                             agaagaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttt                                                                                                     >  8:186547‑M2/1‑146 (MQ=255)
                                                    ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGAtttacccgc                                                                                               <  5:1442‑M2/148‑10 (MQ=255)
                                                     cACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGGCCAATTCCCGGGTAAAACCCGACAAACCGGCCCGGACCGTGCCCGGACC‑‑‑CCCGTAGGAGCGGGTTTAccccc                                                                                               >  2:223560/1‑146 (MQ=255)
                                                       cGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGa                                                                                                                                      >  8:96086/1‑109 (MQ=255)
                                                       cGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGa                                                                                                                                      <  7:96086/109‑1 (MQ=255)
                                                       cGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCggg                                                                                                        <  8:174384/136‑1 (MQ=255)
                                                       cGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCggg                                                                                                        >  7:174384/1‑136 (MQ=255)
                                                       cGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGG‑‑‑GGACCGCGCACGGTCCGCGCCTGTAGGAGCGGAtttacccgcgaag                                                                                           >  6:1442‑M2/1‑136 (MQ=255)
                                                        gcgactataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgcGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAaga                                                                                          >  4:176819‑M1/68‑149 (MQ=255)
                                                            ctataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTAc                                                                                                   <  8:169467‑M1/73‑1 (MQ=255)
                                                            ctataacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTAc                                                                                                   >  7:169467‑M1/64‑136 (MQ=255)
                                                             gggggCGGCGTTCTGTTCGCGCTGGGGCAAGGGGTTGTTCCCCAAAGCGTGATAGCGCCTGCGCGGGCCTCTTGGCGGGTAAATCCGCTCCTACCGGCGGGGACGGTGCGCGTTCC‑‑‑CCCTTTGGAGCGGGTTTACCCGCGAAGAGGCCg                                                                                     <  2:308491/149‑1 (MQ=37)
                                                               taacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGaa                                                                                            >  8:31902‑M1/61‑140 (MQ=255)
                                                               taacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGaa                                                                                            <  7:31902‑M1/80‑1 (MQ=255)
                                                               taacgggtgggtttatacgggttaaagaacaaaaaatgctttggttatgcttctgcctgtGCCGCCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCg                                                                                     >  2:99861‑M1/61‑147 (MQ=255)
                                                                       ttctgttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCACTTAGCGGGTAAATCCGCTCCTAAAGGCGCGGACCGGGCGCGGTGC‑‑‑CCCGTAGGAGCGGGTTTAACCGCGAAGAGACCGGCACAgg                                                                              >  1:269850/1‑146 (MQ=255)
                                                                         ctgttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttt                                                                                                     >  8:35444‑M2/1‑118 (MQ=255)
                                                                         ctgttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttt                                                                                                     <  7:35444‑M2/121‑4 (MQ=255)
                                                                           gttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACTGGCGCGGACCGTGCGCGGGCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggc                                                                                   >  6:251500‑M2/1‑116 (MQ=255)
                                                                           gttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggc                                                                                   <  5:251500‑M2/137‑22 (MQ=255)
                                                                                         aaTGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCACATAGCCGGTAAATACGCTACAAAACGCCCGGACCCTGCCCGGTCC‑‑‑CCCCTACGAGCGGGTTTAACCCCCAAGAGGCCGGCACAAGCAAAAGCAAAACCAAAACaat                                                        >  6:297037/1‑147 (MQ=255)
                                                                                                    tcgcAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCtt                                              <  5:297037/145‑1 (MQ=255)
                                                                                                      ccAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttta                                                                                                    <  5:305389‑M2/93‑5 (MQ=255)
                                                                                                      ccAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGttta                                                                                                    >  6:305389‑M2/1‑89 (MQ=255)
                                                                                                       cAAAGCGTGATAGCGCCTGCGCCGGCCTAATAGCGGGTAAAAACCCACCAACAGGGGCGGGCCCGGGGCGGTGC‑‑‑CCCGGATGGGCGGGGTTTACAGCCAAAAGGCCCGCACAAAGAGAAGCAAAAACCAAAAAATtttttttttata                                             >  8:14087/1‑143 (MQ=255)
                                                                                                       cAAAGCGGTATAGCGCCTGAGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGCCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGTAGCGGGTTTACCCGCTAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTc                                                <  7:14087/143‑1 (MQ=38)
                                                                                                        aaaGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAcc                                         >  7:164824/1‑149 (MQ=255)
                                                                                                                          cgcCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGt                        <  1:223560/148‑1 (MQ=255)
                                                                                                                                GCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAGCCCGTAATAACCCACGCGTTATAGt                  <  5:267789‑M1/149‑1 (MQ=255)
                                                                                                                                            gTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtctt       <  1:319028‑M2/149‑98 (MQ=255)
                                                                                                                                                ATCCGCTCCTACAGGCGCTGACCGTTCGCGGTCC‑‑‑CCCGTAGTAGCGGGTTTACCCGCGAAGAGGCCGTCACAGGCAGAAGCATAACCAAACCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGccc  >  6:228533‑M1/1‑149 (MQ=255)
                                                                                                                                                 tCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGAtttacccgcgaagaggccggcgcaggcgctatcacgctttgggaaacaacccattcccacagcgagaacagaacgccgccccctgctcgtgggtcttcttgt  <  2:168079‑M2/149‑103 (MQ=255)
                                                                                                                                                 |                                                                                                                                                       
CCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCC‑‑‑CCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTATAGTCGCTGGCAATCTGCCC  >  NC_002947/4896255‑4896546

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