Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 515,043 | G→C | 28.4% | C191W (TGC→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,043 | 0 | G | C | 28.4% | 14.2 / 2.5 | 24 | C191W (TGC→TGG) | ydcA | putative rhomboid protease |
Reads supporting (aligned to +/- strand): ref base G (9/8); new base C (2/5); total (11/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.86e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01 |
GTCCTTCCTAGACTTCATGGGTTTTCAATCACGCTGAAAAGATGACAAAATCCTAAAATGAAAACCATTTTAGGATTTTGTCAA‑CAATCTGAAGCGGCGGTTAAT‑A‑CCGCC‑GC‑TC‑TTTTTTTTGCA‑CGC‑ATAAAAACGAAAGAAGAAACCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGTTGATAAACGACATGAGAACGGCAAATGCGAGCAGTGTGATGATCATTTTGGAATGCTCCTGCCC > BS168_reference_genome/514916‑515174 | gTCCTTCCTAGACTTCATGGGTTTTCAATCACGCTGAAAAGATGACAAAATCCTAAAATGAAAACCATTTTAGGATTTTGTCAA‑CAATCTGAAGCGGCGGTTAAT‑A‑CCGCC‑GC‑TG‑TTTTTTTTTGA‑AGC‑Aaaaaaaac > 7:114646/1‑136 (MQ=17) gTCCTTCCTAGACTTCATGGGTTTTCAATCACGCTGAAAAGATGACAAAATCCTAAAATGAAAACCATTTTAGGATTTTGTCAA‑CAATCTGAAGCGGCGGTTAAT‑A‑CCGCC‑CC‑GG‑GTTTTTTTGGAGCGC‑GAAAAAACg > 1:288841/1‑139 (MQ=17) gTCCTTCCTAGACTTCATGGGTTTTCAATCACGCTGAAAAGATGACAAAATCCTAAAATGAAAACCATTTTAGGATTTTGTCAA‑CAATCTGAAGCGGCGGTTAAA‑A‑ACGCC‑GC‑GA‑TTTTTTTTTAA‑CCC‑Aaaaaaaaa > 6:91321/1‑136 (MQ=17) ccttcctAGACTTCATGGGTTTTCAATCACGCTGAAAAGATGACAAAATCCTAAAATGAAAACCATTTTAGGATTTTGTCAA‑CAATCTGAAGCGGCGGTTAAT‑A‑CCCCC‑GC‑GC‑GTTTTTTTGCC‑CGC‑AAAAAAACGaa > 8:316234/1‑138 (MQ=18) ccttcctAGACTTCATGGGTTTTCAATCACGCTGAAAAGATGACAAAATCCTAAAATGAAAACCATTTTAGGATTTTGTCAA‑CAATCTGAAGCGGCGGGTAATAA‑CTGCG‑CC‑GC‑TTTTTTTTGGA‑‑GC‑Aaaaaaaacca > 3:420899/1‑134 (MQ=14) tattGGTTTTACAATCAGGCAGAAATGAGGCCAAAATCAAAAATTCAAAACCTTTTTCTTTTTTTGTCAA‑CAATCTGAAGCGGCGGTTAAT‑A‑CGGCC‑GC‑TC‑TTTTTTTTGCA‑CGC‑ATAAAAACGAAAGAAGAAACccg < 3:217761/135‑1 (MQ=255) cATGGGTTTTCAATCACGCTGAAAAGATGACAAAATCCTAAAATGAAAACCATTTTAGGATTTTGTCAA‑CAATCTGAAGCGGCGGGTAAT‑A‑ACGCC‑GC‑GC‑TTTTTTTTTAA‑CGCAAAAAAAACAAAAAAAAAAAccccc > 1:99673/1‑137 (MQ=17) ttCAATCACGCTGAAAAGATGACAAAATCCTAAAATGAAAACCATTTTAGGATTTTGTCAA‑CAATCTGAAGCGGCGGTTTAT‑A‑ACGCC‑GC‑GCGTTTTTTTTGC‑‑CGC‑CTAAAAAACTAAGAATAAACCCCCCCCATAGta > 7:31661/1‑138 (MQ=17) ttCAATCACGCTGAAAAGATGACAAAATCCTAAAATGAAAACCATTTTAGGATTTTGTCAA‑CAATCTGAAGCGGCGGTTAAA‑A‑CCGCC‑GC‑TC‑TTTTTTTTTCA‑CCC‑ATAAAAACGAAAGAAGAAACCCCCCCCaaaatc > 7:125160/1‑134 (MQ=21) aaaTCACGATGAAAAGATCACAAAATCATAAAATTAAAACAATTTTAGTATTTTGTAAA‑CAATCTGAAGCGGCGGTTAAT‑A‑CCGCC‑GC‑TC‑TTTTTTTTGCA‑CGC‑ATAAAAACGAAAGAAGAAACCCGCCGCATAGT‑Cca < 3:269840/138‑1 (MQ=255) cGCTGAAAAGATGACAAAATCCTAAAATGAAAACCATTTTAGGATTTTGTCAA‑CAATCTGAAGCGGCGGTTAAT‑A‑CCGCC‑GCGGG‑TTTTTTTTGC‑‑CGG‑AAAAAAAACAAATAAGAAACCCCCCCCATAGTGCcaaaaa > 5:144271/1‑136 (MQ=17) aCAAAATCCAAAAATAAAAACCATTTTAGGATTTTGTAAA‑CAATCTGAAGCGGCGGTAAAT‑C‑CCGCC‑GC‑TC‑TTTTTTTTGCA‑CGC‑ATAAAAACGAAAGAAGAACCCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGt < 6:144271/139‑1 (MQ=38) aacaaaaCCAATTTATCCTTTTGTCAACCAACCTTCCCCGGCGGAAAAA‑C‑CCGCC‑GC‑TC‑TTTTTTTACCA‑CGC‑AAAAAACCCAACCACCACCCCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGt < 4:285396/136‑1 (MQ=255) caaa‑caaTCTGACGCGGC‑GTTAAA‑A‑CCCCCGGA‑TT‑TTTTATTACCA‑CCC‑AAAAAACCCAACCACCAACCCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGTTGATAAACGACATGAGAACGGCa < 3:172733/136‑1 (MQ=255) caaa‑caaTCGCACGCGGCGTTAAAA‑C‑CCGCG‑GG‑TT‑TTTTTTTCCCA‑CGA‑ATAAAAACGAACCAACAACCCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGTTGATAAACGACATGAGAACGGCa < 3:166280/136‑1 (MQ=255) acaa‑caaTCTGACGCGGCTGTAAAA‑C‑CCGCG‑GC‑TC‑TTTTTTTTGCC‑CCC‑ATAAAAACCAACCAACAACCCCGCCGCATAGT‑CCAAACAAATGGCCCATCATGTTAATATTGGAGTTGATAAACGACATGAGAACGGCa < 2:33858/137‑1 (MQ=255) tttacgcggcgGTAAAT‑A‑CGGCG‑GA‑TA‑TTTTTTTTGCA‑CGC‑ATAAAAACCAAACACGAACCCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGTTTATAAACGACATGAGAACGGCAAATGCGAg < 4:420899/133‑1 (MQ=255) acgcggcgGTTAAT‑A‑CCGCC‑GC‑TC‑TTTTTTTTCCA‑CGC‑AAAAAAACGAACGAAGAAACCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGTTGATAAACGACATGAGAACGGCAAATGCGAg < 5:174533/133‑1 (MQ=255) atccggggTTAAA‑ACCCGCG‑GC‑TA‑TTTTTTTAGCA‑CCA‑AAAAAAACCAAACAACAAACCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGTTGATAAACGACATGAGAACGGCAAATGCGAg < 3:384285/129‑1 (MQ=255) aaacccgcc‑tC‑TT‑TTTTTTTTCCA‑CCA‑AAAAAAACCAAACACTAACCCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGTTGATAAACGACATGAGAACGGCAAATGCGAGCAGTGTGATGATCa < 1:329103/134‑1 (MQ=255) cc‑gc‑TG‑GTTTTTTTGGA‑CGC‑ATAAAAACGAAAGAAGAAACCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTTATATTTGAGTTGTTAAACGACATGAGAACGGCAAATGCGAGCAGTGTGATGATCAattttgca > 8:409492/1‑135 (MQ=255) cc‑gc‑TA‑TTTTTTTTGCC‑CCA‑AAAAAAACCAAACAAGAACCCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGTTGATAAACGACATGAGAACGGCAAATGCGAGCAGTGTGATGATCATTTTGGaa < 6:187630/139‑1 (MQ=255) ttttCCA‑CGA‑AAAAAAACCAACCAACACACCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGTTGATAAACGACATGAGAACGGCAAATGCGAGCAGTGTGATGATCATTTTGGAATGCTCCTGc < 5:290905/138‑1 (MQ=255) ttGCA‑CGC‑ATAAAAACGAAAGAAGAAACCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGTTGATAAACGACATGAGAACGGCAAATGCGAGCAGTGTGATGATCATTTTGGAATGCTCCTGccc > 7:326060/1‑138 (MQ=255) | GTCCTTCCTAGACTTCATGGGTTTTCAATCACGCTGAAAAGATGACAAAATCCTAAAATGAAAACCATTTTAGGATTTTGTCAA‑CAATCTGAAGCGGCGGTTAAT‑A‑CCGCC‑GC‑TC‑TTTTTTTTGCA‑CGC‑ATAAAAACGAAAGAAGAAACCCGCCGCATAGT‑CCAAACAAATGGGCCATCATGTTAATATTGGAGTTGATAAACGACATGAGAACGGCAAATGCGAGCAGTGTGATGATCATTTTGGAATGCTCCTGCCC > BS168_reference_genome/514916‑515174 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |