Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 3,774,612 | C→A | 27.7% | intergenic (‑51/+43) | ywmE ← / ← ywmD | conserved protein of unknown function/putative exported protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 3,774,612 | 0 | C | A | 27.7% | 92.1 / 6.0 | 50 | intergenic (‑51/+43) | ywmE/ywmD | conserved protein of unknown function/putative exported protein |
Reads supporting (aligned to +/- strand): ref base C (31/4); new base A (2/13); total (33/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.68e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CAAAAAAGAAGGCTGCCAGTTTTAGAAGAACACCAAGCAATATAAATGCAACGGTGGCTATAAGAAAAA‑TCATTCCGAATAATTTCATA‑TGTTCAACCCCTTTTTAGGGTATCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGAACTCATCCGCACTGTTCGCTGAAATATATTCTCCGCCGCCAGC > BS168_reference_genome/3774474‑3774741 | cAAAAAAGAAGGCTGCCAGTTTTAGAAGAACACCAAGCAATATAAATGCAACGGTGGCTATAAGAAAAA‑TCATTCCGAATAATTTCATA‑TGTTCAACCCCTTTTTAGGGTATCATTCCCGCCCAAAC‑TGTTTTTAAAcc < 8:294969/139‑1 (MQ=255) agGCTGCCAGTTTTAGAAGAACACCAAGCAATATAAATGCAACGGTGGCTATAAGAAAAA‑TCATTCCGAATAATTTCATA‑TGTTCAACCCCTTTTTAGGGTATCATTCCCGCCCAAAC‑TGTTGTTAAACCAATaaaa > 3:48996/1‑137 (MQ=255) gCTGCCAGATTTAGAAGAACACCAAGCAATATAAATGCAACGGTGGCTATAAGAAAAA‑TCATTCCGATTAATTTCATA‑TGTTCAACACCTTTTTAGGGTATCATTCCCGCCCAGAC‑TGTTGTTAAACAAATAAAA‑AAgg > 2:51469/1‑139 (MQ=255) ccAGTTTTAGAAGAACACCAAGCAATATAAATGCAACGGTGGCTATAAGAAAAA‑TCATTCCGAATAATTTCATA‑TGTTCAACCCCTTTTTAGGGTATCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGcc > 5:313371/1‑139 (MQ=255) caaCACAAAGCAATATAAATCCAAGGGTGGATAAAAGAAAAA‑TCATACCCAATAATTAAATA‑TGTCAAACCCATTTTTTTGTTATCATCCCCGCCCAAAT‑TTTTTTAAAACAAAAAAAA‑AACGCGCCTGGGGGAGCGcc < 6:313371/138‑1 (MQ=255) aaCACAAATCAATATAATTCCAAGGGGTTATATAATAAAAA‑TAATTCCTAATATTTTCTTT‑TGTACACCGATTTTTTTGGGTATCA‑TCCCCCCCAAACGTTTTTTAAAACAAATAAAA‑ACCGCGCCTGGGGGAGCGCCt < 2:321821/139‑1 (MQ=255) acacCAAGCAATATAAATGCAACGGTGGCTATAAGAAAAA‑TCATTCCGAATAATTTCATA‑TGTTCAACCCCTTTTTAGGGTATCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCt > 8:83915/1‑138 (MQ=255) acacCAAGCAATATAAATGCAACGGTGGCTATAAGAAAAA‑TCATTCCGAATAATTTCATA‑TGTTCAACCCCTTTTTAGGGTATCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTg > 6:34467/1‑139 (MQ=255) accAATATAAATCCACCTGTGGCTATAACAAAAA‑TCATTCCTAATAATTTCATA‑TGTACAACCCCTTTTTCGGGTATCATGCCCGCCCAAAC‑TGTTTTTAACCCAAAAAAA‑ACGGCGCCTGGGGGAGCGCCTGATCttt < 7:44660/137‑1 (MQ=255) acaaaTATAAATGCAAGGGTGCCAAAAAAAAAAA‑AAATTCATAATAATTTAATA‑TGTAAAACCCCTTTTTGGGCTATACTACCCCCCCAATT‑TTTTTTAAAACAAAAAAAA‑ACCGCGCCTGGGGGAGCGCCTGATCttt < 4:52601/136‑1 (MQ=255) aGCAATATAAATGCAACGGTGGCTATAAGAAAAA‑TCATTCCGAATAATTTCATA‑TGTTCAACCCCTTTTTAGGGTATCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGTGGAGCGCCTGATCta > 8:390867/1‑137 (MQ=255) aGCAATATAAATGCAACGGTGGCTATAAGAAAAA‑TCATTCCGAATAATTTCATA‑TGTTCAACCCCTTTTTAGGGTATCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCttt > 3:99368/1‑139 (MQ=255) aaTAAAAATGCAACGGGGGATATAATAAAAATTCATTCCTAATAATTTCATA‑TGTTCAACCCCTTTTTAGGTTATCATTCCCGCCAAAAC‑TTTTTTTAAACAAAAAAAA‑ACCGCGCCTGGGGGAGCGCCTGATCTTTAg < 8:337354/139‑1 (MQ=255) aaaTGCAACGGTGGCTATAAGAAAAA‑TCATTCCGAATAATTTCATA‑TGTTCAACCCCTTTTTAGGGTATCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCt > 3:1234/1‑139 (MQ=255) aTGCAACGTGTGCTATAACAAAAA‑CATTTCCCAATATTTAAAAA‑TTTACAACCCCTTTTTCGGGTATCAGCCCCCCCAAAAA‑TTTTTTAAAACAAAAAAAA‑CCCGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTa < 1:331081/139‑1 (MQ=255) cAACGGTGGCTATAAGAAAAA‑TCATTCCGAATAATTTCATA‑TGTTCAACCCCTTTTTAGGGTATCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCATTACt > 1:172400/1‑138 (MQ=255) cAACGGTGGCTATAAGAAAAA‑TAATTCCGAATAATTTCATA‑TGTTCAACCCCTTTTTAGGGTATCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACt > 4:122935/1‑138 (MQ=255) aCGGTGGCTATAAGAAAAA‑TCATTCCGAATAATTTCATA‑TGTTCAACCCCTTTTTAGGGTATCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTAcctg > 1:185166/1‑135 (MQ=255) aCGGTGGCTATAAGAAAAA‑TCATTCCGAATAATTTCATA‑TGTTCAACCCCTTTTTAGGGTATCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAg > 6:43957/1‑139 (MQ=255) aCGGTGGCTATAAGAAAAA‑TCATTCCGAATAATTTCATA‑TGTTCAACCCCTTTTTAGGGTATCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAg > 6:391874/1‑139 (MQ=255) aCGGTGGCTATAAGAAAAA‑TCATTCCGAATAATTTCATA‑TGTTCAACCCCTTTTTAGGGTATCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAg > 6:349768/1‑139 (MQ=255) aCGGTGGCTATAAGAAAAA‑TCATTCCGAATAATTTCATA‑TGTTCAACCCCTTTTTAGGGTATCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAg > 2:292842/1‑139 (MQ=255) aCGGTGGCTATAAGAAAAA‑ACATTCCGAATAATTTCATA‑TGTTCAACCCCTTTTTAGGGTAACATACCCGCCCAAAC‑TGTTATTAAACAAATAAAA‑AAGGAGCATGGGGGAGCGCCTGATCTTAAGGTACCTTACTCAg > 1:190764/1‑139 (MQ=255) aGAAAAA‑TAATTCATAAAAATTTAATA‑TGTACAACCCATTTTTAGGTTATCTTTCCCGCCCAAAG‑TGTTTTTAAACCAATAAAA‑ACCGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTtgtg < 5:195978/139‑1 (MQ=255) aaa‑aCATTCCCTAT‑ATTTCATTTTGTTCAACCCATTTTATTAGAACCATGCCCCCCCAAAT‑TTTTTTTAAACAAAAAAAA‑ACCGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGcttc < 8:102870/139‑1 (MQ=255) a‑TCATTCCGAATAATTTCATA‑TGTTCAACCCCTTTTTAGGGTATCATTCCCGCCCAAAC‑TGTTTTTAAACCATTAAAA‑AAGTCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGcttctt > 4:134267/1‑139 (MQ=255) ttattttAATA‑TG‑TCAACCCTTTTTTTTGTTAAAAATCCCCCCAAAAT‑TTTTTTTAAAAAAAAAAAACCCCGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGc < 4:48996/135‑1 (MQ=255) ataatTTAATA‑TGTTAACCCCCTTTTTCGGGTATCATTCCCGCCCAAAG‑TGTTTTAAAACAAAAAAAA‑ACCGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGc < 7:213425/139‑1 (MQ=255) aatttaatt‑tGTTAAACCCATTTTTTGGGAAACATCCCCCCCCAACC‑TTTTTTAAAACAAAAAAAA‑ACCGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCt < 8:312529/139‑1 (MQ=255) ta‑tGTTCACCCCATTTTTTGGGTATCATTCCCGCCCAATC‑TGTTTTAAAACAAAAAAAA‑CCCGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGaa < 7:83915/139‑1 (MQ=255) ta‑tGTTCAACCCCTTTTTAGGGTATCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTGTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGaa > 5:143367/1‑139 (MQ=255) ta‑tGTTCAACCCATTTTTGGGTTAAAATACCCCCCCAATC‑TGTTTTAAAACAAAAAAAA‑ACCGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGaa < 5:301905/139‑1 (MQ=255) ccggTTCAACCCATTTTTTGGTTATCATACCCGCCAAAAC‑TTTTTTTAAACAAAAAAAA‑ACCGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGaa < 6:306932/136‑1 (MQ=255) tCAACCCCTTTTTAGGGTATCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGCGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGAACTCAt > 2:179664/1‑139 (MQ=255) ccttttttCGGTTATCAACCCCCCCCAACG‑TTTTTTAAAACAAAAAAAA‑CCCGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGAACTCATCCGCa < 3:122935/136‑1 (MQ=255) cacTTTTTCGGGTATCATACCCGCCAAAAC‑TTTTTTTAAACCAATAAAA‑ACGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGAACTCATCCGCa < 2:281143/137‑1 (MQ=255) ttttAGGGTATCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGAACTCATCCGCAcctg > 5:306932/1‑136 (MQ=255) ggTATCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGCTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGAACTCATCCGCACTGTTCGCt > 5:323098/1‑138 (MQ=255) aTCTTTCCCGCCCAATC‑TGTTTTTATACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTACCTTACTCAGTAAATATTAGAGCTTCTTTTTCCCAAGCCTGTTCGAACTCATCCGCACTGTTCGCTGaaa > 7:126399/1‑139 (MQ=255) tCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGAACTCATCCGCACTGTTCGCTGaaa > 5:360869/1‑138 (MQ=255) tCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGAACTCATCCGCACTGTTCGCTGaaa > 8:332626/1‑138 (MQ=255) ccGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGAACTCATCCGCACTGTTCGCTGAAATATATtc > 3:45876/1‑139 (MQ=255) ccGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGAACTCATCCGCACTGTTCGCTGAAATATATtc > 4:246658/1‑139 (MQ=255) gCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGAACTCATCCGCACTGTTCGCTGAAATATATtctc > 2:402312/1‑139 (MQ=255) gCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGCCGCCTGTGGGAGCGCCTCATCTTTATGTTAATTAATCACTAAATATTTGTGCTACTTTTTCCCAAGCCTGTAAGAACACATCCGCACAGTTCCCTGAAATAAATtccc > 2:199425/1‑137 (MQ=255) gCCCAAAC‑TGTTTATAAACCAATAAAA‑ACCGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGAACTCATCCGCACTGTTCGCTGAAATATATtctc < 3:329598/139‑1 (MQ=255) aaC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGAACTCATCCGCACTGTTCGCTGAAATATATTCTccgccg > 2:19109/1‑139 (MQ=255) aaC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGAACTCATCCGCACTGTTCGCTGAAATATATTCTccgccg > 3:249525/1‑139 (MQ=255) aC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGAACTCATCCGCACTGTTCGCTGAAATATATTCTCcgccgc > 5:39483/1‑139 (MQ=255) tttttAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGAACTCATCCGCACTGTTCGCTGAAATATATTCTCCGCCGCCAGc > 7:178924/1‑139 (MQ=255) tttttAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGAACTCATCCGCACTGTTCGCTGAAATATATTCTCCGCCGCCAGc > 5:198114/1‑139 (MQ=255) | CAAAAAAGAAGGCTGCCAGTTTTAGAAGAACACCAAGCAATATAAATGCAACGGTGGCTATAAGAAAAA‑TCATTCCGAATAATTTCATA‑TGTTCAACCCCTTTTTAGGGTATCATTCCCGCCCAAAC‑TGTTTTTAAACCAATAAAA‑AAGGCGCCTGGGGGAGCGCCTGATCTTTAGGTTCCTTACTCAGTAAATTTTTGTGCTTCTTTTTCCCAAGCCTGTTCGAACTCATCCGCACTGTTCGCTGAAATATATTCTCCGCCGCCAGC > BS168_reference_genome/3774474‑3774741 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |