Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 1,593,713 | G→A | 14.2% | G4S (GGT→AGT) | divIB → | cell‑division initiation protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 1,593,713 | 0 | G | A | 14.2% | 72.2 / 3.4 | 28 | G4S (GGT→AGT) | divIB | cell‑division initiation protein |
Reads supporting (aligned to +/- strand): ref base G (14/10); new base A (2/2); total (16/12) | |||||||||||
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 |
AGTTCTGACTGAAGCTGTTCATATGATATACTGTAAGCAAACGACAAACGGCATCATAGTATGCCGTTTGTTTTGGAATAGACAGACTTTTAACAGCTGTTTCATTTGAATGAGGTGAACAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACAAAAGGCAAACCGCCGTCTCATTTCATTTATTATGCTCTTTTTTATTATGGTGCTGATCATTGTGTACCTG > BS168_reference_genome/1593579‑1593850 | aGTTCTGACTGAAGCTGTTCATATGATATACTGTAAGCAAACGACAAACGGCATCATAGTATGCCGTTTGTTTTGGATTAGACAGACTTTTAACAGCAGTTTCATTTGAATGACGTGAACAGGCAATGAACCCGGGta > 6:221397/1‑137 (MQ=255) atGATATACTGTAAGCAAACGACAAACGGCATCATATTATGCCGTTTTTTTTGGAATAGACAGACTTTTAACAGCTGTTTCATTTGAATGAGGTGAACAGGCAATGAACCCGGGTAAAGGCCGACAAAAAATCGTTaa > 2:92358/1‑138 (MQ=255) atatACTGTAAGCAAACGACAAACGGCATCATAGTATGCCGTTTGTTTTGGAATAGACAGACTTTTAACAGCTGTTTCATTTGAATGAGGTGAACAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGTTAACAtt > 3:98353/1‑139 (MQ=255) aGCAAACGACAAACGGCATCATAGTATGCCGTTTGTTTTGGAATAGACAGACTTTTAACAGCTGTTTCATTTGAATGAGGTGAACAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGa > 7:102303/1‑139 (MQ=255) gCAAACGACAAACGGCATCATAGTATGCCGTTTGTTTTGGAATAGACAGACTTTTAACAGCTGTTTCATTTGAATGAGGTGAACAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGAt > 4:22476/1‑139 (MQ=255) cGTTTGTTTTGGAATAGACAGACTTTTAACAGCTGTTTCATTTGAATGAGGTGAACAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGt > 4:140075/1‑93 (MQ=255) cGTTTGTTTTGGAATAGACAGACTTTTAACAGCTGTTTCATTTGAATGAGGTGAACAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGt < 3:140075/93‑1 (MQ=255) tttgtttTGGAATAGACAGACTTTTAACAGCTGTTTCATTTGAATGAGGTGAACAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACaa > 3:29743/1‑137 (MQ=255) tttgtttTGGAATAGACAGACTTTTAACAGCTGTTTCATTTGAATGAGGTGAACAGGCAATGAACCCGAGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACaa > 2:146645/1‑137 (MQ=255) tttgtttTGGAATAGACAGACTTTTAACAGCTGTTTCATTTGAATGAGGTGAACAGGCAATGAACCCGAGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACaa > 2:143591/1‑137 (MQ=255) ttgtttTGGAATAGACAGACTTTTAACAGCTGTTTCATTTGAATGAGGTGAACAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGt > 6:162320/1‑90 (MQ=255) ttgtttTGGAATAGACAGACTTTTAACAGCTGTTTCATTTGAATGAGGTGAACAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGt < 5:162320/90‑1 (MQ=255) tttGGAATAGACAGACTTTTAACAGCTGTTTCATTTGAATGAGGTGAACAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACAAAAGGCa > 8:172908/1‑138 (MQ=255) tacTTTTAACAGCTGTTTCATTTGAATGAGGTGAACAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACGAAAGGCAAACCGCCGTCtcat < 1:92358/138‑1 (MQ=255) gacTTTTAACAGCTGTTTCATTTGAATGAGGTGAACAGGCAATGAACCCGAGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACAAAAGGCAAACCGCCGTCtcat < 1:146645/139‑1 (MQ=255) gacTTTTAACAGCTGTTTCATTTGAATGAGGTGAACAGGCAATGAACCCGAGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACAAAAGGCAAACCGCCGTCtcat < 1:143591/139‑1 (MQ=255) acTTTTAACAGCTGTTTCATTTGAATGAGGTGAACAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACAAAAGGCAAACCGCCGTCtcatt > 2:84018/1‑139 (MQ=255) gCTGTTTCATTTGAATGAGGTGAACAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACAAAAGGCAAACCGCCGTCTCAATTCATttatta > 3:173986/1‑139 (MQ=255) aTGAGGTGAACAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACAAAAGGCAAACCGCCGTCTCATTTCATTTATTATGCTCTTTTttatt > 8:162192/1‑139 (MQ=255) tGAGGTGAACAGGCAATGAACCCGGGTCCAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACAAAAGGCAAACCGCCGTCTCATTTCATTTATTATGCTCTTTTttatta < 5:54824/139‑1 (MQ=255) tGAGGTGAACAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACAAAAGGCAAACCGCCGTCTCATTTCATTTATTATGCTCTTTTttatta < 8:102303/139‑1 (MQ=255) tGAGGTGAACAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACAAAAGGCAAACCGCCGTCTCATTTCATTTATTATGCTCTTTTttatta < 1:84018/139‑1 (MQ=255) aGGTGAACAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACAAAAGGCAAACCGCCGTCTCATTTCATTTATTATGCTCTTTTTTATTATg < 4:29743/139‑1 (MQ=255) ttGAACTGGGAATGTTCCCGGGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGTTCCGTAATATTAAAGAGCAGAGGAAACAAAAGGCATACCGCCGTCTCATTTCATTTATTATGCTCTTTTTTATTATGGt < 4:173986/138‑1 (MQ=255) ggAACAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACAAAAGGCAAACCGCCGTCTCATTTCATTTATTATGCTCTTTTTTATTATGGTg < 3:88511/138‑1 (MQ=255) aaCAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACAAAAGGCAAACCGCCGTCtcatttcat < 1:123572/111‑1 (MQ=255) aaCAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACAAAAGGCAAACCGCCGTCtcatttcat > 2:123572/1‑111 (MQ=255) aaTGAACCCGGGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACAAAAGGCAAACCGCCGTCTCATTTCATTTATTATGCTCTTTTTTATTATGGTGCTGATCAtt > 8:67175/1‑139 (MQ=255) gggTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACAAAAGGCAAACCGCCGTCTCATTTCATTTATTATGCTCTTTTTTATTATGGTGCTGATCATTGTGTACCTg > 3:43787/1‑139 (MQ=255) | AGTTCTGACTGAAGCTGTTCATATGATATACTGTAAGCAAACGACAAACGGCATCATAGTATGCCGTTTGTTTTGGAATAGACAGACTTTTAACAGCTGTTTCATTTGAATGAGGTGAACAGGCAATGAACCCGGGTCAAGACCGAGAAAAAATCGTTAACATTGAAGAGCGGATCCCTAAAATTAAAGAGCAGAGGAAACAAAAGGCAAACCGCCGTCTCATTTCATTTATTATGCTCTTTTTTATTATGGTGCTGATCATTGTGTACCTG > BS168_reference_genome/1593579‑1593850 |
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: ‑ |