| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 1,204,942 | T→C | 16.1% | pseudogene (539/1134 nt) | aaaE → | pseudogene, phage terminase protein A family, e14 prophage;Phage or Prophage Related |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,204,942 | 0 | T | C | 16.1% | 163.2 / 22.8 | 74 | pseudogene (539/1134 nt) | aaaE | pseudogene, phage terminase protein A family, e14 prophage;Phage or Prophage Related |
| Reads supporting (aligned to +/- strand): ref base T (37/25); new base C (5/7); total (42/32) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.42e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01 | |||||||||||
GCGAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCGAGAAAAATGAAG > NC_000913/1204843‑1205040 | gCGAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACt > 1:2050052/1‑100 (MQ=255) cGAGGATAAATCCCTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACtt > 1:1743276/1‑100 (MQ=255) cGAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACtt < 2:489300/100‑1 (MQ=255) cGAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACtt < 1:982676/100‑1 (MQ=255) gATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGAACTGGCGCGTAAGCTGGAAATGAAAAGCATAGCGCGACTTtata > 1:604007/1‑100 (MQ=255) aTAAATCACTGACCCTTGAGCAGTTCGAGGTGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGATGAATATAAACAGCATGGCGCGACCATATAc > 1:295326/1‑100 (MQ=255) tAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATAcc < 1:763735/100‑1 (MQ=255) aaaTCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATAccc > 2:1707809/1‑100 (MQ=255) aaTCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCcg > 1:1494471/1‑100 (MQ=255) cACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCga < 1:850408/100‑1 (MQ=255) aCTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCgag > 1:1191586/1‑100 (MQ=255) aCTGACCCTGGAGCAGTTCTAGGGGCGTCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCgag < 2:2004316/100‑1 (MQ=255) cTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGATCAGCATGGCGCGACTTTATACCCGCgaga < 2:448341/100‑1 (MQ=255) tGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGAt > 2:1606654/1‑100 (MQ=255) cccTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTg < 2:2169815/99‑1 (MQ=255) cccTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGa > 1:788183/1‑100 (MQ=255) ccTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGAc < 2:1943416/100‑1 (MQ=255) ttGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACgg > 2:116555/1‑100 (MQ=255) gAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTa > 2:1301635/1‑100 (MQ=255) gAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTa > 2:2545591/1‑100 (MQ=255) gCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTaaa > 1:737301/1‑100 (MQ=255) cAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTaaaa > 2:675177/1‑100 (MQ=255) aGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAAc < 1:2293758/100‑1 (MQ=255) aGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAAc < 1:2184903/100‑1 (MQ=255) aGTTCGAGGGGAAGCAGTGCATTCTGTCCTTTGATCTGGTGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTaaa < 1:2566357/98‑1 (MQ=255) gTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACg < 1:2469690/100‑1 (MQ=255) ttCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGc < 2:1208330/100‑1 (MQ=255) tCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCa > 2:1134497/1‑100 (MQ=255) cGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCAt > 2:685912/1‑100 (MQ=255) gAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCAtt > 2:1375397/1‑100 (MQ=255) gAGGGGCAGCCGTGCATTCTGGCCTTTGACCTCGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCAtt > 1:2519513/1‑100 (MQ=255) aGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATta < 2:406830/100‑1 (MQ=255) aGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATta < 1:1467977/100‑1 (MQ=255) ggggCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATtac < 2:1087987/100‑1 (MQ=255) ggggCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATtac < 2:2342600/100‑1 (MQ=255) gggAAGCCGTGCATTCTGGCCTTTGACCTGGCGCGGAAGCGGGATATTAATAGCATGGCGCGACCTTATACGCGCGAGATTGACGGTAAAGCGAATtac < 2:885152/99‑1 (MQ=255) gtcAGCCGTTCATTCTGGCCTTTGTCCTTGCTCGTAAGCTTGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATtac < 1:2099605/96‑1 (MQ=255) gCAGCCGTGCATTTTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATtactac < 1:2633173/100‑1 (MQ=255) cAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATAACCGAGAGATTAAAGGTAAAACGCATTACTACa > 1:404713/1‑100 (MQ=255) gCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAgt > 2:727826/1‑100 (MQ=255) gCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAgt > 2:424445/1‑100 (MQ=255) cGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAgtgt > 2:2432561/1‑100 (MQ=255) cGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAgtgt > 1:2211213/1‑100 (MQ=255) tGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTgg < 2:2339916/100‑1 (MQ=255) aTTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGcc > 2:2300816/1‑99 (MQ=255) ggCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCcgcg > 2:241489/1‑99 (MQ=255) ggCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGt < 1:1325609/100‑1 (MQ=255) ccTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGttt > 2:2005252/1‑100 (MQ=255) cTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTc > 2:1210464/1‑100 (MQ=255) ttGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTg < 2:561890/100‑1 (MQ=255) ttGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTg < 1:865747/100‑1 (MQ=255) gACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTggg < 2:737294/100‑1 (MQ=255) aCCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATAACCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCACGCGCTTCTGGGt > 2:1132573/1‑100 (MQ=255) tGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAGCGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTAcc < 2:1396343/100‑1 (MQ=255) ggCGCGTAAGCTGGATATGAACAGCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCg > 2:2316869/1‑100 (MQ=255) gcgTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTAt > 1:64006/1‑100 (MQ=255) cgTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATg < 1:2545621/100‑1 (MQ=255) gTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGT‑AAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGac > 2:2056288/1‑100 (MQ=255) tAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGac > 2:1702489/1‑100 (MQ=255) tAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGac < 1:504571/100‑1 (MQ=255) tGGATATGAGCAGCAGGGCGCGACTTTATACCCGCGGGATTGACGGTAAAACGGATTACTACAGTGTGGCCCCGCGGTTCTGGGTACCGTATGACACgg > 2:217769/1‑99 (MQ=255) tGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGgt > 2:229377/1‑100 (MQ=255) atatGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTa > 2:79805/1‑100 (MQ=255) atGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACa < 1:2300844/100‑1 (MQ=255) aaCAGCATGGCGCGACTTTATACCCGCGTGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGc < 2:361957/99‑1 (MQ=255) aCAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGt > 2:1642411/1‑100 (MQ=255) aGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCg < 2:64006/100‑1 (MQ=255) gCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCga > 2:695940/1‑100 (MQ=255) gCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCga > 1:1202688/1‑100 (MQ=255) gCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCATTTCTGGGTACCGTATGACACGGTGTACAGCGTCga > 1:2018881/1‑100 (MQ=255) gCATGGCGCGACCTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCg > 2:946039/1‑99 (MQ=255) ggCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCGAGaaa < 2:2146249/100‑1 (MQ=255) ggCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCGAGaaa > 1:933028/1‑100 (MQ=255) gcgcGGCTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTCTCTGGGTACCGTATGAAATGGTGTACAGCGTCGAGaaaa > 2:2238472/1‑100 (MQ=255) gcgACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCGAGAAAAAt < 1:241491/100‑1 (MQ=255) gcgACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCGAGAAAAAt < 1:156852/100‑1 (MQ=255) cgACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCGAGAAAAAt > 1:107849/1‑99 (MQ=255) cTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCGAGAAAAATGAAg > 1:1074259/1‑100 (MQ=255) | GCGAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTTCTGGGTACCGTATGACACGGTGTACAGCGTCGAGAAAAATGAAG > NC_000913/1204843‑1205040 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |