Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 1,427,990 | G→A | 43.6% | intergenic (+296/‑178) | ydaY → / → ynaA | hypothetical protein/predicted tail protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 1,427,990 | 0 | G | A | 43.6% | 15.1 / 46.9 | 39 | intergenic (+296/‑178) | ydaY/ynaA | hypothetical protein/predicted tail protein |
Reads supporting (aligned to +/- strand): ref base G (9/13); new base A (11/6); total (20/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.00e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01 |
GATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGAT > W3110S.gb/1427930‑1428059 | gattgaGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAgtgtgtg < 1:1238708/70‑1 (MQ=255) ttgaGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAg > 1:1031051/1‑71 (MQ=255) ttgaGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAg > 1:767582/1‑71 (MQ=255) ttgaGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAg > 1:1815883/1‑71 (MQ=255) ttgaGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAg > 1:1298056/1‑71 (MQ=255) cgcgATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAgtgtgtg < 1:1775133/63‑1 (MQ=255) cgcgATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAgtgtgtg < 1:741045/63‑1 (MQ=255) cgcgATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAgtgtgtg < 1:611062/63‑1 (MQ=255) cgcgATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAgtgtgtg < 1:1878008/63‑1 (MQ=255) gATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAAt < 1:719464/71‑1 (MQ=255) aTGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAAt < 1:152910/70‑1 (MQ=255) aTCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa > 1:1317945/1‑69 (MQ=255) aTCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa > 1:922812/1‑69 (MQ=255) aTCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa > 1:1100802/1‑69 (MQ=255) aTCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa > 1:759047/1‑69 (MQ=255) aTCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa > 1:1967175/1‑69 (MQ=255) aTCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa > 1:196071/1‑69 (MQ=255) aTCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa > 1:1301047/1‑69 (MQ=255) aTCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa > 1:1383226/1‑69 (MQ=255) aTCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa > 1:1521280/1‑69 (MQ=255) aTCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa > 1:1622458/1‑69 (MQ=255) aTCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAAATACGTACACTCTGTTa > 1:1270123/1‑69 (MQ=255) cTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACa > 1:967614/1‑69 (MQ=255) cTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACa > 1:1651673/1‑69 (MQ=255) gCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAAt < 1:391709/47‑1 (MQ=255) cACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGAtgt < 1:1713816/71‑1 (MQ=255) cACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACGGAGAtgt < 1:633867/71‑1 (MQ=255) cACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGAtgt < 1:1835214/71‑1 (MQ=255) ggATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTacac > 1:371707/1‑51 (MQ=255) ggATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTacac > 1:1265126/1‑51 (MQ=255) ggATGACCAGTATTTTAAGGAGCAGTGTGTGTAGTTGTTAATTACGTacac > 1:151193/1‑51 (MQ=255) taCCAGTATTTTCAGGAGCAGTGTGTGTAGTTGGTAAt < 1:701642/37‑1 (MQ=38) aGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCgat < 1:1903621/71‑1 (MQ=255) aGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCgat < 1:204970/71‑1 (MQ=255) aGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCgat < 1:1211732/71‑1 (MQ=255) aGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCgat < 1:1175391/71‑1 (MQ=255) aGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCgat < 1:1021478/71‑1 (MQ=255) aGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCgat < 1:929269/71‑1 (MQ=255) aGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCgat < 1:175363/71‑1 (MQ=255) | GATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGAT > W3110S.gb/1427930‑1428059 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |