| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913_3_bme_pgi | 1,529,203 | C→A | pseudogene (1302/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,203 | 0 | C | A | 92.3% | 19.1 / ‑4.3 | 13 | pseudogene (1302/2037 nt) | rhsE | pseudogene, Rhs family |
| Reads supporting (aligned to +/- strand): ref base C (0/1); new base A (8/4); total (8/5) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.85e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01 | |||||||||||
GAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGG > NC_000913_3_bme_pgi/1529153‑1529348 | gAATACACACCGGCGCGAAACGATCTACTTTATCACTGCGAACCCTGCGGCGTCCTGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCt < 1:716237/149‑1 (MQ=11) ccGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGcc > 1:38508/1‑149 (MQ=11) ctCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCa > 1:286750/2‑149 (MQ=11) tttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGAGGCGAGGGATGTATGATGAATGGGGCAACCAGATGAAAAAGAAGAAGCGGCCTCAACTGCACAAGCCTTAACGTCAGCCACGgcgg > 1:335409/1‑147 (MQ=2) aCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGc > 1:98310/1‑123 (MQ=255) aCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGc < 2:98310/123‑1 (MQ=255) gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCAtgatga < 1:498066/149‑1 (MQ=255) accaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTGAATGAGGAGAACCCGCATCACCTGCCCCCGCCGTACCGCCTGCCAGGGCAGCAGCATGATgag < 1:572774/149‑1 (MQ=14) caccGGGGACTGCCGCTGTCGCTGATCAGCGAAGACGGCAATACGGCGTGGCGCGCGGAGTATGATGAATGGGGCAACTAGCTTAATGAGGAGAACTCGCATCACCTGCACCAGCCGTACCGTTAGCACAGGCAGCAGCATGATgagga > 1:675467/1‑149 (MQ=11) caccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGAGTGGCGCGGGGAGTATGATACATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTATGCCAGGGCACCAGCATGATgagga > 1:661771/1‑149 (MQ=14) accGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATgaggag > 2:679372/1‑149 (MQ=255) gggACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgg > 1:262037/1‑149 (MQ=255) gggACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGCGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgg < 1:202178/149‑1 (MQ=255) | GAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGG > NC_000913_3_bme_pgi/1529153‑1529348 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |