| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,907,968 | Δ1 bp | 27.3% | coding (55/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_002947 | 5,907,968 | 0 | C | . | 27.3% | 45.8 / 17.9 | 22 | coding (55/1095 nt) | spuE | spermidine/putrescine ABC transporter substrate‑binding protein |
| Reads supporting (aligned to +/- strand): ref base C (11/5); new base . (3/3); total (14/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.24e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01 | |||||||||||
TTCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCTTCGATGCCGCAA > NC_002947/5907828‑5908116 | ttCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGGCATG‑CAGCGTcca > 1:439900/1‑148 (MQ=255) tCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTTCACCGTCGGCGCCGCTTGGACGCTGGCATG‑CAGCGTccaa < 2:439900/150‑4 (MQ=255) gTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCCCGAGCGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCCCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCaa > 1:245973/1‑150 (MQ=255) aaCACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGcc > 1:197435/1‑150 (MQ=255) ttGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGaa < 2:197435/150‑1 (MQ=255) ttGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATa < 1:175560/142‑1 (MQ=255) ttGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATa > 2:175560/1‑142 (MQ=255) gTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTgg > 2:464657/1‑150 (MQ=255) gCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGGCATG‑CAGCGTccaa < 2:86403/99‑4 (MQ=255) gCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGGCATG‑CAGCGTccaa > 1:86403/1‑96 (MQ=255) cgcgAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGGCATG‑CAGCGTccaa < 1:474768/85‑4 (MQ=255) cgcgAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGGCATG‑CAGCGTccaa > 2:474768/1‑82 (MQ=255) cgcgAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCaa < 2:202564/102‑1 (MQ=255) cgcgAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCaa > 1:202564/1‑102 (MQ=255) tAGTCGGACCAGTTATAGAAGTGCACCGTCGGTGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGc < 1:483719/88‑1 (MQ=255) tAGTCGGACCAGTTATAGAAGTGCACCGTCGGTGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGc > 2:483719/1‑88 (MQ=255) gTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAAAAATAGACACGTGGGGGTCCTCACAATCAGTGTCTTGGTAGCGACAATGCTGCCg > 2:155497/1‑150 (MQ=255) ccAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACg > 1:99567/1‑150 (MQ=255) ccAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACg > 1:275551/1‑150 (MQ=255) cgcTTGGACGCT‑GCATGCCGGCGTCATGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGc > 2:427260/1‑150 (MQ=255) gCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCTTCGAt < 1:464657/150‑1 (MQ=255) gCCAGCGTCAGGCCCGCGCCAGCCCTGGAGTCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCTTCGATGCCGCaa > 1:395431/1‑150 (MQ=255) | TTCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCTTCGATGCCGCAA > NC_002947/5907828‑5908116 |
| 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: ‑ |