| 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 | 90.9% | 65.6 / 0.3 | 22 | intergenic (‑109/+45) | pstS/pbpA | phosphate ABC transporter (phosphate binding lipoprotein)/transpeptidase (penicillin‑binding protein 2A) |
| Reads supporting (aligned to +/- strand): ref base . (0/2); new base T (17/3); total (17/5) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.33e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| 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. | |||||||||||
CAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑G‑TTTTTTGT‑ATAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCT > BS168_reference_genome/2581600‑2581850 | cAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑G‑tttttttgtat > 5:66462/1‑133 (MQ=255) aTTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GTTTTTTTGT‑ATaa > 6:223339/1‑138 (MQ=255) tAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GTTTTTTTGT‑ATAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTTAGTTAtaat > 7:164694/1‑134 (MQ=255) gTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GTTTTTTTGT‑ATAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTTAGTTATCAGAAGACGtt > 5:165978/1‑139 (MQ=255) tataGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GTTTTTTTGT‑ATAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTg > 7:169192/1‑139 (MQ=255) ataGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GTTTTTTTGT‑ATAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTTAGTTTTCAGAAGACGTTGTGTTTTCTg > 3:117962/1‑138 (MQ=255) ataGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACAGAATGTTAAGGTTTTGTAAATCAGTGC‑GTTTTTTTGT‑ATAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTTAGTTATGAGAAGGCGTTGTGTTTtgtgg > 3:123228/1‑138 (MQ=255) aaGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GTTTTTTTGT‑ATAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTttctt > 7:63830/1‑139 (MQ=255) aGGTTCTATGTAAAAAGGGTTTGACCAGAATGATAAGGTTTTGTAAATCAGTGC‑GTTTTTTTGT‑ATAAAAAAAACGCACACATCATGTGAGCG‑TTTTTTTTTTTTTATAAGAATCCGTTGTGTTTTCTGCTTTttc > 8:62593/1‑136 (MQ=255) aGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GTTTTTTTGT‑ATAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTttctt > 8:124553/1‑138 (MQ=255) cTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GTTTTTTTGT‑ATAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAAttt > 1:154728/1‑139 (MQ=255) cTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GTTTTTTTGT‑ATAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAAttt > 5:152565/1‑139 (MQ=255) cTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GTTTTTTTGT‑ATAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAAttt > 6:152679/1‑139 (MQ=255) tGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GTTTTTTTGT‑AAAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAAttt < 5:152679/136‑1 (MQ=255) aaaGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GTTTTTTTGT‑ATAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTa > 1:134992/1‑139 (MQ=255) gTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GTTTTTTTGT‑ATAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCtt > 1:214697/1‑139 (MQ=255) gTTTCAACCTATTTTTATGGTTTTGAAAATCCGTGCTT‑TTTTTTTT‑ATAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCtt < 6:152565/139‑1 (MQ=255) agAATGTTAAGGTTTTGTAAATCCGGGC‑GTTTTTTTGT‑AAAAAAAACACGCCCACATCATGTGAGCG‑TTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAg < 6:117643/138‑1 (MQ=255) gAATGTTAAGGTTTTGTAAATCAGTGC‑GTTTTTTTGT‑ATAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTTATTTATCAGAAGAAGTTGTTTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGat > 3:171781/1‑137 (MQ=255) aaTGTTAAGGTTTTGTAAATCAGTGC‑GTTTTTTTGT‑ATAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTa > 1:116114/1‑139 (MQ=255) aTGTTACGGTTTTTTAAATCAGGGC‑GTTTTTTTGT‑AAAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTTAGTTATCAGAAGACGTGGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTAc < 5:52029/139‑1 (MQ=255) tGTAAATCAGTGC‑GTTTTTTTGT‑ATAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTcc > 6:115681/1‑139 (MQ=255) ttAAATCAGTGC‑T‑TTTTTTTTAAAAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCt < 4:117962/137‑1 (MQ=255) | CAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑G‑TTTTTTGT‑ATAAAAAAAACGCCCACATCATGTGAGCG‑TTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCT > BS168_reference_genome/2581600‑2581850 |
| 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: ‑ |