Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,907,962:1 | +G | 13.8% | 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_002947 | 5,907,962 | 1 | . | G | 13.8% | 91.5 / 10.7 | 29 | coding (61/1095 nt) | spuE | spermidine/putrescine ABC transporter substrate‑binding protein |
Reads supporting (aligned to +/- strand): ref base . (15/10); new base G (2/2); total (17/12) | |||||||||||
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 |
CAGCAGCTTGCCTTCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCTTCGAT > NC_002947/5907816‑5908109 | cagcagCTTGCCTTCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCTGTCATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCa > 1:494125/1‑150 (MQ=255) ttCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAg > 2:114854/1‑150 (MQ=255) ccAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGc < 1:114854/150‑1 (MQ=255) ggTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCcgc < 2:279174/150‑1 (MQ=255) tCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTACGCGAGTGTGGTGGGGCATATGTAGTCGTACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCgcca > 1:383016/1‑150 (MQ=255) gAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGGCATG‑CAGCGTccaagcg < 1:432348/133‑7 (MQ=255) gAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGGCATG‑CAGCGTccaagcg > 2:432348/1‑127 (MQ=255) gAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGc > 1:405385/1‑150 (MQ=255) gAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGc > 2:480394/1‑150 (MQ=255) gacacATGTTGCACAGGCATGATCCCCGTGGCTTTCTTCAAGGCGGCGAGAGTGGTGGGGCCGGTGGAGCCGGAACAGTTCTAGAAGTGCACCGTCGGCGCCGCTTGGACGCG‑GCATGGCAGCGTCAAGCCCGCGCCAGCCCTCaatggc > 2:182626/2‑146 (MQ=255) acaTCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTg > 1:439318/1‑150 (MQ=255) ttGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATa < 2:383016/150‑1 (MQ=255) ttttCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGt < 2:290016/149‑1 (MQ=255) tttCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTc > 2:507883/1‑150 (MQ=255) aaGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACaa > 2:432113/1‑150 (MQ=255) aGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAAt < 1:410800/150‑1 (MQ=255) gtgtGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTg < 2:405385/150‑1 (MQ=255) tggtggGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGt > 1:311924/1‑150 (MQ=255) gCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTCGGTAGCGaga > 2:477993/1‑148 (MQ=255) cGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAAt < 2:311924/150‑1 (MQ=255) cGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAAt < 1:182626/150‑1 (MQ=255) cGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAAt > 1:240152/1‑150 (MQ=255) tGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAAtgct < 1:532199/150‑1 (MQ=255) aGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGt < 1:145722/112‑1 (MQ=255) aGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGt > 2:145722/1‑112 (MQ=255) ctgACGCTGGCATG‑CAGCGTCCGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGc < 2:42006/98‑1 (MQ=255) ctgACGCTGGCATG‑CAGCGTCCGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGc > 1:42006/3‑100 (MQ=255) gCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACGGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCTTCGAt > 1:478010/1‑150 (MQ=255) gCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCTTCGAt > 2:520205/1‑150 (MQ=255) | CAGCAGCTTGCCTTCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCT‑GCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCTTCGAT > NC_002947/5907816‑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: ‑ |