| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 2,994,415 | G→A | pseudogene (34/60 nt) | pbl → | pseudogene, peptidoglycan‑binding enzyme family |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,994,415 | 0 | G | A | 98.5% | 202.8 / ‑6.6 | 65 | pseudogene (34/60 nt) | pbl | pseudogene, peptidoglycan‑binding enzyme family |
| Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (38/26); minor base C (0/1); total (38/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.15e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGA > NC_000913/2994320‑2994550 | gggACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCt > 2:257211/1‑110 (MQ=255)gggACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCt < 1:257211/110‑1 (MQ=255) ttATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct < 1:587887/130‑1 (MQ=255) ttATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct > 2:587887/1‑130 (MQ=255) ttATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTtat > 1:147069/1‑139 (MQ=255) tGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctc > 2:1972/1‑123 (MQ=255) tGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctc < 1:1972/123‑1 (MQ=255) aaaaGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGtttt > 1:574385/1‑139 (MQ=255) gaCAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACATCCATACGTAAAGCTGTTTTGATCaa < 2:750918/139‑1 (MQ=255) aGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTc > 1:110251/1‑139 (MQ=255) atatGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATa < 2:574385/139‑1 (MQ=255) cAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACtgtg < 2:110251/139‑1 (MQ=255) tatatatTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACAt > 1:33344/1‑139 (MQ=255) tGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAAttt > 1:678818/1‑139 (MQ=255) tatataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACTACTGTGAACATCTAATttct > 2:218839/1‑139 (MQ=255) tatataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATAACACTGTGAACATCTAATttct > 1:25252/1‑139 (MQ=255) tatataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTATATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATttct > 2:48942/1‑139 (MQ=255) atataCCAGGCTTAATGAACTAAACACTCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct > 1:59675/1‑62 (MQ=255) atataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct < 2:59675/62‑1 (MQ=255) atataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTtata < 1:379062/72‑1 (MQ=255) atataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTtata > 2:379062/1‑72 (MQ=255) atataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaaa > 2:183369/1‑81 (MQ=255) atataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaaa < 1:183369/81‑1 (MQ=255) tataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct < 1:97167/61‑1 (MQ=255) tataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct > 2:97167/1‑61 (MQ=255) ataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCa < 2:128072/139‑1 (MQ=255) ataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCa > 2:245738/1‑139 (MQ=255) ccAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGAt < 1:832124/101‑1 (MQ=255) ccAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGAt > 2:832124/1‑101 (MQ=255) cAGGCTTAATGAACTACACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct < 1:612639/56‑1 (MQ=255) cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct > 2:84053/1‑56 (MQ=255) cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct < 1:771649/56‑1 (MQ=255) cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct < 1:84053/56‑1 (MQ=255) cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct > 2:771649/1‑56 (MQ=255) cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct > 2:612639/1‑56 (MQ=255) cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCtttttt < 1:12314/61‑1 (MQ=255) cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCtttttt > 2:12314/1‑61 (MQ=255) cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGat > 1:47065/1‑139 (MQ=255) aGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATc < 2:198623/100‑1 (MQ=255) aGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATc > 1:198623/1‑100 (MQ=255) aGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCaa < 2:252700/102‑1 (MQ=255) aGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCaa > 1:252700/1‑102 (MQ=255) gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct < 2:531648/53‑1 (MQ=255) gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct > 1:23829/1‑53 (MQ=255) gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct > 2:408136/1‑53 (MQ=255) gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct < 1:144876/53‑1 (MQ=255) gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct > 1:531648/1‑53 (MQ=255) gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct > 1:109734/1‑53 (MQ=255) gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct < 2:23829/53‑1 (MQ=255) gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct < 2:109734/53‑1 (MQ=255) gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct > 2:144876/1‑53 (MQ=255) gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct < 1:408136/53‑1 (MQ=255) gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACt < 2:175496/102‑1 (MQ=255) gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACt > 1:175496/1‑102 (MQ=255) gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt > 1:491022/1‑139 (MQ=255) gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt > 2:665205/1‑139 (MQ=255) gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATAAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACt > 1:263550/1‑102 (MQ=255) gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATAAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACt < 2:263550/102‑1 (MQ=255) aaTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctc > 1:41957/1‑48 (MQ=255) aaTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctc < 2:41957/48‑1 (MQ=255) aaTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCGGTTTTGATCAACTCGATAGCACTGTGAACATCTAATTTCTTCATGATATTGAGt > 1:506847/1‑139 (MQ=255) aTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTc > 2:505247/1‑139 (MQ=255) tGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCt > 1:88521/1‑139 (MQ=255) aCTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATg < 1:505247/139‑1 (MQ=255) cTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAGCTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGa > 1:343280/1‑139 (MQ=255) | GGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGA > NC_000913/2994320‑2994550 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |