| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 1,325,850 | C→T | 48.4% | R369C (CGC→TGC) | yciQ → | enhancer of membrane protein expression; putative inner membrane protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,325,850 | 0 | C | T | 48.4% | 5.1 / 80.6 | 62 | R369C (CGC→TGC) | yciQ | enhancer of membrane protein expression; putative inner membrane protein |
| Reads supporting (aligned to +/- strand): ref base C (19/13); new base T (15/15); total (34/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.10e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01 | |||||||||||
GAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGG > NC_000913/1325751‑1325949 | gAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGAc < 2:412762/100‑1 (MQ=255) aGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGc < 1:1759056/100‑1 (MQ=255) aGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGc < 1:1379164/100‑1 (MQ=255) aGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGc > 1:1819000/1‑100 (MQ=255) ggCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGc < 2:1903972/99‑1 (MQ=255) cgACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACag < 2:1818979/100‑1 (MQ=255) aCCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACagag < 1:1921951/100‑1 (MQ=255) aCCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACagag < 1:1663228/100‑1 (MQ=255) ccAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGt > 2:2359320/1‑100 (MQ=255) taatGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGc > 1:2252574/1‑100 (MQ=255) aatGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCTTTCTTACGGTGAGACGCAACAGAGTCTTACGc < 2:1342966/99‑1 (MQ=255) atGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGAAg > 1:2057326/1‑100 (MQ=255) gCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAggg < 1:52579/100‑1 (MQ=255) cGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAgg > 2:1091704/1‑98 (MQ=255) cGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGc > 1:1637558/1‑100 (MQ=255) tATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCAtt > 1:392077/1‑100 (MQ=255) aCTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTa > 1:1956668/1‑100 (MQ=255) cTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTAc > 1:1581442/1‑100 (MQ=255) gCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACgg < 1:1500774/100‑1 (MQ=255) gCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACgg > 2:2363661/1‑100 (MQ=255) gCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACgg > 2:1209674/1‑100 (MQ=255) gCCCGCAGAGCAGTTATTACTGGAAGCGGTATTTCGTAAAGGCGCTAAGGTCGTTCTTACGGTGAGATGCAACAGAGTCTTATGCAGGGCATTTTTACg > 1:34559/1‑99 (MQ=255) cGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGAt > 1:1387681/1‑99 (MQ=255) cAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCa < 2:1590395/100‑1 (MQ=255) gagCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGa > 2:467007/1‑100 (MQ=255) gagCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGa < 2:2099197/100‑1 (MQ=255) gagCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGa > 1:2308523/1‑100 (MQ=255) gCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGaaa < 2:531001/100‑1 (MQ=255) gCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGaaa < 2:388157/100‑1 (MQ=255) gTTATTACTGGGAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAAtt < 2:1423778/99‑1 (MQ=255) tattaCTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGGTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGTCATTTTTACGGATGCTGAAATTTTa > 1:265721/1‑100 (MQ=255) cTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGc > 1:1373958/1‑100 (MQ=255) cTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGc < 2:1871403/100‑1 (MQ=255) ggAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCcg > 1:1154946/1‑100 (MQ=255) ggAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGcc > 2:51240/1‑99 (MQ=255) ggAAGCGTTATTTCGTAAAGGCGATAAAGTCGTTCTTACGGGGAGACGCAACAGAGGCTTACGCAGGGCATCTTTACGGATGCAGAAATTTTATCAGCcg > 2:775499/1‑100 (MQ=255) gAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCcgc > 2:342903/1‑100 (MQ=255) gAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCcgc < 2:1740888/100‑1 (MQ=255) gAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCcg < 1:528952/99‑1 (MQ=255) ncGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTTTGTATGCTGTAATTTTATCTGCCGCGt > 2:1707152/2‑99 (MQ=255) gCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTa < 1:342906/100‑1 (MQ=255) gTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTaaa > 2:210044/1‑100 (MQ=255) gTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTaaa > 2:935423/1‑100 (MQ=255) gTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTaaa > 1:1422731/1‑100 (MQ=255) cGTAAAGGCGATAAGGTCGTTATTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTAACTGAAGAATAAAAAGTCtt < 2:2325174/100‑1 (MQ=255) aaaGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGt > 1:768114/1‑100 (MQ=255) gATAAGATCGTTCTTACTGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTAcn < 2:1498242/100‑2 (MQ=255) tAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCgg > 2:1937001/1‑100 (MQ=255) tAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCgg > 2:388702/1‑100 (MQ=255) aGGTCGTTCTTACTGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTc < 1:1469481/100‑1 (MQ=255) aGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTc > 1:1660599/1‑100 (MQ=255) aGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTc < 1:1909213/100‑1 (MQ=255) ggTCGTTCTTACGGGGAGACGGAACAGGGTCTAATGGAGGGCAGTGATACGGATGTAGAAATTTTATCTGCTGCGTAAAAAGttttgggtgtag < 2:23055/94‑9 (MQ=255) cTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCCGATACtt < 1:1258857/100‑1 (MQ=255) tACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACtttt > 1:760415/1‑100 (MQ=255) tACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACtttt > 2:1790631/1‑100 (MQ=255) ggggAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTg < 1:280129/100‑1 (MQ=255) ggggAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTg < 2:58231/100‑1 (MQ=255) ggggAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTg < 2:149078/100‑1 (MQ=255) gggAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGc < 1:1316657/100‑1 (MQ=255) gggAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGc > 1:2032707/1‑100 (MQ=255) gaCGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATg > 1:1839157/1‑100 (MQ=255) gaCGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATg > 2:2151659/1‑100 (MQ=255) cGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATggg > 2:1416984/1‑100 (MQ=255) | GAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGCAATGGG > NC_000913/1325751‑1325949 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |