| 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 | 95.0% | 41.0 / NA | 20 | pseudogene (1317/2037 nt) | rhsE | pseudogene, Rhs family |
| Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (9/10); total (9/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACC > NC_000913_3_bme_pgi/1529151‑1529363 | tggAATACACACCGGCTCGAAAAGTTTATCTTTCTCACTGGGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGCGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAc < 2:400801/148‑1 (MQ=11) cacCGGCGCCAAAAGCTCATCTTTATCACTGCGACCACCGGGGACGGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAg < 1:329192/149‑1 (MQ=11) ggCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGt > 1:608419/1‑149 (MQ=11) gAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGACGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGTACCCGCATCACCTGCACCAGCCGTACCGt > 1:272033/1‑149 (MQ=11) gAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGt < 2:20509/149‑1 (MQ=14) gctCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCTAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGcc > 1:557350/3‑149 (MQ=11) tatcACTGCGACCACCGGGGACGGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcagca < 2:435956/146‑1 (MQ=255) atcACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGGGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcagcag > 1:349088/3‑149 (MQ=255) tcACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcagcagc > 1:681081/2‑149 (MQ=255) aCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCat > 1:652779/1‑149 (MQ=255) cTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCatg > 1:435451/1‑149 (MQ=255) gCGACCACCGGGGACTGCCGCTGGCTCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCatgat < 1:112424/149‑1 (MQ=255) cGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCatgatg < 1:576209/149‑1 (MQ=255) gACCCCCGGGGACGGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCAtgatga < 1:429357/149‑1 (MQ=255) gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCAtgatga > 1:666448/1‑149 (MQ=255) gagcGGGGACTGCCGCTGGGGGTTATCAGGGAAGACGGCAATATGGCGGGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATgagga < 2:217098/146‑1 (MQ=14) accGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGGACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGCCAGCCGCATGATgaggag > 2:476615/1‑149 (MQ=14) ggacTGCCGGTGCCGCTTATCAGCGAAGACGGCCATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCggg < 1:537873/146‑1 (MQ=14) gcTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAAcc < 1:773381/149‑1 (MQ=255) gcTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAAcc < 1:644571/149‑1 (MQ=255) | CGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACC > NC_000913_3_bme_pgi/1529151‑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: ‑ |