| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 603,395 | C→G | 28.3% | intergenic (+60/+21) | pheP → / ← ybdG | phenylalanine transporter/mechanosensitive channel protein, miniconductance |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 603,395 | 0 | C | G | 28.3% | 81.5 / 40.1 | 46 | intergenic (+60/+21) | pheP/ybdG | phenylalanine transporter/mechanosensitive channel protein, miniconductance |
| Reads supporting (aligned to +/- strand): ref base C (16/17); new base G (5/8); total (21/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.44e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-01 | |||||||||||
ATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAAT > NC_000913/603249‑603532 | atggatgATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCg < 3:1137466/148‑1 (MQ=255)atggatgATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCg < 3:643027/148‑1 (MQ=255)atggatgATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCg > 4:900033/1‑148 (MQ=255) tggatgATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAGACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGa < 3:900033/148‑1 (MQ=255) tgatATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATgag < 3:246926/148‑1 (MQ=255) tgatATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATgag < 4:151776/148‑1 (MQ=255) atatGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGagag < 2:788926/148‑1 (MQ=255) cGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATgag > 1:328456/1‑131 (MQ=255) cGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATgag < 2:328456/131‑1 (MQ=255) aTCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGa < 3:933559/125‑1 (MQ=255) aTCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGa > 4:933559/1‑125 (MQ=255) aTTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg > 4:615111/1‑109 (MQ=255) aTTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg < 3:615111/109‑1 (MQ=255) cATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGagag > 2:982190/1‑94 (MQ=39) cATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGagag < 1:982190/94‑1 (MQ=255) cATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGagag < 1:267094/94‑1 (MQ=39) cATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGagag > 2:267094/1‑94 (MQ=39) cATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCg < 2:849802/87‑1 (MQ=255) cATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCg > 1:849802/1‑87 (MQ=255) ttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGagag < 3:941026/91‑1 (MQ=255) ttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGagag > 4:941026/1‑91 (MQ=255) cTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGa < 2:803897/77‑1 (MQ=255) cTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGa > 1:803897/1‑77 (MQ=255) gAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCg < 2:454213/68‑1 (MQ=255) gAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCg > 1:454213/1‑68 (MQ=255) aGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg < 3:408318/65‑1 (MQ=255) aGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg > 4:408318/1‑65 (MQ=255) ggCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg < 4:565491/64‑1 (MQ=255) ggCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg > 3:565491/1‑64 (MQ=255) ttCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCg > 4:847922/1‑57 (MQ=255) ttCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCg < 3:847922/57‑1 (MQ=255) ctctCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaa < 3:691746/44‑1 (MQ=255) ctctCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaa > 4:691746/1‑44 (MQ=255) ctctCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGa > 1:687867/1‑130 (MQ=255) ctctCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGa < 2:687867/130‑1 (MQ=255) tctcATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCAt > 1:150844/1‑70 (MQ=37) tctcATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCAt < 2:150844/70‑1 (MQ=37) ctcATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGAAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAAt > 1:800263/1‑148 (MQ=255) cATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGa > 1:446585/1‑80 (MQ=255) cATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGa < 2:446585/80‑1 (MQ=255) tCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAat > 2:558228/1‑130 (MQ=255) tCCATTCGAAGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCAGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAat < 1:558228/130‑1 (MQ=255) ccATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTg < 1:142962/66‑1 (MQ=255) ccATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTg > 2:142962/1‑66 (MQ=255) ccATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc < 2:150820/90‑1 (MQ=255) ccATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc > 1:150820/1‑90 (MQ=255) | ATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAAT > NC_000913/603249‑603532 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |