| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 4,002,709 | A→G | 11.1% | W31R (TGG→CGG) | yigF ← | DUF2628 family putative inner membrane protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,002,709 | 0 | A | G | 11.1% | 105.6 / 4.9 | 36 | W31R (TGG→CGG) | yigF | DUF2628 family putative inner membrane protein |
| Reads supporting (aligned to +/- strand): ref base A (16/16); new base G (2/2); total (18/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01 | |||||||||||
AAATATATCCAGGAACAGAAGAAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTT > NC_000913/4002603‑4002821 | aaataTATCCAGGAACAGAAGAAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCa > 4:442106/1‑114 (MQ=255)aaataTATCCAGGAACAGAAGAAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCa < 3:442106/114‑1 (MQ=255) atatatCCAGGAACAGAAGAAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCa < 1:820842/112‑1 (MQ=255) atatatCCAGGAACAGAAGAAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCa > 2:820842/1‑112 (MQ=255) ggAACAGAAGAAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATcc > 4:420930/1‑102 (MQ=255) ggAACAGAAGAAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATcc < 3:420930/102‑1 (MQ=255) cTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCAt > 2:927638/1‑99 (MQ=255) cTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCAt < 1:927638/99‑1 (MQ=255) tttAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATaaaaa > 4:979398/1‑75 (MQ=255) tttAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATaaaaa < 3:979398/75‑1 (MQ=255) tttAAATGCTGCTTTATATTCTGGTGTTGCCCCCCGGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCAt < 2:1080556/104‑1 (MQ=255) tttAAATGCTGCTTTATATTCTGGTGTTGCCCCCCGGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCAt > 1:1080556/1‑104 (MQ=255) tttAAATGCTGCTTTATATTCTGGTGTTGCCCCCCGGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCAt < 4:894266/104‑1 (MQ=255) tttAAATGCTGCTTTATATTCTGGTGTTGCCCCCCGGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCAt > 3:894266/1‑104 (MQ=255) tttAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCg > 4:146290/1‑106 (MQ=255) tttAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCg < 3:146290/106‑1 (MQ=255) aaaTGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGaa < 1:962086/97‑1 (MQ=255) aaaTGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATATTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGaa > 2:962086/1‑97 (MQ=255) aTGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCtt > 2:862950/1‑115 (MQ=255) aTGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCtt < 1:862950/115‑1 (MQ=255) ctgctTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATAtt > 3:1039548/1‑108 (MQ=255) ctgctTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATAtt < 4:1039548/108‑1 (MQ=255) tCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATaaaaa > 2:129700/1‑48 (MQ=255) tCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATaaaaa < 1:129700/48‑1 (MQ=255) tCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGttt < 1:419313/129‑1 (MQ=255) tCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGttt > 2:419313/1‑129 (MQ=255) tGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGtt > 1:214683/1‑120 (MQ=255) tGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGtt < 2:214683/120‑1 (MQ=255) cccccAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCAt < 1:140692/96‑1 (MQ=255) cccccAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCAt > 2:140692/1‑96 (MQ=255) cccAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATa < 1:649311/106‑1 (MQ=255) cccAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATa > 2:649311/1‑106 (MQ=255) ccAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACt < 3:965377/100‑1 (MQ=255) ccAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACt > 4:965377/1‑100 (MQ=255) ccAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATa < 3:851290/105‑1 (MQ=255) ccAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATa > 4:851290/1‑105 (MQ=255) | AAATATATCCAGGAACAGAAGAAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTT > NC_000913/4002603‑4002821 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |