| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | AP009048 | 3,294,438 | T→C | 10.6% | L395S (TTG→TCG) | yraM → | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | AP009048 | 3,294,438 | 0 | T | C | 10.6% | 175.7 / 12.5 | 66 | L395S (TTG→TCG) | yraM | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base T (27/32); new base C (4/3); total (31/35) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.98e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01 | |||||||||||
TAAGCCAGGTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAATATCGAAAATCGCG > AP009048/3294369‑3294503 | tAAGCCAGGTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTg > 1:1572546/1‑71 (MQ=255)tAAGCCAGGTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTg > 1:333228/1‑71 (MQ=255) cAGGTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGaa < 1:943094/71‑1 (MQ=255) aGGTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAag < 1:77216/71‑1 (MQ=255) ggTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGaa > 1:2329198/1‑69 (MQ=255) gTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGtt < 1:1810743/62‑1 (MQ=255) ttCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGc < 1:1702360/71‑1 (MQ=255) ttCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTCGCTGAAGAGc < 1:1417300/71‑1 (MQ=255) tCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGc > 1:4340069/1‑62 (MQ=255) tCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGc > 1:4334840/1‑62 (MQ=255) caggagGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAaga < 1:2298645/67‑1 (MQ=255) agcagGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGaa < 1:3118853/64‑1 (MQ=255) gATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTcc < 1:1616955/71‑1 (MQ=255) gcgcGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGc < 1:2163795/59‑1 (MQ=255) gAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGATCg < 1:2032248/71‑1 (MQ=255) aGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTCGCTGAAGAGCAACACTCCGCTGAACGt < 1:2429191/71‑1 (MQ=255) gTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTa < 1:166508/71‑1 (MQ=255) tATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGtt < 1:2470767/42‑1 (MQ=255) ttGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCt > 1:2868854/1‑43 (MQ=255) ttGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACt < 1:1976112/70‑1 (MQ=255) ttGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTg < 1:3341617/71‑1 (MQ=255) tgtgGTCGGTCTGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAagag > 1:119606/1‑48 (MQ=255) gtgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTgg > 1:3293590/1‑70 (MQ=255) gtgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTgg < 1:2640447/70‑1 (MQ=255) tgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGc < 1:2451477/70‑1 (MQ=255) tgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGc < 1:4522395/70‑1 (MQ=255) tgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGc < 1:2598185/70‑1 (MQ=255) tcggtcCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGc < 1:3463361/67‑1 (MQ=255) cggtcCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACt > 1:1726086/1‑69 (MQ=255) ggtcCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGa < 1:4176299/70‑1 (MQ=255) ggtcCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGa > 1:27152/1‑70 (MQ=255) gtcCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCt < 1:146362/52‑1 (MQ=255) tcCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAAc > 1:1006567/1‑70 (MQ=255) tcCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAAc > 1:4304310/1‑70 (MQ=255) cGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTg > 1:924032/1‑59 (MQ=255) cGTTGCTGAAAAATAACGTTGAAGAGTCGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAg > 1:2906663/1‑71 (MQ=255) gTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGa < 1:1730537/65‑1 (MQ=255) gTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGc < 1:2169685/71‑1 (MQ=255) ttGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGcc > 1:1955400/1‑71 (MQ=255) ttGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGcc < 1:288238/71‑1 (MQ=255) ttGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGcc > 1:2023554/1‑71 (MQ=255) ttGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGcc < 1:1332797/71‑1 (MQ=255) ttGCTGAAAAATAACGTTGAAGAGTCGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAAc > 1:971584/1‑66 (MQ=255) ttGCTGAAAAATAACGTTGAAGAGTCGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGcc > 1:976957/1‑71 (MQ=255) gCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCg > 1:3555336/1‑70 (MQ=255) tGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTgg > 1:3244128/1‑54 (MQ=255) tGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCgg < 1:2460516/69‑1 (MQ=255) gAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGgag > 1:1687017/1‑70 (MQ=255) gAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGgag > 1:1912328/1‑70 (MQ=255) aaaaaTAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACt > 1:1540741/1‑50 (MQ=255) aaaaaTAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGgag < 1:1713540/69‑1 (MQ=255) aaaaaTAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGaa > 1:4314521/1‑71 (MQ=255) aaaaTAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCa < 1:1394849/61‑1 (MQ=255) aaaaTAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAat < 1:2269904/71‑1 (MQ=255) aaaTAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGt < 1:3308640/45‑1 (MQ=255) aaaTAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGaa > 1:836327/1‑57 (MQ=255) aaaTAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAata > 1:3387809/1‑71 (MQ=255) aTAACGTTGAAGAGTCGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAATATc > 1:4608042/1‑71 (MQ=255) tAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACg < 1:2069577/41‑1 (MQ=255) aaCGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGgaga > 1:522272/1‑64 (MQ=255) aCGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGaa > 1:4137259/1‑37 (MQ=255) aCGTTGAAGAGTCGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAATATCGaa < 1:4540947/71‑1 (MQ=255) gTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAATATCGaaaa < 1:1872369/71‑1 (MQ=255) gTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAATATCGaaaa < 1:4421539/71‑1 (MQ=255) gAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGgag < 1:4318029/57‑1 (MQ=255) gAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAATATCGAAAATc < 1:4328613/70‑1 (MQ=255) agagTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGgaga > 1:3257322/1‑56 (MQ=255) agagTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAATATCGAAAATcgcg > 1:2137646/1‑71 (MQ=255) | TAAGCCAGGTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAATATCGAAAATCGCG > AP009048/3294369‑3294503 |
| 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: ‑ |