Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,907,962:1 +G 28.6% coding (61/1095 nt) spuE ← spermidine/putrescine ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,907,9621.G28.6% 32.5 / 22.7 21coding (61/1095 nt)spuEspermidine/putrescine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (10/5);  new base G (3/3);  total (13/8)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCAGCAGCTTGCCTTCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCTTCGAT  >  NC_002947/5907813‑5908109
                                                                                                                                                      |                                                                                                                                                   
ggCCAGCAGCTTGCCTTCCAGGGTTTCGTTGGAGTCGAAGACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCt                                                                                                                                                      <  1:387080/150‑1 (MQ=255)
               ttCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGGCATGCAGCGTcca                                                                                                                                      >  1:439900/1‑148 (MQ=255)
                tCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTTCACCGTCGGCGCCGCTTGGACGCTGGCATGCAGCGTccaa                                                                                                                                     <  2:439900/150‑4 (MQ=255)
                                 gTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCCCGAGCGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCCCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCaa                                                                                                                    >  1:245973/1‑150 (MQ=255)
                                     aaCACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGcc                                                                                                                >  1:197435/1‑150 (MQ=255)
                                             ttGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGaa                                                                                                        <  2:197435/150‑1 (MQ=255)
                                                       ttGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATa                                                                                                      <  1:175560/142‑1 (MQ=255)
                                                       ttGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATa                                                                                                      >  2:175560/1‑142 (MQ=255)
                                                                gTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTgg                                                                                     >  2:464657/1‑150 (MQ=255)
                                                                   gCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGGCATGCAGCGTccaa                                                                                                                                     <  2:86403/99‑4 (MQ=255)
                                                                   gCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGGCATGCAGCGTccaa                                                                                                                                     >  1:86403/1‑96 (MQ=255)
                                                                    cTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCt                                                                                                                                                      <  2:246475/82‑1 (MQ=255)
                                                                    cTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCt                                                                                                                                                      >  1:246475/1‑82 (MQ=255)
                                                                                 cgcgAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGGCATGCAGCGTccaa                                                                                                                                     <  1:474768/85‑4 (MQ=255)
                                                                                 cgcgAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGGCATGCAGCGTccaa                                                                                                                                     >  2:474768/1‑82 (MQ=255)
                                                                                 cgcgAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCaa                                                                                                                    <  2:202564/102‑1 (MQ=255)
                                                                                 cgcgAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCaa                                                                                                                    >  1:202564/1‑102 (MQ=255)
                                                                                                       tAGTCGGACCAGTTATAGAAGTGCACCGTCGGTGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGc                                                                                                            <  1:483719/88‑1 (MQ=255)
                                                                                                       tAGTCGGACCAGTTATAGAAGTGCACCGTCGGTGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGc                                                                                                            >  2:483719/1‑88 (MQ=255)
                                                                                                         gTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAAAAATAGACACGTGGGGGTCCTCACAATCAGTGTCTTGGTAGCGACAATGCTGCCg                                            >  2:155497/1‑150 (MQ=255)
                                                                                                               ccAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACg                                      >  1:99567/1‑150 (MQ=255)
                                                                                                               ccAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACg                                      >  1:275551/1‑150 (MQ=255)
                                                                                                                                          cgcTTGGACGCTGCATGCCGGCGTCATGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGc           >  2:427260/1‑150 (MQ=255)
                                                                                                                                                   gCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCTTCGAt  <  1:464657/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GGCCAGCAGCTTGCCTTCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCTTCGAT  >  NC_002947/5907813‑5908109

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: