| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913_3_bme_pgi | 1,529,218 | C→T | pseudogene (1317/2037 nt) | rhsE → | pseudogene, Rhs family |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913_3_bme_pgi | 1,529,218 | 0 | C | T | 96.4% | 64.5 / ‑5.0 | 28 | pseudogene (1317/2037 nt) | rhsE | pseudogene, Rhs family |
| Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (13/14); minor base G (1/0); total (14/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01 | |||||||||||
ACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACC > NC_000913_3_bme_pgi/1529159‑1529363 | acacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCa > 1:140636/1‑149 (MQ=11) gAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGAATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCgc > 1:317245/1‑148 (MQ=11) ctCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCa > 1:357322/2‑149 (MQ=255) tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTGTGA‑‑TGAAGGGGGCGATGAGCTTAATGAGGAGCAGCAGCATCACCTGCACCAGCCGTACCGTCTGCCAg > 1:488900/1‑149 (MQ=11) tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTCCCGccagccag > 1:59143/1‑149 (MQ=11) atctttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAggg < 2:451167/146‑1 (MQ=255) ctttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgca > 1:128245/2‑149 (MQ=255) tttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcag < 2:604746/149‑1 (MQ=255) ttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcagc < 1:233610/149‑1 (MQ=255) ttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcagc < 1:417731/149‑1 (MQ=255) ttAGCACTGCGACCACCGTGGACGGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcagc < 1:618301/149‑1 (MQ=14) tatcACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcagca < 2:557120/146‑1 (MQ=255) tcACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCACTACGGCGCGGCGCGGGGAGTATGA‑‑TGAATGGGGCAAAAAGCTTAACGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCACGGgcagcggc > 1:351501/2‑149 (MQ=14) tcACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcagcagc > 1:75864/2‑149 (MQ=255) aCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGAATGA‑‑TGAATGGGGCAACCAGCTAAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGGCGGCACGGgcagcagcag > 2:748752/1‑148 (MQ=14) gCGACCACCGGGGCCTGTCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCatgat < 1:749588/149‑1 (MQ=255) gCGACCACCGGGGCCTGCGGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCatgat < 2:617601/149‑1 (MQ=255) gCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCatgat < 1:5597/149‑1 (MQ=255) gCGACCACCGGGGACTGCCGCCGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCatgat < 1:76945/149‑1 (MQ=255) gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACAGCAATCCGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCCACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCATCTCCCAGGGCAGAAGCgtgagaa > 2:125413/1‑146 (MQ=14) accaccGGGGACTGCCGCGGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATgag < 1:287953/149‑1 (MQ=255) caccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTACTGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATgagga > 1:683402/1‑149 (MQ=255) gacggCCGCTGGAGCGTACCAGCGAGTACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAAGGAGGAGAACCCGGATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgggg < 2:344737/145‑1 (MQ=17) tGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATAATCTCAATCGGG‑‑ACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTg < 1:371582/149‑1 (MQ=17) cGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTAc > 2:58661/1‑149 (MQ=255) cGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTAc > 1:252595/1‑149 (MQ=255) gtcttcgcTGATAAGCGCAGAAGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGAGCAACCAGATTAATGAGGAGAACCCGCATAACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGCGGAGTCGGGCCTGTACt > 2:306967/6‑149 (MQ=17) gcTTATCAGCGATGAGGGAAATGCGGGGGTGCCCCGGGAGTAGGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAAcc < 2:316908/149‑1 (MQ=17) | ACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGA‑‑TGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACC > NC_000913_3_bme_pgi/1529159‑1529363 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |