Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | USA300TCH1516_ALE | 665,392 | T→G | 100% | L79* (TTA→TGA) | paiA_1 → | Spermidine/spermine N(1)‑acetyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 665,392 | 0 | T | G | 97.2% | 68.6 / ‑6.0 | 21 | L79* (TTA→TGA) | paiA_1 | Spermidine/spermine N(1)‑acetyltransferase |
Reads supporting (aligned to +/- strand): ref base T (0/0); major base G (14/6); minor base C (1/0); total (15/6) | |||||||||||
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 | |||||||||||
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. |
GACAACAGTTAATCACTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGATAATAAAAATGAAGAA > USA300TCH1516_ALE/665255‑665528 | gACAACAGTTAATCACTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATgata > 1:637616/1‑141 (MQ=255) gACAACAGTTAATCACTTGTTAGAAGCTTAATATAATGAACATCATTTATAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAGTgata > 1:157202/1‑141 (MQ=255) acaacaGTTCATCACTTGTTAGAAGCTTCATCTCTGGTGCAGCTGGTCAAGAAAAGACTTCAAGGACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGatat > 1:596676/1‑141 (MQ=255) caacaGTTAATCACTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATAtt > 1:502704/1‑141 (MQ=255) aCTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCAtt < 1:282485/141‑1 (MQ=255) aCTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCAtt < 1:770591/141‑1 (MQ=255) cTTGTTAGAAGCTTCGTATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTAAttt > 1:603413/1‑140 (MQ=255) ttAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGtt > 1:177992/1‑141 (MQ=255) ttAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACGTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGtt > 1:133129/1‑141 (MQ=255) aGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAAcc < 1:366057/141‑1 (MQ=255) tcatTTAAAGAAAAGAATTCAAGAACAATTATTCTTAGGCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAgg > 1:455146/1‑141 (MQ=255) gaaaagaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTAc < 1:807400/141‑1 (MQ=255) aaGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCCTTGTTaa > 1:403543/1‑141 (MQ=255) aaTTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTACCTTTATTTACGGTGACGAAACTTTCTTATCAGCTAATTATGTTAATCCAGAATCGCAACATACAGTTTATGGTACAGCATTGTTAAATGAAg > 1:801615/1‑141 (MQ=255) attagtCGTAGTCGGTGAAGAAGTTAATGACATCGTTGGCTTTGCTAACTTTATTTACCGTGAAGAATGAGATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCAGTGTTAAATGAAgg < 1:480291/141‑1 (MQ=255) gTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGt > 1:424115/1‑141 (MQ=255) aagGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATa > 1:145291/1‑141 (MQ=255) gACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGaa > 1:716548/1‑141 (MQ=255) tatttaCGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAGTTAACGCTGACAAGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGataa < 1:354994/141‑1 (MQ=255) tttaCGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGataata > 1:7040/1‑141 (MQ=255) aCGGCGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCCCCACACACCGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGATAATaaaa > 1:555803/1‑141 (MQ=255) aagaaTGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGGTAATAAAAATgaag > 1:314954/1‑141 (MQ=255) gaaTGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGATAATAAAAATgaagaa > 1:721107/1‑141 (MQ=255) | GACAACAGTTAATCACTTGTTAGAAGCTTCATATAATGAACATCATTTAAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGATAAATTTGAAGGTGTTTACTTAGAAGTAGATAATAAAAATGAAGAA > USA300TCH1516_ALE/665255‑665528 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |