Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,483,882 C→G 13.8% intergenic (+20/‑48) yhdW → / → yhdX amino acid ABC transporter‑binding protein YhdW/amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,483,8820CG13.8% 69.9 / 5.5 29intergenic (+20/‑48)yhdW/yhdXamino acid ABC transporter‑binding protein YhdW/amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (14/11);  new base G (2/2);  total (16/13)
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.38e-01

AAGCAAGTGGGTAACTATGGCGAAGTGTTCGAGAAGAACCTGGGCCAGAGCACTGACTTGAAGATCGACCGTGGCATGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCCGGGGCACTTCATGCAAAATCAAATCGGCGCACCCAAGGGCTTGTCCCTTAACGATCCGCGTGTGCGCGCGTGGCTGTTCCAGATCCTCACGATCGTCTTCGTGGTCG  >  NC_002947/1483734‑1484026
                                                                                                                                                    |                                                                                                                                                 
aagcaagTGGGTAACTATGGCGAAGTGTTCGAGAAGAACCTGGGCCAGAGCACTGACTTGAAGATCGACCGTGGCATGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCa                                                                                                                                                  >  2:739112/1‑150 (MQ=255)
  gcaagTGGGTAACTATGGCGAAGTGTTCGAGAAGAACCTGGGCCAGAGCACTGACTTGAAGATCGACCGTGGCATGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATc                                                                                                                                                <  2:588430/150‑1 (MQ=255)
           tAACTATGGCGAAGTGTTCGAGAAGAACCTGGGCCAGAGCACTGACTTGAAGATCGACCGTGGCATGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCccgccgccg                                                                                                                                       >  2:130389/1‑150 (MQ=255)
             aCTATGGCGAAGTGTTCGAGAAGAACCTGGGCCAGAGCACTGACTTGAAGATCGACCGTGGCATGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCa                                                                                                                                     <  2:424773/150‑1 (MQ=255)
                  ggCGAAGTGTTCGAGAAGAACCTGGGCCAGAGCACTGACTTGAAGATCGACCGTGGCATGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTg                                                                                                                                <  1:739112/150‑1 (MQ=255)
                           ttCGAGAAGAACCTGGGCCAGAGCACTGACTTGAAGATCGACCGTGGCATGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGGATGccgccgcagg                                                                                                                                      <  1:638321/135‑4 (MQ=255)
                           gtCGAGAAGAACCTGGGCCAGAGCACTGACTTGAAGATCGACCGTGGCATGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGGATGccgccgcagg                                                                                                                                      >  2:638321/2‑132 (MQ=255)
                            tCGAGAAGAACCTGGGCCAGAGCACTGACTTGAAGATCGACCGTGGCATGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATccc                                                                                                                      <  1:148524/150‑1 (MQ=255)
                              gagaAGAACCTGGGCCAGAGCACTGACTTGAAGATCGACCGTGGCATGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTa                                                                                                                    >  2:465758/1‑150 (MQ=255)
                              gagaAGAACCTGGGCCAGAGCACTGACTTGAAGATCGACCGTGGCATGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTa                                                                                                                    <  1:130389/150‑1 (MQ=255)
                                       gtGGGCCAGAGCACTGACTTGAAGATCGACCGTGGCATGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCCggg                                                                                                           >  2:250772/2‑150 (MQ=255)
                                            ccAGAGCACTGACTTGAAGATCGACCGTGGCATGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCTCGCCGCCGCAGGCTGTTTGAATCCCTACTGTCCGGGGCACt                                                                                                     >  1:394818/1‑150 (MQ=255)
                                             cAGAGCACTGACTTGAAGATCGACCGTGGCATGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCCGGGGCACtt                                                                                                    >  2:540171/1‑150 (MQ=255)
                                                         ttGAAGATCGACCGTGGCATGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGAGGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCGGGGCACCTTCATCCCaagtca                                                                                        >  2:150871/1‑150 (MQ=255)
                                                                           aTGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGGATGccgccgcag                                                                                                                                       <  1:649146/86‑3 (MQ=255)
                                                                           aTGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGGATGccgccgcag                                                                                                                                       >  2:649146/1‑84 (MQ=255)
                                                                                      gtgGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTa                                                                                                                    <  1:695431/94‑1 (MQ=255)
                                                                                      gtgGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTa                                                                                                                    >  2:695431/1‑94 (MQ=255)
                                                                                           acaacGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCCGGGGCACTTCATGCAAAATCAAATCGGCGCACCCAAGGGCTTGTCCCTTAACGAt                                                      >  2:549801/1‑150 (MQ=255)
                                                                                              acGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATGCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCCGGGGCACTTCATGCAAAATCAAATCGGCGCACCCAAGGGCTTGTCCCTTAACGATCcg                                                   >  1:507586/1‑150 (MQ=255)
                                                                                               gggcggcATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCCGGGGCACTTCATGCAAAATCAAATCGGCGCACCCAAGGGCTTGTCCCTTAACGATCcgc                                                  <  1:264552/149‑1 (MQ=255)
                                                                                                             aCGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCCGGGGCACTTCATGCAAAATCAAATCGGCGCACCCAACGGCTTGTCCCTTAACGATCCGCGTGGGCGCGCGTgg                                    <  2:521968/150‑1 (MQ=255)
                                                                                                               gcgcCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCCGGGGCACTTCATGCAAAATCAAATCAGCGCGCCCAAGGGCTTGTCCCTTAACGATCCGCGTGTGCGCGCGTGGCt                                  >  1:406352/1‑150 (MQ=255)
                                                                                                                 gcCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCCGGGGCACTTCATGCAAAATCAAATCGGCGCACCCAAGGGCTTGTCCCTTAACGATCCGCGTGTGCGCGCGTGGCTGt                                >  2:558694/1‑150 (MQ=255)
                                                                                                                                aTGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCCGGGGCACTTCATGCAAAATCAAATCGGCGCACCCAAGGGCTTGTCCCTTAACGATCCGCGTGTGCGCGCGTGGCTGTTCCAGATCCTCACGa                 <  1:549801/150‑1 (MQ=255)
                                                                                                                                     tGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCCGGGGCACTTCATGCAAAATCAAATCGGCGCACCCAAGGGCTTGTCCCTTAACGATCCGCGTGTGCGCGCGTGGCTGTTCCAGAACCTCACGATc               >  2:393715/1‑147 (MQ=255)
                                                                                                                                     tGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCCGGGGCACTTCATGCAAAATCAAATCGGCGCACCCAAGGGCTTGTCCCTTAACGATCCGCGTGTGCGCGCGTGGCTGTTCCAGAACCTCACGATc               <  1:393715/147‑1 (MQ=255)
                                                                                                                                          ccTGAGGCGGCATCCCGCCGCCGCGGGATGTGCGAATCCCTACTGGCCGGGGCACTTCATGCAAAATCAAATCGGCGCCCCCCAGGGCTTGTGCCTTCACGCCGCTCGCATGTCCGCGGTTCTGTTCCAGATCCTCACGATCGTCTTCGt       >  2:612000/1‑150 (MQ=255)
                                                                                                                                           cTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCCGGGGCACTTCATGCAAAATCAAATCGGCGCACCCAAGGGCTTGTCCCTTAACGATCCGCGTGTGCGCGCGTGGCTGTTCCAGATCCTCACGATCGTCTTCGTg      <  1:250772/150‑1 (MQ=255)
                                                                                                                                               ggcggcATCCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCCGGGGCACTTCATGCAAAATCAAATCGGCGCACCCAAGGGCTTGTCCCTTAACGATCCGCGTGTGCGCGCGTGGCTGTTCCAGATCCTCACGATCATCTTCGTGGTCg  <  1:150871/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AAGCAAGTGGGTAACTATGGCGAAGTGTTCGAGAAGAACCTGGGCCAGAGCACTGACTTGAAGATCGACCGTGGCATGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCCGGGGCACTTCATGCAAAATCAAATCGGCGCACCCAAGGGCTTGTCCCTTAACGATCCGCGTGTGCGCGCGTGGCTGTTCCAGATCCTCACGATCGTCTTCGTGGTCG  >  NC_002947/1483734‑1484026

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