| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,896,266 | C→T | 12.9% | intergenic (‑35/+150) | cysU ← / ← sbp‑II | sulfate ABC transporter permease CysT/sulfate ABC transporter substrate‑binding protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 5,896,266 | 0 | C | T | 12.9% | 72.4 / 3.9 | 31 | intergenic (‑35/+150) | cysU/sbp‑II | sulfate ABC transporter permease CysT/sulfate ABC transporter substrate‑binding protein |
| Reads supporting (aligned to +/- strand): ref base C (19/8); new base T (2/2); total (21/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.77e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01 | |||||||||||
CCAGCGGTATCAGCACGATCAGGCTGAGGTACACCAAGGTGTAGCCCAGCGTCAGCCCGAAGCCGGGTATGACGGGGGAGATACGACGCGACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCGCTTTCCT > NC_002947/5896138‑5896408 | ccAGCGGTATCAGCACGATCAGGCTGAGGTACACCAAGGTGTAGCCCAGCGTCAGCCCGAAGCCGGGTATGACGGGGGAGATACGACGCGACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCt > 2:428391/1‑150 (MQ=255) aTCAGGCTGAGGTACACCAAGGTGTAGCCCAGCGTCAGCCCGAAGCCGGGTATGACGGGGGAGATACGACGCGACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGgctg > 1:720034/1‑150 (MQ=255) tgggctACACCTTGGTGTAGCCCAGCGTCAGCCCGAAGCCGGGTATGACGGGGGAGATACGACGCGACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTgggg < 1:452878/135‑1 (MQ=255) gTACACCAAGGTGTAGCCCAGCGTCAGCCCGAAGCCGGGTATGACGGGGGAGATACGACGCGACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGccc > 1:636309/1‑150 (MQ=255) cacCAAGGTGTAGCCCAGCGTCAGCCCGAAGCCGGGTATGACGGGGGAGATACGACGCGACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCAtt > 1:675507/1‑150 (MQ=255) acCAAGGTGTAGCCCAGCGTCAGCCCGAAGCCGGGTATGACGGGGGAGATACGACGCGACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTc > 1:607457/1‑150 (MQ=255) ccAAGGTGTAGCCCAGCGTCAGCCCGAAGCCGGGTATGACGGGGGAGATACGACGCGACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTcg < 2:675507/150‑1 (MQ=255) cAAGGTGTAGCCCAGCGTCAGCCCGAAGCCGGGTATGACGGGGGAGATACGACGCGACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTcgc > 1:83066/1‑150 (MQ=255) tgtAGCCCAGCGTCAGCCCGAAGCCGGGTATGACGGGGGAGATACGACGCGACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGacac > 1:555640/1‑150 (MQ=255) gCCCAGCGTCAGCCCGAAGCCGGGTATGACGGGGGAGATACGACGCGACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAgg > 1:103657/1‑150 (MQ=255) cccGAAGCCGGGTATGACGGGGCAGATACGACGCGACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAg < 1:531200/150‑1 (MQ=255) ccGAAGCCGGGTATGACGGGGGAGATACGACGCGACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAgg > 2:243579/1‑150 (MQ=255) gggggAGATACGACGCGACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTtgt > 2:114433/1‑150 (MQ=255) gggggAGATACGACGCGACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTtgt > 2:114431/1‑150 (MQ=255) gggggAGATACGACGCGACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAAGCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTtgt < 2:83066/150‑1 (MQ=255) gcgACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTtgt > 2:293547/1‑150 (MQ=255) gACATAAATATCCCTGGGTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGc > 2:385508/1‑150 (MQ=255) cccTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGc < 2:636309/150‑1 (MQ=255) cccTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGc < 2:720034/150‑1 (MQ=255) cTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCg > 2:520050/1‑150 (MQ=255) cTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTACTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCg > 1:364195/1‑150 (MQ=255) aaaagcttgcaccgGCCCGGGATTGTTTCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTtgt < 1:664125‑M1/104‑1 (MQ=255) aaaagcttgcaccgGCCCGGGATTGTTTCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTtgt > 2:664125‑M1/15‑118 (MQ=255) ttGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCAcg > 1:365537/1‑150 (MQ=255) aaagcttgcaccgGCCCGGGATTGTTTCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTc > 1:556461‑M1/14‑112 (MQ=255) aaagcttgcaccgGCCCGGGATTGTTTCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTc < 2:556461‑M1/99‑1 (MQ=255) aCGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGc < 2:607457/150‑1 (MQ=255) aCTAAGCCCGGAAACAATCCCGGGCCGGTGCACGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCg > 1:93684/1‑150 (MQ=255) aCTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAg > 1:465196/1‑61 (MQ=255) aCTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAg < 2:465196/61‑1 (MQ=255) CCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATAATCGTCCCTTGGGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCGCTTTCCt > 2:441947‑M1/1‑150 (MQ=255) | CCAGCGGTATCAGCACGATCAGGCTGAGGTACACCAAGGTGTAGCCCAGCGTCAGCCCGAAGCCGGGTATGACGGGGGAGATACGACGCGACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCGCTTTCCT > NC_002947/5896138‑5896408 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |