Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 2,581,726:1 | +T | 100% | intergenic (‑109/+45) | pstS ← / ← pbpA | phosphate ABC transporter (phosphate binding lipoprotein)/transpeptidase (penicillin‑binding protein 2A) |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 2,581,726 | 1 | . | T | 95.0% | 125.6 / ‑0.2 | 40 | intergenic (‑109/+45) | pstS/pbpA | phosphate ABC transporter (phosphate binding lipoprotein)/transpeptidase (penicillin‑binding protein 2A) |
Reads supporting (aligned to +/- strand): ref base . (1/1); new base T (29/9); total (30/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.42e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
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. |
CATAAGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑G‑TTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCC > BS168_reference_genome/2581591‑2581849 | cATAAGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑G‑ttt > 8:67956/1‑139 (MQ=255) ccTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAA‑GGGTTTGAACCGATTGTTAAGGTGTTGTAAATCA‑GGGC‑GTTTTTTTGTATAAAAAAAACTCCCACATCATGTGAGCGtgttt > 1:269243/1‑139 (MQ=255) tataATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAg > 8:193959/1‑139 (MQ=255) aaTGTGTAATTTGTCCTACACGTAAAGAGTAAAGGTTCTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTATTATTATAGTTa > 7:212077/1‑139 (MQ=255) tgtgTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTAt > 3:311671/1‑138 (MQ=255) ttAATTTTTCATACACTTATTGATTAAAGGTTTTATTTAAAAC‑GGTTTTGAACCTAATGTTAAGGTTTTTTAAATAA‑GTGG‑GTTTTTTTTTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCag < 8:184414/137‑1 (MQ=255) tatttggTCCAACACGTATATATAAATGGTTCTA‑GTAAAAACGGGTTCGACCCTATTTTTATGGTTTTTTAAACCT‑GGGG‑GTTTTTTTTTAAAAAAACAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCaga < 6:162406/132‑1 (MQ=255) tAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCag > 7:102941/1‑137 (MQ=255) ccaacacGTATAGAGTAAAGGTTCTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACTCCCAAATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGCCGTtgt > 1:80333/4‑139 (MQ=255) ccTACACGTATAGAGTAAAGGTTCTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTtgt > 7:184414/1‑139 (MQ=255) ataGAGTAAAGGTTCTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGCGG‑TTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTg > 1:55729/1‑138 (MQ=255) ataGAGTAAAGGTTCTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTg > 4:204962/1‑138 (MQ=255) aaGGTTCTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGt > 8:143591/1‑106 (MQ=255) aaGGTTCTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTTAAAGAACACGTTGTGTTTAGTGCAttttat > 2:298843/1‑136 (MQ=255) aaGGTTCTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTttctt > 3:63297/1‑139 (MQ=255) aaGGTTCTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTttct > 7:138751/1‑138 (MQ=255) aaGGTTCTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATACGAAGCCGTTGTGTTTTCTGCATTttct > 6:200783/1‑138 (MQ=255) aaGGTTCTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTAAAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGt < 7:143591/106‑1 (MQ=255) cTATGTAAAAA‑GGGTTTGAACCGATTGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTGTAGTTATCAGAAGACGGTGTGTGTGCTGCATTTTCTTGAAttt > 6:112543/1‑139 (MQ=255) cTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAAttt > 1:128553/1‑139 (MQ=255) cTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAAttt > 1:88310/1‑139 (MQ=255) cTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAAttt > 8:111884/1‑139 (MQ=255) cTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAAttt > 3:256834/1‑139 (MQ=255) cTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAAtttt > 1:205383/1‑139 (MQ=255) aaaa‑GGGTTTGAACCGAATGTTAAGGTTTTGTACATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCTTAAGACGTTGTGTTTTCGGCATTTTGTTCAATTTTGGCttt > 1:107744/1‑139 (MQ=255) aaa‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTa > 4:51481/1‑139 (MQ=255) gTTTGAACTGAATGTTAAGGTTTTGTAAATCA‑GTCC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTGTTTTTTTAGTTATAAGAAGACGTTGTGTTTTCTGAATTTTATTCAATTTTGTCGTTTTTAtt > 6:33600/1‑139 (MQ=255) gTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCtt < 7:13764/132‑1 (MQ=255) gTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCtt > 8:13764/1‑132 (MQ=255) gTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCtt > 4:91010/1‑139 (MQ=255) gTTTGAACCGAATGTTAAGGTTTTGTAAATAA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAACACGTTGTGTTTTCTGCATTTTCTACAATTTTGTCTTTATTCtt > 3:214789/1‑139 (MQ=255) aaCCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTTTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCt < 7:67956/139‑1 (MQ=255) gAATGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAg > 8:109270/1‑138 (MQ=255) aaTGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTg > 7:2136/1‑132 (MQ=255) aaTGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTg < 8:2136/132‑1 (MQ=255) aTGTTAAGTTTTGTTAAATCAGGGGT‑T‑TTTTTTTTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTAc < 2:88310/139‑1 (MQ=255) aTGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTTTAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTAc < 7:193959/139‑1 (MQ=255) aTGTTAAGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAACACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTAc < 5:159362/139‑1 (MQ=255) aGGTTTTGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACcttct < 8:102941/138‑1 (MQ=255) ttGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTcc > 8:66472/1‑139 (MQ=255) ttGTAAATCA‑GTGC‑GTTTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTcc > 4:253448/1‑139 (MQ=255) | CATAAGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAA‑GGGTTTGAACCGAATGTTAAGGTTTTGTAAATCA‑GTGC‑G‑TTTTTTGTATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCC > BS168_reference_genome/2581591‑2581849 |
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: ‑ |