Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 2,027,468 | Δ1 bp | 100% | intergenic (+118/‑180) | yodD → / → yedP | hypothetical protein/conserved hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 2,027,468 | 0 | T | . | 81.2% | 82.8 / 21.8 | 32 | intergenic (+118/‑180) | yodD/yedP | hypothetical protein/conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (3/3); new base . (15/11); total (18/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. |
AAGAAGCACTTAAAAAATTCGTTACACCAGGAAATCTGATGTGTTCATCACCTTATCCGCAATTTTTTTCGCTGACAAGAAAATATTCGCCATATGAATGATTAACGCTACTTATGAATAGAAATGTGA > W3110S.gb/2027406‑2027534 | aagaagCACTTAAAAAATTCGTTACACCAGGAAATCTGATGTGTTCATCACCTTATACGCAAttttttcgc > 1:294646/1‑68 (MQ=255) agCACTTAAAAAATTCGTTACACCAGGAAATCTGATGTGTTCATCACCTTATCCGCAA‑TTTTTTCGCTGa > 1:1265284/1‑70 (MQ=255) agCACTTAAAAAATTCGTTACACCAGGAAATCTGATGTGTTCATCACCTTATCCGCAA‑TTTTTTCGCTGa > 1:2297778/1‑70 (MQ=255) cTTAAAAAATTCGTTACACCAGGAAATCTGATGTGTTCATCACCTTATCCGCAA‑TTTTTTCGCTGACAAGa < 1:2270606/71‑1 (MQ=255) tAAAAAATTCGTTACACCAGGAAATCTGATGTGTTCATCACCTTATCCGCAA‑TTTTTTCGCTGACAAGaa > 1:1238161/1‑70 (MQ=255) aaTTCGTTACACCAGGAAATCTGATGTGTTCATCACCTTATCCGCAAtttt < 1:899705/51‑1 (MQ=255) ttCGTTACACCAGGAAATCTGATGTGTTCATCACCTTATCCGCAA‑TTTTTTCGCTGACAAGAAAATATTCg > 1:1396876/1‑71 (MQ=255) cAGGAAATCTGATGTGTTCATCACCTTATCTGCAA‑TTTTTTCGCTGACAAGAAAATATTCGCCATatgaat < 1:1852425/71‑1 (MQ=255) aGGAAATCTGATGTGTTCATCACCTTATCCGCAA‑TTTTTTCGCTGACAAGAAAATATTCGCCATatgaa > 1:1102506/1‑69 (MQ=255) ggAAATCTGATGTGTTCATCACCTTATCCGCAA‑TTTTTTCGCTGACAAGAAAATATTCGCCa < 1:1096271/62‑1 (MQ=255) aaTCTGATGTGTTCATCACCTTATCCGCAATTTTTTTCGCTGACAAGAAAATATTCGCCATATGAATGAtt < 1:1911731/71‑1 (MQ=255) aaTCTGATGTGTTCATCACCTTATCCGCAATTTTTTTCGCTGACAAGAAAATATTCGCCATATGAATGAt > 1:930887/1‑70 (MQ=255) aaTCTGATGTGTTCATCACCTTATCCGCAA‑TTTTTTCGCTGACAAGAAAatat > 1:2472013/1‑53 (MQ=255) aaTCTGATGTGTTCATCACCTTATCCGCAA‑TTTTTTCGCTGACAAGAAAATATTCGCCATATGAATGAtt > 1:1093778/1‑70 (MQ=255) cTGATGTGTTCATCACCTTATCCGCAATTTTTTTCGCTGa > 1:15268/1‑40 (MQ=255) aTGTGTTCATCACCTTATCCGCAA‑TTTTTTCGCTGACAAGAAAATATTCGCCATATGAATGATTAACGc < 1:562102/69‑1 (MQ=255) tgtTCATCACCTTATCCGCAATTTTTTTCGCTGACAAGAAAATATTCGCCATATGAATGATTAACGCTACt < 1:496174/71‑1 (MQ=255) tgtTCATCACCTTATCCGCAA‑TTTTTTCGCTGACAAGAAAATATTCGCCATatga < 1:2084821/55‑1 (MQ=255) tgtTCATCACCTTATCCGCAA‑TTTTTTCGCTGACAAGAAAATATTCGCCATATGAATGATTAACGCTACt < 1:1355923/70‑1 (MQ=255) ttCATCACCTTATCCGCAATTTTTTTCGCTGACAAGaaaa > 1:2113595/1‑40 (MQ=255) ttCATCACCTTATCCGCAA‑TTTTTTCGCTTACAAGAAAATATTCGCCATATGAATGATTAACGCTACTTAt > 1:2691712/1‑71 (MQ=255) ttCATCACCTTATCCGCAA‑TTTTTTCGCTGACAAGAAAATATTCGCCATATGAATGATTaa > 1:2276428/1‑61 (MQ=255) ttCATCACCTTATCCGCAA‑TTTTTTCGCTGACAAGAAAATATTCGCCATATGAATGATTAACGCTACt < 1:1402535/68‑1 (MQ=255) ttCATCACCTTATCCGCAA‑TTTTTTCGCTGACAAGAAAATATTCGCCATATGAATGATTAACGCTACTTAt > 1:784548/1‑71 (MQ=255) catcaCCTTATCCGCAA‑TTTTTTCGCTGACAAGAAAATATTCGCCa > 1:842874/1‑46 (MQ=39) aCCTTATCCGCAA‑TTTTTTCGCTGACAAGAAAATATTCGCCATATGAATGATTAACGCTACTTATGaa > 1:1412500/1‑68 (MQ=255) ccTTATCCGCAA‑TTTTTTCGCTGACAAGAAAATATTCGCCATATGAATGATTAACGCTACTTATGAATAGa < 1:1155005/71‑1 (MQ=255) cTTATCCGCAATTTTTTTCGCTGACAAGAAAATATTCGCCATATGAATGATTAACGCTACTTATGAATAGa < 1:1198308/71‑1 (MQ=255) ttATCCGCAA‑TTTTTTCGCTGACAAGAAAATATTCGCCATATGAATGATTAACg > 1:2464420/1‑54 (MQ=255) tacGCAA‑TTTTTTCGCTGACAAGAAAATATTCGCCATATGAATGAtt < 1:1443169/45‑1 (MQ=38) tacGCAA‑TTTTTTCGCTGACAAGAAAATATTCGCCATATGAATGAtt < 1:1446334/45‑1 (MQ=38) tCCGCAA‑TTTTTTCGCTGACAAGAAAATATTCGCCATATGAATGAtt < 1:1327471/47‑1 (MQ=39) cGCAA‑TTTTTTCGCTGACAAGAAAATATTCGCCATATGAATGATTAACGCTACTTATGAATAGAAATGTGa > 1:694574/1‑71 (MQ=255) cGCAA‑TTTTTTCGCTGACAAGAAAATATTCGCCATATGAATCATTAACGCTACTTATGAATAGAAATGTGa > 1:2117093/1‑71 (MQ=255) | AAGAAGCACTTAAAAAATTCGTTACACCAGGAAATCTGATGTGTTCATCACCTTATCCGCAATTTTTTTCGCTGACAAGAAAATATTCGCCATATGAATGATTAACGCTACTTATGAATAGAAATGTGA > W3110S.gb/2027406‑2027534 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |