| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | REL606 | 666,130:1 | +T | 100% | pseudogene (358/387 nt) | gltL ← | ATP‑binding protein of glutamate/aspartate transport system; b0652_1 |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 666,130 | 1 | . | T | 96.7% | 97.0 / ‑2.9 | 30 | pseudogene (358/387 nt) | gltL | ATP‑binding protein of glutamate/aspartate transport system; b0652_1 |
| Reads supporting (aligned to +/- strand): ref base . (1/0); new base T (14/15); total (15/15) | |||||||||||
| 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.63e-01 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
| Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. | |||||||||||
CTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA‑T‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACAGCTCGAAATGCTGGAATACCATCCCGACGCGGGAACGCAGCTTTGCC > REL606/665998‑666272 | cTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCa > 2:182506/1‑150 (MQ=255) gAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCgg > 2:446563/1‑150 (MQ=255) atatTGAATACTACAACAGCAGAAGAATTAGCCTGAATTTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA‑TGTTCAGGGCTTTTTCACGAGCCGGCGCt > 2:408989/1‑150 (MQ=255) atatTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCt < 1:181981/150‑1 (MQ=255) atTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCttt > 2:140619/1‑150 (MQ=255) ctacAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTACCGCGtt > 1:275438/1‑150 (MQ=255) cagcagAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCa < 2:638079/150‑1 (MQ=255) gaaTTAGCCTGGAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCAcc < 1:104836/150‑1 (MQ=255) gaaTTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCAcc < 2:8900/150‑1 (MQ=255) gaaTTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCAcc > 2:46963/1‑150 (MQ=255) gaaTTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCAcc > 1:8871/1‑150 (MQ=255) aaaTTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGt < 1:382018/150‑1 (MQ=255) aTTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCa > 2:550319/1‑150 (MQ=255) ggTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTc > 2:123181/1‑150 (MQ=255) gTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCg > 2:557372/1‑150 (MQ=255) ccAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGa > 1:611250/1‑150 (MQ=255) gAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATg < 1:140217/150‑1 (MQ=255) tCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGa < 1:255955/150‑1 (MQ=255) ggAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAAc < 2:302056/150‑1 (MQ=255) gACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACAGCTCGa > 2:461993/1‑150 (MQ=255) aCCTATATGCCTCGTGTTTAACTGTCCAACTCTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCTACAGATGGGGGAACAGCTCGaa < 2:431206/150‑1 (MQ=38) tataTGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACAGCTCGAAATg > 1:620840/1‑150 (MQ=255) cttgtTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACATCTCGAAATGCTGGAATAc < 1:555763/148‑1 (MQ=255) cGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACAGCTCGAAATGCTGGAATAc < 2:622660/150‑1 (MQ=255) gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACAGCTCGAAATGCTGGAATAcc > 1:355695/1‑150 (MQ=255) gtTTAACTGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACAGCTCGAAATGCTGGAATACCAt < 2:637663/150‑1 (MQ=255) tGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGTTTTTTCACTAGCCGGGGCTTTATCGCGTTTAAGCACTTTAGCCTGCTCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACAGCTCGAAATGCTGGAATACCATACCGAcg < 1:102739/150‑1 (MQ=255) tGTCCAACTTTTTGGGGTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACAGCTCGAAATGCTGGAATACCATCCCGAcg < 2:613042/150‑1 (MQ=255) ggTCAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACAGCTCGAAATGCTGGAATACCATCCCGACGCGGGAACGCAGCttt > 2:65896/1‑150 (MQ=255) cAGTACATT‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACAGCTCGAAATGCTGGAATACCATCCCGACGCGGGAACGCAGCTTTGcc < 1:122822/150‑1 (MQ=255) | CTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA‑T‑TTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACAGCTCGAAATGCTGGAATACCATCCCGACGCGGGAACGCAGCTTTGCC > REL606/665998‑666272 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |