Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,483,885 C→G 15.5% intergenic (+23/‑45) 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,8850CG15.5% 68.8 / 6.0 31intergenic (+23/‑45)yhdW/yhdXamino acid ABC transporter‑binding protein YhdW/amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (14/12);  new base G (3/2);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

GCAAGTGGGTAACTATGGCGAAGTGTTCGAGAAGAACCTGGGCCAGAGCACTGACTTGAAGATCGACCGTGGCATGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCCGGGGCACTTCATGCAAAATCAAATCGGCGCACCCAAGGGCTTGTCCCTTAACGATCCGCGTGTGCGCGCGTGGCTGTTCCAGATCCTCACGATCGTCTTCGTGGTCGGCCTGG  >  NC_002947/1483736‑1484032
                                                                                                                                                     |                                                                                                                                                    
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                                                                                                                          <  1:130389/150‑1 (MQ=255)
                            gagaAGAACCTGGGCCAGAGCACTGACTTGAAGATCGACCGTGGCATGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTa                                                                                                                          >  2:465758/1‑150 (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                                                                                                                                             >  2:649146/1‑84 (MQ=255)
                                                                         aTGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGGATGccgccgcag                                                                                                                                             <  1:649146/86‑3 (MQ=255)
                                                                                    gtgGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTa                                                                                                                          >  2:695431/1‑94 (MQ=255)
                                                                                    gtgGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTa                                                                                                                          <  1:695431/94‑1 (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)
                                                                                                                                                   aTCCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCCGGGGCACTTCATGCAAAATCAAATCGGCGCACCCAAGGGCTTGTCCCTTAACGATCCGCGTGTGCGCGCGTGGCTGTTCCAGATCCTCACGATCGTCTTCGTGGTCGGCCTgg  >  1:588999/1‑150 (MQ=255)
                                                                                                                                                    tCCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCCGGGGCACTTCATGCAAAATCAAATCGGCGCACCCAAGGGCTTGTCCCTTAAc                                                               >  2:200529/1‑88 (MQ=255)
                                                                                                                                                    tCCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCCGGGGCACTTCATGCAAAATCAAATCGGCGCACCCAAGGGCTTGTCCCTTAAc                                                               <  1:200529/88‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCAAGTGGGTAACTATGGCGAAGTGTTCGAGAAGAACCTGGGCCAGAGCACTGACTTGAAGATCGACCGTGGCATGAACGCCCTGTGGAACAACGGCGGCATCCAGTACGCGCCTCCAGTGCGCTGATGGCTGCACCCTGCGGCGGCATCCCGCCGCCGCAGGCTGTTCGAATCCCTACTGTCCGGGGCACTTCATGCAAAATCAAATCGGCGCACCCAAGGGCTTGTCCCTTAACGATCCGCGTGTGCGCGCGTGGCTGTTCCAGATCCTCACGATCGTCTTCGTGGTCGGCCTGG  >  NC_002947/1483736‑1484032

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: