Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,896,422 C→T 10.8% intergenic (+10/‑198) 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,4220CT10.8% 94.2 / 3.7 37intergenic (+10/‑198)PP_4304/sbp‑IVIC family cation transporter/sulfate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (17/16);  new base T (2/2);  total (19/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 = 7.25e-01

TCTTTACCGCCGAGCTGGCCAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTA  >  NC_002947/4896276‑4896571
                                                                                                                                                  |                                                                                                                                                     
tCTTTACCGCCGAGCTGGCCAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCC                                                                                                                                                    >  1:45730‑M1/1‑150 (MQ=255)
             tctGGCCAATGCCATGCGTGGGGAGCGGCTGCAGCATGATTGCACGACCTGCCACAAGAAGACCCACGCGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCGCAAAGCGTGATAGCGCCTGCGCCGACCTCGGCGCGG                                                                                                                                       <  2:67073‑M1/149‑1 (MQ=255)
                 gCCAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAA                                                                                                                                   >  1:299051‑M1/1‑150 (MQ=255)
                              gTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACa                                                                                                                      >  1:456969/1‑150 (MQ=255)
                                           gcagcaTGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTg                                                                                                         >  2:135693/1‑150 (MQ=255)
                                                      gCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTccccc                                                                                              <  2:45730/150‑1 (MQ=255)
                                                              tGCCACAAGAAGACCCACGCGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGc                                                                                      >  1:144782/1‑150 (MQ=255)
                                                              tGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGc                                                                                      >  1:353246/1‑150 (MQ=255)
                                                                ccACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCgg                                                                                    <  2:166565/150‑1 (MQ=255)
                                                                 cacaAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAtccgctcctacgggggaccgcgcacggtccgcgcctgtaggagcgga                                                                                   <  2:299051‑M1/150‑48 (MQ=255)
                                                                         gaCCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGtttacccg                                                                           <  2:173499‑M2/150‑9 (MQ=255)
                                                                            ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTc                                                                                                                           <  2:427080/99‑1 (MQ=255)
                                                                            ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTc                                                                                                                           >  1:427080/1‑99 (MQ=255)
                                                                            ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACAGTGCGCGGTCCCCCGTAGGAGCgg                                                                                    >  2:338770/1‑138 (MQ=255)
                                                                            ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACAGTGCGCGGTCCCCCGTAGGAGCgg                                                                                    <  1:338770/138‑1 (MQ=255)
                                                                            ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAA                                                                                                                                   <  1:223216‑M1/91‑1 (MQ=255)
                                                                            ccACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAA                                                                                                                                   >  2:223216‑M1/1‑91 (MQ=255)
                                                                             cACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGCGCACGGTccg                                                                                                <  2:145466/125‑2 (MQ=255)
                                                                             cACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGCGCACGGTccg                                                                                                >  1:145466/1‑124 (MQ=255)
                                                                              aCGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAcccgctcctacgggggaccgcgcacggtccgcgcctgtaggagcggatttacccgcgaa                                                                       <  2:366562‑M1/149‑60 (MQ=255)
                                                                              aCGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAcccgctcctacgggggaccgcgcacggtccgcgcctgtaggagcggatttacccgcgaa                                                                       >  1:366562‑M1/1‑90 (MQ=255)
                                                                                     gggggCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGtttacccgcgaagaggccgg                                                               >  1:135693‑M2/1‑130 (MQ=255)
                                                                                      ggggCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGTGCCGGCCTCTTCGCGGGTAAAcccgctcctacgggggaccgcgcacggtccgcgcctgtaggagcggatttacccgcgaagaggccggc                                                              >  2:452999‑M1/1‑82 (MQ=255)
                                                                                              gttctgttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAg                                                      >  1:158673/1‑150 (MQ=255)
                                                                                                   gttctCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGCCGCGGACCGGGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCa                                                 <  2:250841/150‑1 (MQ=255)
                                                                                                           tgtgGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGTGCCGGCCTCTTCGCGGGTAAAcccgctcctacgggggaccgcgcacggtccgcgcctgtaggagcggatttacccgcgaagaggccggcgcaggcgctatcacgctttgg                                         <  1:452999‑M1/150‑90 (MQ=255)
                                                                                                           tgtgGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAcccgctcctacgggggaccgcgcacggtccgcgcctgtaggagcgga                                                                                   <  1:408324‑M1/108‑48 (MQ=255)
                                                                                                           tgtgGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAcccgctcctacgggggaccgcgcacggtccgcgcctgtaggagcgga                                                                                   >  2:408324‑M1/1‑61 (MQ=255)
                                                                                                                gAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAcccgctcctacgggggaccgcgcacggtccgcgcctgtaggagcggatttacc                                                                             >  2:124205‑M1/1‑56 (MQ=255)
                                                                                                                gAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAAcccgctcctacgggggaccgcgcacggtccgcgcctgtaggagcggatttacc                                                                             <  1:124205‑M1/109‑54 (MQ=255)
                                                                                                                              ccAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGtttacccgcgaagaggccggcacaggcgctatcacgctttgggaaacaacccattcccaca                      <  1:452999‑M2/150‑62 (MQ=255)
                                                                                                                               cAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAcc                     >  1:269647/1‑150 (MQ=255)
                                                                                                                                aaaGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAAccc                    <  2:265573/150‑1 (MQ=255)
                                                                                                                                     gTGATAGCGCCTGTGCCGGCCTCTTCGCGGGTAAAcccgctcctacgggggaccgcgcacggtccgcgcctgtaggagcggat                                                                                  >  2:175848‑M1/1‑35 (MQ=255)
                                                                                                                                     gTGATAGCGCCTGTGCCGGCCTCTTCGCGGGTAAAcccgctcctacgggggaccgcgcacggtccgcgcctgtaggagcggat                                                                                  <  1:175848‑M1/83‑49 (MQ=255)
                                                                                                                                        aTAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAac            <  1:173499/150‑1 (MQ=255)
                                                                                                                                              ccTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAacccaccc      >  2:501834/1‑150 (MQ=255)
                                                                                                                                                  cgcCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTta  <  1:501834/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TCTTTACCGCCGAGCTGGCCAATGCCATGCGTGGGGAGCAGCTGCAGCATGATTGCCCGACCTGCCACAAGAAGACCCACGAGCAGGGGGCGGCGTTCTGTTCTCGCTGTGGGAATGGGTTGTTTCCCAAAGCGTGATAGCGCCTGCGCCGGCCTCTTCGCGGGTAAATCCGCTCCTACAGGCGCGGACCGTGCGCGGTCCCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAGAAGCATAACCAAAGCATTTTTTGTTCTTTAACCCGTATAAACCCACCCGTTA  >  NC_002947/4896276‑4896571

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