Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 515,053 | G→C | 38.7% | S188W (TCG→TGG) | ydcA ← | putative rhomboid protease |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 515,053 | 0 | G | C | 38.7% | 9.3 / 2.2 | 20 | S188W (TCG→TGG) | ydcA | putative rhomboid protease |
Reads supporting (aligned to +/- strand): ref base G (3/7); new base C (2/5); total (6/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 |
GTCCTTCCTAGACTTCATGGGTTTTCAATCACGCTGAAAAGATGACAAAATCCTAAAATGAAAACCATTTTAGGATTTTGTCAACAATCTG‑AAGCGGCGGTTAATACCGCCGC‑TC‑TTTTTTTTGCACGCATAAAAA‑CGAAAGAAGAAACCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGTTGATAAACGACATGAGAACGGCAAATGCGAGCAGTGTGATGATCATTTTGGAATGCTCCTGCCCAATT > BS168_reference_genome/514916‑515178 | ggcctgcctAGACTTCATGGGTTTTCAATCCCGCTGAAAAGATCAAAAAAGCCTAAAATGAAAACCATTTTCGGATTTTGTCAACAATCTG‑AAGCGGCGGTTAATACCGCCGC‑TC‑TTTTTTTTGCACGCATAAAAA‑CGa < 6:138315/137‑1 (MQ=21) gtcttGGGTTTTCAAAACCGAATAAAAGATGAAAAAAACCAAAATTAAAAACAATTTTCGGATTTTTCAAACAATCTC‑AAGCGGGGGTTATAACCGCCGC‑TC‑TTTTTTTTGCCCGCATAAAAA‑CCAAAGAAGAACCccg < 8:410440/135‑1 (MQ=255) cATGGGTTTTCAATCACGCTTAAAAGATGACAAAATCCTAAAATTAAAACAATTTTTGTATTTTGTCAACAATCTG‑AAGCGGCGGTTAATACCGCCGC‑TC‑TTTTTTTTGCACGCATAAAAA‑CGAAAGAAGAAACccg < 5:391911/137‑1 (MQ=25) cATGGGTTTTCAATCACGCTGAAAAGATGACAAAATCCTAAAATGAAAACCATTTTAGGATTTTGTCAACAATCTG‑AAGCGGCGGTTAATACCGCCGCGGT‑TTTTTTTTGCA‑GCATAAAAA‑ACAAATGAGAAAAccccc > 3:343019/1‑139 (MQ=17) ggTTTTAATACACGATTAAAACATCAAAAAAACATAAAATCAAACCAATTTTTGTATTTTGTCAACAATCTG‑AGGCGGCGGTTAATACCCCCGC‑TC‑TTTTTTTTGCACGCATAAAAA‑CGAAATAAGAAACCCGCCGCa < 7:317049/138‑1 (MQ=255) ttCAATCACGCTGAAAAGATGACAAAATCCTAAAATGAAAACCATTTTAGGATTTTGTCAACAATCTG‑AAGCGGCGGGTTATAACGCCGC‑TCTTTTTTTTTGTACACAAATAAA‑TAAAAAAAAAAACCCGCCGCaagta > 5:339185/1‑135 (MQ=17) aaaTCACGATTAAAAGATAAAAAAAACATAATATCAAAAAATTTTTTGTTTTTTTTAAACAATCTT‑AAGCGGCGGTTAATACCGCCGC‑TC‑TTTTTTTTGCACGCATAAAAA‑CGAAAGAAGAAACCCGCCGCATAGT‑Cca < 5:141006/138‑1 (MQ=255) acaaaaacaaaaAATAAAAACATTTTTAGTTTTTTGTAAACAATCTG‑AAGCGGCGGTTAATACCGCCGC‑TC‑TTTTTTTTGCACGCATAAAAA‑CGAAAGAAGAAACCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGt < 5:147850/139‑1 (MQ=255) aaaaaCCCAAAAAAAAAAAACATTTTAGTTTTTTTTCAAACATCCGCACGCGGGGTTAAATCCCGCCGC‑TA‑TTTTTTTTCCACCAAAAAAAA‑CCAAAGACGAAACCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGt < 1:96730/138‑1 (MQ=255) taTTTTGAAAACAATCTG‑ACGCGGCGGTTAATACCGCCGC‑TC‑TTTTTTTTCCACGCATAAAAA‑CGAAAGAAGAAACCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGTTGATAAACGACATgag < 5:147457/137‑1 (MQ=255) ttGTCAACAATCTG‑AAGCGGCGGTTAATACCGCCGC‑GC‑GTTTTTTTTGGCGCCAAAAAA‑ACAAAGAAGAAAACCCCCCCCGAGTGCCAAAAAAATGGGCCGACA‑GTTTATATTTGTGTTGTTAAAAGACAAGAGAACgg > 2:95188/1‑139 (MQ=255) taaataaTCTT‑AGGCGGGGTTTATTCCCGCGGC‑TT‑TTTTTTTTCCACTAAAAAAAA‑CTAAATAACAACCCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGTTGATAAACGACATGAGAACGGCa < 7:165206/133‑1 (MQ=255) aaTCTC‑AAGCGGCGTTAAATACCGCCTC‑TC‑TTTTTTTTCCACCCATAAAAA‑CCAACTAACAAACCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGTTGATAAACGACATGAGAACGGCAAATGCg < 7:93762/139‑1 (MQ=255) gcggaggcgtcgTTTAATCCCGCCGC‑TG‑TTTTTTTTGAACGAAAAAAAA‑CGAAATAACAACCCCGCCGCATAGT‑CCAAACAATTGGGCCATCATGTTAATATTGGAGTTGATAAACGACATGAGAACGGCAAATGCGAg < 8:130619/133‑1 (MQ=255) ggTTAATACCGCCGC‑TC‑TTTTTTTTGCCCG‑ATAAAAACCTAAATATGAAACCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATTTTGGAGTTGATAAACGACATGAGAACGGCAAATGCGAGCAGTGtgatga < 8:196369/139‑1 (MQ=255) ccgc‑TT‑TTTTTTTGCCACAAAAAAAAC‑CCAAACAACAACCCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGTTGATAAACGACATGAGAACGGCAAATGCGAGCAGTGTGATGATCATTTTGGaa < 3:304915/139‑1 (MQ=255) ccgc‑TG‑GTTTTTTTTGACGCATAAAAA‑CGGAAGTAGAAACCCGCCGCATAGT‑CCAAACAAATGGGCCAAGATGTTAATATTGGAGTTGATAAACGACATGAGAACGGGAAATGCGAGGAGGGTGATGATCATTTTGGaa > 2:266452/1‑139 (MQ=255) ccgc‑TC‑TTTTTTTTGCACGCATAAAAA‑CGAAAGAAGAAACCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGTTGATAAACGACATGAGAACGGCAAATGCGAGCAGTGTGATGATCATTTTGGaa > 1:54922/1‑139 (MQ=255) ttttttCCACCCATAAAAA‑CCAAAGAAGAAACCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGTTGATAAACGACATGAGAACGGCAAATGCGAGCAGTGTGATGATCATTTTGGAATGCTCCt < 3:175873/138‑1 (MQ=255) aCGCATAAAAA‑CGAAAGAAGAAACCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGTTGATAAACGACATGAGAACGGCAAATGCGAGCAGTGTGATGATCATTTTGGAATGCTCCTGCCCAAtt > 1:409601/1‑138 (MQ=255) | GTCCTTCCTAGACTTCATGGGTTTTCAATCACGCTGAAAAGATGACAAAATCCTAAAATGAAAACCATTTTAGGATTTTGTCAACAATCTG‑AAGCGGCGGTTAATACCGCCGC‑TC‑TTTTTTTTGCACGCATAAAAA‑CGAAAGAAGAAACCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGTTGATAAACGACATGAGAACGGCAAATGCGAGCAGTGTGATGATCATTTTGGAATGCTCCTGCCCAATT > BS168_reference_genome/514916‑515178 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |