Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 1,902,269:1 | +T | 100% | coding (14/342 nt) | yebY ← | hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,902,269 | 1 | . | T | 92.0% | 78.3 / 0.7 | 25 | coding (14/342 nt) | yebY | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base . (1/1); new base T (12/11); total (13/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 = 8.95e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. |
GCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑TTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTCGAGTTTAGTCATCCACACAAACAAGGGTGCCGCACGATCTTCATC > REL606/1902125‑1902405 | gCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑tttt‑t > 1:596292/1‑150 (MQ=255) cATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATc > 2:497155/1‑150 (MQ=255) cATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTACAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATc > 2:438805/1‑150 (MQ=255) tctTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAgg < 2:598050/150‑1 (MQ=255) cccATTTGTCTTTACCCACTTCAAAACGGCTAACGGTAATCCCCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGAGAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCACCAAAACGATTCGAGAGTctacagc > 1:164484/1‑144 (MQ=255) tGTCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCtg > 2:76429/1‑150 (MQ=255) gTCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCtga > 2:474915/1‑150 (MQ=255) gTCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGACtga > 1:347241/1‑150 (MQ=255) tCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCtgat > 1:571646/1‑150 (MQ=255) tCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCtgat > 2:614072/1‑150 (MQ=255) tttACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCtgatga < 2:510422/150‑1 (MQ=255) ttCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGta < 1:495747/150‑1 (MQ=255) aaCGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTAcc < 1:76204/150‑1 (MQ=255) aCGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCg < 2:348235/150‑1 (MQ=255) gCTGACAGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCAc < 1:295687/150‑1 (MQ=255) gTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGc > 1:125304/1‑150 (MQ=255) ccgccAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑TTTT‑TTCTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTCGAGTTTAGTc < 1:248317/150‑1 (MQ=255) cTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTCGAGTTTAGTCATCCACacaaac > 2:23479/1‑150 (MQ=255) cTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTCGAGTTTAGTCATCCACacaaac > 2:259637/1‑150 (MQ=255) cTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑TTTTGTT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTCGACTTTACTCATCCACgcaaac > 1:431862/1‑150 (MQ=255) gAACTGGTGAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTCGAGTTTAGTCATCCACACAAACAAggg < 1:258897/150‑1 (MQ=255) agagacaggAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTCGAGTTTAGTCATCCACACAAACAAggg < 1:313451/142‑1 (MQ=255) aaCTGGTGAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTCGAGTTTAGTCATCCACACAAACAAGGGt < 1:23407/150‑1 (MQ=255) gAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTCGAGTTTAGTCATCCACACAAACAAgggctgtctct > 2:314339/1‑142 (MQ=255) gAGTAACAGAAACGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTCGAGTTTAGTCATCCACACAAACAAGGGTGCCGCAc < 2:573316/150‑1 (MQ=255) aaCGCCAGAATACTTTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTCGAGTTTAGTCATCCACACAAACAAGGGTGCCGCACGATCTTCATc < 2:143774/150‑1 (MQ=255) | GCATTACCCGGTCGGCAAGTCAGCATCACCTCTTCGCGATTAAACGCCCATTTGTCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑TTTT‑TT‑TCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTCGAGTTTAGTCATCCACACAAACAAGGGTGCCGCACGATCTTCATC > REL606/1902125‑1902405 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |