| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 1,607,310 | C→T | intergenic (+21/‑36) | lsrG → / → tam | (4S)‑4‑hydroxy‑5‑phosphonooxypentane‑2,3‑dione isomerase/trans‑aconitate 2‑methyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,607,310 | 0 | C | T | 95.0% | 59.7 / ‑3.1 | 20 | intergenic (+21/‑36) | lsrG/tam | (4S)‑4‑hydroxy‑5‑phosphonooxypentane‑2,3‑dione isomerase/trans‑aconitate 2‑methyltransferase |
| Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (13/6); total (13/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.50e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01 | |||||||||||
AACCACCTGGGCTCTGTACAGGAAGAAGGCAATTTGCGCTTCGATGTCTTACAGGACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAACGCGCTGACGGGAGTAAAAAAATGTCTGACTGGAACCCCTCTTTATATCTACACTTTTCCGCTGAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCT > NC_000913/1607071‑1607552 | aaCCACCTGGGCTCTGTACAGGAAGAAGGCAATTTGCGCTTCGATGTCTTACAGGACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATTTACAATTgggg > 1:55775/1‑251 (MQ=255) aGGAAGAAGGCAATTTGCGCTTCGATGTCTTACAGGACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATTTACAATTGGGGTAACGCGCTGACGGGAGt > 2:96770/1‑250 (MQ=255) aGGAAGAAGGCAATTTGCGCTTCGATGTCTTACAGGACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATTTACAATTGGGGTAACGCGCTGACGGGAGt > 2:129738/1‑250 (MQ=255) tttGCGCTTCGATGTCTTACAGGACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATTTACAATTGGGGTAACGCGCTGACGGGAGTAAAAAAATGTCTg > 1:233908/1‑250 (MQ=255) ccGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATTTACAATTGGGGTAACGCGCTGACGGGAGTAAAAAAATGTCTGACTGGAACCCCTCTTTATATCTacac < 1:120981/251‑1 (MQ=255) aCGAAGACTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATTTACAATGGGGGTAACGCGCTGACGGGAGTAAAAAAATGTCTGACTGGAACCCCTCTTTATATCTACACTTTTCCGCTGAACGATCGCGTCcggcgg < 1:172286/251‑1 (MQ=255) aaaGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATTTACAATTGGGGTAACGCGCTGACGGGAGTAAAAAAATGTCTGACTGGAACCCCTCTTTATATCTACACTTTTCCGCTGAACGATCGCGTCCGGCGGTGGAGCTGCtt < 2:25693/251‑1 (MQ=255) gATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAACGCGCTGACGGGAGTAAAAAAATGTCTGACTGGAACCCCTCTTTATATCTACACTTTTCCGCTGAACGATCGCGTCCGGCGGTGGAGCTGCTTGcc < 2:113659/251‑1 (MQ=255) ggCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATTTACAATTGGGGTAACGCGCTGACGGGAGTAAAAAAATGTCTGACTGGAACCCCTCTTTATATCTACACTTTTCCGCTGAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTgg < 1:78983/250‑1 (MQ=255) gtCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATTTACAATTGGGGTAACGCGCTGACGGGAGTAAAAAAATGTCTGACTGGAACCCCTCTTTATATCTACACTTTTCCGCTGAACGAGCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGccc < 2:25350/250‑1 (MQ=255) aCCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATTTACAATTGGGGTAACGCGCTGACGGGAGTAAAAAAATGTCTGACTGGAACCCCTCTTTTTATCTACACTTTTCCGCTGAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTaca > 2:172128/1‑250 (MQ=255) ccGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATTTACAATTGGGGTAACGCGCTGACGGGAGTAAAAAAATGTCTGACTGGAACCCCTCTTTATATCTACACTTTTCCGCTGAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTacaac > 1:46538/1‑251 (MQ=255) cgcgTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATTTACAATTGGGGTAACGCGCTGACGGGAGTAAAAAAATGTCTGACTGGAACCCCTCTTTATATCTACACTTTTCCGCTGAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGtt > 1:61730/1‑251 (MQ=255) cgTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATTTACAATTGGGGTAACGCGCTGACGGGAGTAAAAAAATGTCTGACTGGAACCCCTCTTTATATCTACACTTTTCCGCTGAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTgg > 1:35374/1‑251 (MQ=255) aaaaCGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATTTACAATTGGGGTAACGCGCTGACGGGAGTAAAAAAATGTCTGACTGGAACCCCTCTTTATATCTACACTTTTCCGCTGAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGc < 1:96770/251‑1 (MQ=255) aCGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATTTACAATTGGGGTAACGCGCTGACGGGAGTAAAAAAATGTCTGACTGGAACCCCTCTTTATATCTACACTTTTCCGCTGAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGc > 1:96141/1‑251 (MQ=255) gCCGTGAGGCGAATTTATCAATTTTATTTACAATTGGGGTAACGCGCTGACGGGAGTAAAAAAATGTCTGACTGGAACCCCTCTTTATATCTACACTTTTCCGCTGAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATACCAGGCATCGACTc > 2:84166/1‑251 (MQ=255) tttATCAATTTTATTTACAATTGGGGTAACGCGCTGACGGGAGTAAAAAAATGTCTGACTGGAACCCCTCTTTATATCTACACTTTTCCGCTGAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCgatg > 1:222611/1‑251 (MQ=255) tATCAATTTTATTTACAATTGGGGTAACGCGCTGACGGGAGTAAAAAAATGTCTGACTGGAACCCCTCTTTATATCTACACTTTTCCGCTGAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCgatgat > 1:90575/1‑251 (MQ=255) aaTTTTATTTACAATTGGGGTAACGCGCTGACGGGAGTAAAAAAATGTCTGACTGGAACCCCTCTTTATATCTACACTTTTCCGCTGAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCCTCGACTCGTCTCCGGCGCTGATTGCt > 2:228266/1‑251 (MQ=255) | AACCACCTGGGCTCTGTACAGGAAGAAGGCAATTTGCGCTTCGATGTCTTACAGGACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAAGACGCAGTGGCGTTCCATAAAACCACGCCCCACTACAAAACCTGTGTCGCGAAACTGGAATCTTTAATGACCGGGCCGCGTAAAAAACGTCTGTTCAATGGTTTGATGCCGTGAGGCGAATTTATCAATTTTATCTACAATTGGGGTAACGCGCTGACGGGAGTAAAAAAATGTCTGACTGGAACCCCTCTTTATATCTACACTTTTCCGCTGAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCT > NC_000913/1607071‑1607552 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |