Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 1,803,104 | A→G | 20.6% | pseudogene (11/426 nt) | arpB → | pseudogene, ankyrin repeats |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,803,104 | 0 | A | G | 20.6% | 163.7 / 41.8 | 85 | pseudogene (11/426 nt) | arpB | pseudogene, ankyrin repeats |
Reads supporting (aligned to +/- strand): ref base A (27/39); new base G (6/11); total (34/51) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.85e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01 |
ATAACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCG > NC_000913/1803006‑1803201 | ataACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAAc < 1:1296527/100‑1 (MQ=255) aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGat < 1:1555215/100‑1 (MQ=255) aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGat > 2:88981/1‑100 (MQ=255) aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGat > 1:235545/1‑100 (MQ=255) aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGat < 1:2304960/100‑1 (MQ=255) aaTTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGatat < 1:1213012/100‑1 (MQ=255) atgtatCGTGCATATTTCGTTGTGGCGTTAGTGATTAATCATCATCCTCAAAACCTCCATATATATATAAGGCTGTCGGATGATGAGTCAAAACGatata < 2:875426/95‑2 (MQ=255) aTTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATc < 2:306125/100‑1 (MQ=255) aTTTTTCCTGCACATTTCGTTATGGCGCTAATAATCAATCATCATCCTCAAAGCCTCCAGATATATATAAGGCTAACAGATTATGAGTCAAACCGATATc < 2:2106365/100‑1 (MQ=255) tttttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCa < 2:165432/100‑1 (MQ=255) tttttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCa < 2:165510/100‑1 (MQ=255) tttttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCa < 2:165290/100‑1 (MQ=255) tttttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCa < 2:1542352/100‑1 (MQ=255) ttttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCAt < 1:88981/100‑1 (MQ=255) tttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCAtt < 1:987017/100‑1 (MQ=255) tCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTAt < 2:900750/100‑1 (MQ=255) gTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCa < 2:423029/100‑1 (MQ=255) tGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAg < 2:2234025/100‑1 (MQ=255) gCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGa < 2:64905/100‑1 (MQ=255) acaTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAAc < 1:2091795/100‑1 (MQ=255) acaTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAAc < 1:115860/100‑1 (MQ=255) caTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACt > 2:841486/1‑100 (MQ=255) tttCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCa > 2:1520941/1‑100 (MQ=255) tttCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCa > 2:2017997/1‑100 (MQ=255) tttCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCa > 2:871904/1‑100 (MQ=255) cGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATta < 1:572272/100‑1 (MQ=255) gTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATtat > 1:2274738/1‑100 (MQ=255) ttCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATAGATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATtata < 1:1626622/100‑1 (MQ=255) tgtgtcGATAATAATTAATCATCAAACTAAAAACCTCTAGATATATATAAGGCTAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATtata < 2:2227386/94‑1 (MQ=255) cTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAg < 2:2037772/100‑1 (MQ=255) tGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTAATTATAAGt > 2:1181279/1‑100 (MQ=255) ggCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGtc < 2:2153712/100‑1 (MQ=255) gCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGtct > 2:1082539/1‑100 (MQ=255) ataatTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCAt > 2:1798456/1‑100 (MQ=255) taattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATa > 1:1954833/1‑100 (MQ=255) aattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAg > 2:1706440/1‑100 (MQ=255) aattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATCATAAGTCTCCTCATAg > 2:1376739/1‑100 (MQ=255) taatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGAtt < 2:1035263/100‑1 (MQ=255) taatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGAt < 1:300080/99‑1 (MQ=255) tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca > 1:1076520/1‑100 (MQ=255) tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca < 2:788562/100‑1 (MQ=255) tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca < 1:1742261/100‑1 (MQ=255) catcatcCTCAAAACCTCCATATATATATATAGCGTACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcac < 2:1705142/100‑1 (MQ=255) catcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcac > 2:273868/1‑100 (MQ=255) tcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCACCTCATAGATTGCACAt > 1:216724/1‑100 (MQ=255) catcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATc < 2:2106461/100‑1 (MQ=255) ccTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGAt > 2:1936143/1‑100 (MQ=255) ccTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGAt > 2:152521/1‑100 (MQ=255) tCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAg > 2:709275/1‑100 (MQ=255) tCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAg > 1:576171/1‑100 (MQ=255) aaaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCg < 2:389354/100‑1 (MQ=255) aaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGa < 1:2091299/100‑1 (MQ=255) aaaCCACCCGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTACCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGa < 1:106060/100‑1 (MQ=255) aCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACa < 1:1099590/100‑1 (MQ=255) ccTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACAt > 1:903779/1‑100 (MQ=255) ctacAGATATATATTAGTCGACCAGATTATGTGTCAAAACGATATCATTATTAAAATTCATTATAAGTCTCATCATATATTGCACATCGATAGCTTCATa > 2:1005904/4‑100 (MQ=255) ccAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAcc < 1:1462929/100‑1 (MQ=255) ccAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAcc < 2:2076364/100‑1 (MQ=255) ccAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTAATAGATTGCACATCGATAGCGACATAcc > 2:1448867/1‑100 (MQ=255) cAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGTACATCGATAGCGACATACCa < 2:1442188/100‑1 (MQ=255) gATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac < 2:1462087/100‑1 (MQ=255) gATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac < 1:1171083/100‑1 (MQ=255) gATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac < 1:711130/100‑1 (MQ=255) atatatatAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAaca < 1:1739516/100‑1 (MQ=255) atatatatAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAaca < 1:2327910/100‑1 (MQ=255) tatatataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTACTCATAGATTGCACATCGATAGCGACATACCAaca > 2:1090640/1‑99 (MQ=255) atataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCtt > 1:1599310/1‑100 (MQ=255) tataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTGTAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtc > 2:1343199/1‑100 (MQ=255) ataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtca > 2:1139268/1‑100 (MQ=255) taAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcat > 1:267392/1‑100 (MQ=255) aaGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATGGATTGCACATCGATAGCGACATACCAACACCTtcatc < 2:1255638/100‑1 (MQ=255) aGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcatca < 2:332503/100‑1 (MQ=255) ggCGAACAGATTATGAGTCAAAACGATATCATTATCCGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCag > 1:1656635/1‑100 (MQ=255) gCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCaga < 1:1490542/100‑1 (MQ=255) gAACAGATTATGAGTCAAATCGATATTATTATCTGAACTCATTATATGTCGCTTCTTATATTGCACATCGATATCTATTTACTTTCTCCTTCATCAGAGc > 2:783317/1‑100 (MQ=255) aaCAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGcc < 1:1761252/100‑1 (MQ=255) cAGATTACGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTa > 1:1952529/1‑100 (MQ=255) aGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAt < 1:1192074/100‑1 (MQ=255) gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt < 1:2241611/100‑1 (MQ=255) tATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAAt < 2:576168/100‑1 (MQ=255) tGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTACTCATAGATTGCACATCGATAGCGACATACCAACACATTCATCAGAGCCTATTAATCa > 2:1145162/1‑100 (MQ=255) aGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCa < 2:1936201/98‑1 (MQ=255) gTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAtt > 2:2194977/1‑100 (MQ=255) tCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAttt > 1:17719/1‑100 (MQ=255) cAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg < 1:2145389/100‑1 (MQ=255) aaaCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACCTCGATAGCGACATACCAACACCTTCATCAGAGACTATTAATCAATTTgcg > 1:889637/1‑100 (MQ=255) aaaCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgcg < 2:2268460/100‑1 (MQ=255) aaaCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgcg > 1:393927/1‑100 (MQ=255) aaaCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgcg > 1:858259/1‑100 (MQ=255) | ATAACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCG > NC_000913/1803006‑1803201 |
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: ‑ |