Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,489,365 | Δ1 bp | coding (689/831 nt) | glpG ← | predicted intramembrane serine protease |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,489,364 | 0 | A | . | 96.9% | 146.7 / ‑1.1 | 32 | coding (690/831 nt) | glpG | predicted intramembrane serine protease |
Reads supporting (aligned to +/- strand): ref base A (1/0); new base . (21/10); total (22/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01 |
ATGTGTGCTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATAATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCA > REL606/3489290‑3489434 | aTGTGTGCTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATat < 1:1471347/76‑2 (MQ=255) gtgtgCTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATatta < 2:1644004/76‑4 (MQ=255) gtgtgCTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATatta < 1:1453540/76‑4 (MQ=255) cTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGAT‑ATTAACCCa > 2:739813/1‑76 (MQ=255) cTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGAT‑ATTAACCCa < 2:274623/76‑1 (MQ=255) tCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGAT‑ATTAACCCAc > 1:373500/1‑76 (MQ=255) gTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGAT‑ATTAACCCACGtt > 2:1417735/1‑76 (MQ=255) tCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGAT‑ATTAACCCACGTTGc < 1:1067327/76‑1 (MQ=255) ccATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGAT‑ATTAACCCACGTTGCAgg < 1:508127/76‑1 (MQ=255) cATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGAT‑ATTAACCCACGTTGCAGGt > 2:779583/1‑76 (MQ=255) tCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGAT‑ATTAACCCACGTTGCAGGTaa > 1:1303523/1‑76 (MQ=255) cGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGAT‑ATTAACCCACGTTGCAGGTaaa > 2:1026013/1‑76 (MQ=255) gACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGAT‑ATTAACCCACGTTGCAGGTAAAt < 1:1125554/76‑1 (MQ=255) aTCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGAT‑ATTAACCCACGTTGCAGGTAAATGcc > 1:666837/1‑76 (MQ=255) aTCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGAT‑ATTAACCCACGTTGCAGGTAAATGcc > 1:77047/1‑76 (MQ=255) ccAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGAT‑ATTAACCCACGTTGCAGGTAAATGCCACt > 1:1271642/1‑76 (MQ=255) aaacaaaTCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGAT‑ATTAACCCACGTTGCAGGTAAATGCCACttt < 2:1367418/76‑1 (MQ=255) cATCCGGCGACAATCCAGATCAGCGCAAAGAT‑ATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTc > 1:1435149/1‑76 (MQ=255) aTCCGGCGACAATCCAGATCAGCGCAAAGAT‑ATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCg > 2:1210834/1‑76 (MQ=255) tCCGGCGACAATCCAGATCAGCGCAAAGAT‑ATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGc < 1:194319/76‑1 (MQ=255) ccGGCGACAATCCAGATCAGCGCAAAGATAATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGc > 1:1140220/1‑76 (MQ=255) cGACAATCCAGATCAGCGCAAAGAT‑ATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGt < 2:332877/76‑1 (MQ=255) aCAATCCAGATCAGCGCAAAGAT‑ATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAg < 1:1043665/76‑1 (MQ=255) aTCCAGATCAGCGCAAAGAT‑ATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCa > 2:796936/1‑76 (MQ=255) cAGATCAGCGCAAAGAT‑ATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGAc > 1:1048305/1‑76 (MQ=255) aGATCAGCGCAAAGAT‑ATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACg > 1:815173/1‑76 (MQ=255) cAGCGCAAAGAT‑ATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGc > 1:1089649/1‑76 (MQ=255) gcgcAAAGAT‑ATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGccc < 1:1140167/76‑1 (MQ=255) gcgcAAAGAT‑ATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGccc > 1:1028550/1‑76 (MQ=255) gcAAAGAT‑ATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTCGCCCac > 2:71634/1‑75 (MQ=255) gcAAAGAT‑ATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCAt > 2:116515/1‑76 (MQ=255) gcAAAGAT‑ATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCAt > 2:134500/1‑76 (MQ=255) cAAAGAT‑ATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATc < 1:296749/76‑1 (MQ=255) aaaGAT‑ATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCa > 1:1656595/1‑76 (MQ=255) aaaGAT‑ATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCa > 2:781445/1‑76 (MQ=255) | ATGTGTGCTCCGTTCGCCATCGACATCCCAAACAAATCAAACCATCCGGCGACAATCCAGATCAGCGCAAAGATAATTAACCCACGTTGCAGGTAAATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCA > REL606/3489290‑3489434 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |