| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 3,735,429 | G→A | 31.4% | A323T (GCA→ACA) | bglH → | carbohydrate‑specific outer membrane porin, cryptic |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 3,735,429 | 0 | G | A | 31.4% | 46.4 / 26.0 | 35 | A323T (GCA→ACA) | bglH | carbohydrate‑specific outer membrane porin, cryptic |
| Reads supporting (aligned to +/- strand): ref base G (10/14); new base A (8/3); total (18/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.46e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AGATATCCCGTTATGGGATAAAGCCACCTTAATGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCGCATCGGAAAAAGATAATCAGGATAATAACGGGTATTATGACTGGAAAGATACCTGGATGTTTGGCACATC > W3110S.gb/3735362‑3735499 | aGATATCCCGTTATGGGATAAAGCCACCTTAATGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCGCAt > 1:1229807/1‑71 (MQ=255) aTGGGATAAAGCCACCTTAATGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCGCATCGGAAAAAGAt > 1:1058451/1‑70 (MQ=255) aTGGGATAAAGCCACCTTAATGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCACATCGGAAAAAGAt < 1:1424818/70‑1 (MQ=255) aTGGGATAAAGCCACCTTAATGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCACATCGGAAAAAGAt > 1:478234/1‑70 (MQ=255) tAAAGCCACCTTAATGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCGCATCGGAAAAAGATAATCAgg < 1:1145762/71‑1 (MQ=255) tAAAGCCACCTTAATGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCGCATCGGAAAAAGATAATCAg > 1:863963/1‑70 (MQ=255) tAAAGCCACCTTAATGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCACATCGGAAAAAGATAATCAg > 1:707899/1‑70 (MQ=255) ccACCTTAATGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCGCATCGGAAAAAGATAATCAGGataat < 1:592500/71‑1 (MQ=255) cTTAATGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCGCATCGGAAAAAGATAATCAGGataat < 1:222872/67‑1 (MQ=255) cTTAATGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCGCATCGGAAAAAGATAATCAGGATAATAACg < 1:585836/71‑1 (MQ=255) ttAATGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCGCATCGGAAAAAGATAATCa > 1:161671/1‑59 (MQ=255) ttAATGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCACATCGGaaaaa > 1:1343175/1‑51 (MQ=255) ttAATGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCACATCGGAAAAAGa > 1:78201/1‑53 (MQ=255) tAATGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCACATCGGAAAAAGATaa > 1:18780/1‑55 (MQ=255) tAATGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCACATCGGAAAAAGATAATCAGGATAAAAACggg > 1:1313042/1‑71 (MQ=255) aaTGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCGCATCGGAAAAAGATAATCAGGATAATAACgg < 1:1518545/69‑1 (MQ=255) aTGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCGCATCGGAAAAAGATAATCAGGATAATAACGGGta > 1:1340103/1‑71 (MQ=255) tGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCGCATCGGAAAAAGATAAt < 1:1632119/53‑1 (MQ=255) tGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCGCATCGGAAAAAGATAATCAGGATAATAACGGGtat > 1:1738975/1‑71 (MQ=255) gTAAGTGGTCGTTATGTCACGGCAAACGAAATCGCATCGGAAAAAGATAATCAGGATAATAACGGGtatta > 1:1044424/1‑71 (MQ=255) tAAGTGGTCGTTATGTCACGGCAAACGAAAGCACATCGGAAAAAGATAATCAGGATAATAACGGGtatta < 1:1497002/70‑1 (MQ=255) aGTGGTCGTTATGTCACGGCAAACGAAAGCGCATCGGAAAAAGATAATCAGGATAATAACGGGTATTATGa > 1:1474549/1‑71 (MQ=255) aGTGGTCGTTATGTCACGGCAAACGAAAGCACATCGGAAAAAGATAATCAGTATAATAACGGGTATTATg < 1:696087/70‑1 (MQ=255) gTGGTCGTTATGTCACGGCAAACGAAAGCGCATCGGAAAAAGATAATCAGGATAATAACGGGtattat < 1:464629/68‑1 (MQ=255) gTGGTCGTTATGTCACGGCAAACGAAAGCGCATCGGAAAAAGATAATCAGGATAATAACGGGtattat < 1:1173493/68‑1 (MQ=255) tGGTCGTTATGTCACGGCAAACGAAAGCGCATCGGAAAAAGATAATCAGGATAATAACGGGTATTATGACt > 1:533374/1‑71 (MQ=255) tCGTTATGTCACGGCAAACGAAAGCGCATCGGAAAAAg < 1:1382668/38‑1 (MQ=255) tCGTTATGTCACGGCAAACGAAAGCGCATCGGAAAAAGATAATCAGGATAATAACGGGtattat < 1:447853/64‑1 (MQ=255) ttATGTCACGGCAAACGAAAGCGCATCGGAAAAAGATAATCAGGATAATAACGGGtat < 1:1628228/58‑1 (MQ=255) aTGTCACGGCAAACGAAAGCGCATCGGAAAAAGATAATCAGGATAATAACGGGtat > 1:1531943/1‑56 (MQ=255) aTGTCACGGCAAACGAAAGCACATCGGAAAAAGATAATCAGGATAATAACGGGTATTATGACTGGAAAGa > 1:1613036/1‑70 (MQ=255) tCACGGCAAACGAAAGCGCATCGGAAAAAGATAATCAGGata < 1:778019/42‑1 (MQ=255) aCGGCAAACGAAAGCGCATCGGAAAAAGATAATCAGGATAATAACGGGTATTATGACTGGAAAGATACCTg < 1:866808/71‑1 (MQ=255) gCAAACGAAAGCGCATCGGAAAAAGATAATCAGGATAATAACGGGTATTATGACTGGAAAGATAcc < 1:1490182/66‑1 (MQ=255) aaCGAAAGCACATCGGAAAAAGATAATCAGGATAATAACGGGTATTATGACTGGAAAGATACCTGGATGtt > 1:1760369/1‑71 (MQ=255) gcATCGGAAAAAGATAATCAGGATAATAACGGGTATTATGACTGGAAAGATACCTGGATGTTTGGCACATc > 1:577002/1‑71 (MQ=255) | AGATATCCCGTTATGGGATAAAGCCACCTTAATGGTAAGTGGTCGTTATGTCACGGCAAACGAAAGCGCATCGGAAAAAGATAATCAGGATAATAACGGGTATTATGACTGGAAAGATACCTGGATGTTTGGCACATC > W3110S.gb/3735362‑3735499 |
| 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: ‑ |