| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 4,174,520 | T→G | 8.6% | Y119S (TAT→TCT) | gspO ← | bifunctional prepilin leader peptidase and methylase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 4,174,520 | 0 | T | G | 8.6% | 133.2 / 8.3 | 58 | Y119S (TAT→TCT) | gspO | bifunctional prepilin leader peptidase and methylase |
| Reads supporting (aligned to +/- strand): ref base T (32/21); new base G (2/3); total (34/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.40e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01 | |||||||||||
ACACCACAGCACCCCATAGCCCGCTACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGAGTTTGTCGGGCAGGAGCTGAGTGC > W3110S.gb/4174456‑4174589 | acacCACAGCACCCCATAGCCCGCTACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTgcg > 1:3814282/1‑71 (MQ=255) caccacAGCACCCCATAGCCCGCTACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTgcg > 1:551854/1‑70 (MQ=255) caccacAGCACCCCATAGCCCGCTACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTgcg > 1:1757594/1‑70 (MQ=255) cacaGCACCCCATAGCCCGCTACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCAtt < 1:3693120/71‑1 (MQ=255) cacaGCACCCCATAGCCCGCTACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCAtt < 1:3442005/71‑1 (MQ=255) acaGCACCCCATAGCCCGCTACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTa < 1:70696/71‑1 (MQ=255) cACCCCATAGCCCGCTACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTaaa > 1:1167914/1‑69 (MQ=255) aCCCCATAGCCCGCTACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCa > 1:3236108/1‑63 (MQ=255) ccccATAGCCCGCTACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATAcc < 1:3813911/71‑1 (MQ=255) cATAGCCCGCTACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATAccaag > 1:1595436/1‑71 (MQ=255) aGCCCGCTACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATa > 1:1249895/1‑49 (MQ=255) aGCCCGCTACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATa > 1:238429/1‑49 (MQ=255) gCCCGCTACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATAccaagcc > 1:1199974/1‑69 (MQ=255) gCCCGCTACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATAccaagcc > 1:229997/1‑69 (MQ=255) cGCTACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATACCaagccaagcc < 1:837840/71‑1 (MQ=255) gCTACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATa > 1:2762561/1‑58 (MQ=255) gCTACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATACCAagccaagcca < 1:18190/71‑1 (MQ=255) tACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATACcaagccaagc < 1:637202/67‑1 (MQ=255) tACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATACcaagccaagc < 1:2953333/67‑1 (MQ=255) aCCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTaaa > 1:3419430/1‑53 (MQ=255) aCCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTaaa > 1:332701/1‑53 (MQ=255) aCCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCagag > 1:2602275/1‑71 (MQ=255) ccGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATACCAagccaagcca > 1:3172861/1‑67 (MQ=255) gcCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCg > 1:2465035/1‑71 (MQ=255) ccGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCgg > 1:1294512/1‑71 (MQ=255) cGTAAACCGGATCTTGTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGt < 1:1807864/71‑1 (MQ=255) cGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGt < 1:3728195/71‑1 (MQ=255) cGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGt < 1:2791386/71‑1 (MQ=255) gTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTa < 1:1883502/71‑1 (MQ=255) tAAACCGCATCATGTAAATCAATTAATCCATACTGCGCAtt < 1:403189/41‑1 (MQ=255) tAAACCGCATCATGTAAATCAATTAATCCATACTGCGCAtt < 1:1670137/41‑1 (MQ=255) tAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGt > 1:3423734/1‑69 (MQ=255) aaCCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAgg > 1:1690828/1‑71 (MQ=255) aaCCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAgg > 1:2583342/1‑71 (MQ=255) aCCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGt < 1:1563409/71‑1 (MQ=255) aCCGCATCATGTAAATCAATTAATCCAGACTGCGCATTAAATAccaagc > 1:2124628/1‑49 (MQ=255) aCCGAATCATGTAAATCAATTAATCCAGACTGCGCATTAAATAccaagc > 1:1400364/1‑49 (MQ=255) ccGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGt < 1:1578342/70‑1 (MQ=255) cGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGt > 1:3382746/1‑69 (MQ=255) cGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGa < 1:3576913/71‑1 (MQ=255) atGTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGAGTTTGt > 1:2978177/1‑71 (MQ=255) atGTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGAGTTTGt > 1:2788308/1‑71 (MQ=255) atGTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGAGTTTGt > 1:1720163/1‑71 (MQ=255) gTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTa < 1:1990910/57‑1 (MQ=255) gTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGAGTTTg < 1:73400/68‑1 (MQ=255) aaTCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCaga > 1:1450736/1‑45 (MQ=255) aaTCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGAGTTTGTCggg < 1:769872/70‑1 (MQ=255) tCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGAGTTTGTCGGGCAg > 1:3078522/1‑71 (MQ=255) aTTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTg < 1:3163104/54‑1 (MQ=255) aaTCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAgg > 1:473802/1‑49 (MQ=255) aaTCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGAGTTTGTCGGGCa > 1:3530675/1‑64 (MQ=255) aTCCAGACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGAGTTTGTCGGGCAGGAGCTGa < 1:2587224/71‑1 (MQ=255) aTCCAGACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGAGTTTGTCGGGCAGGAGCTGa < 1:3486955/71‑1 (MQ=255) aTCCAGACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGAGTTTGTCGGGCAGGAGCTGa < 1:3145654/71‑1 (MQ=255) cATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGAGTTTGTCGGGCAGGAGCTGAGTg > 1:3387424/1‑71 (MQ=255) cATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGAGTTTGTCGGGCAGGAGCTGAGTg > 1:204712/1‑71 (MQ=255) cATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGAGTTTGTCGGGCAGGAGCTGAGTg > 1:3693815/1‑71 (MQ=255) cATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGAGTTTGTCGGGCAGGAGCTGAGTg > 1:1587517/1‑71 (MQ=255) tACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGAg < 1:3686610/47‑1 (MQ=255) tACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGAGTTTGTCGGGCAGGAGCt < 1:465460/64‑1 (MQ=255) tACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGAGTTTGTCGGGCAGGAGCTGa > 1:3775154/1‑66 (MQ=255) tACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGAGTTTGTCGGGCAGGAGCTGAGTGc > 1:3770639/1‑70 (MQ=255) | ACACCACAGCACCCCATAGCCCGCTACCGCGCCGTAAACCGCATCATGTAAATCAATTAATCCATACTGCGCATTAAATACCAAGCCAAGCCAGAGCAGCGGTAAGGTGAGTTTGTCGGGCAGGAGCTGAGTGC > W3110S.gb/4174456‑4174589 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |