| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | USA300TCH1516_ALE | 644,035 | C→T | 23.6% | V207I (GTA→ATA) | USA300HOU_RS03070 ← | Putative heme‑dependent peroxidase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | USA300TCH1516_ALE | 644,035 | 0 | C | T | 23.6% | 24.0 / 4.7 | 17 | V207I (GTA→ATA) | USA300HOU_RS03070 | Putative heme‑dependent peroxidase |
| Reads supporting (aligned to +/- strand): ref base C (4/9); new base T (2/2); total (6/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.84e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01 | |||||||||||
AGAAATCGCAAAGAATTGATCGAATTCGTTTGTGTTAATAATATGTCCTACAAAGAAACTACCGAATTCACCGTATCGTGCTGTTGTTTCATCAAAGCGCATTTCGTATACAATTTTTTTGAATTGTAATACGTCATCTGAGAACAATGTTACGCCCCATTCGAAATCATCAAACCCTACAGAACCAGTAATAAATTGTTTGATTTTGCCAGCATATTTTCTACCAATCATACCATGGTCATACATTAATTTTTGGCGT > USA300TCH1516_ALE/643904‑644162 | aGAAATCGCAAAGAATTGATCGAATTCGTTTGTGTTAATAATATGTCCTACAAAGAAACTACCGAATTCACCGTATCGTGCTGTTGTTTCATCAAAGCGCATTTCGTATACAATTTTTTTGCATTGTAATACGTCATCTga > 2:90125/1‑141 (MQ=255) aGAATTGATCGAATTCGTTTGTGTTAATAATATGTCCTACAAAGAAACTACCGAATTCACCGTATCGTGCTGTTGTTTCATCAAAGCGCATTTCGTATACAATTTTTTTGAATTGTAATACGTCATCTGAGAACAATGttc > 2:88039/1‑140 (MQ=255) ttGTGTTAATAATATGTCCTACAAAGAAACTACCGAATTCACCGTATCGTGCTGTTGTTTCATCAAAGCGCATTTCGTATACAATTTTTTTGAATTGTAATACGTCATCTGAGAACAATGTTACGCCCCATTCGAAatcat < 1:108150/141‑1 (MQ=255) taataaTATGTCCTACAAAGAAACTACCGAATTCACCGTATCGTGCTGTTGTTTCATCAAAGCGCATTTCGTATACAATTTTTTTGAATTGTAATACGTCATCTGAGAACAATGTTACGCCCCATTCGAAATCATCAAAcc < 1:215992/141‑1 (MQ=255) taataTGTCCTACAAAGAAACTACCGAATTCACCGTATCGTGCTGTTGTTTCATCAAAGCGCATTTCGTATACAATTTTTTTGAATTGTAATACGTCATCTGAGAACAATGTTACGCCCCATTCGAAATCATCAAACCCTa < 2:210090/141‑1 (MQ=255) cTACAAAGAAACTACCGAATTCACCGTATCGTGCTGTTGTTTCATCAAAGCGCATTTCGTATACAATTTTTTTGAATTGTAATACGTCATCTGAGAACAATGTTACGCCCCATTCGAAATCATCAAACCCTACAGAACCAg > 1:184187/1‑141 (MQ=255) aaGAAACTACCGAATTCACCGTATCGTGCTGTTGTTTCATCAAAGCGCATTTCGTATACAATTTTTTTGAATTGTAATATGTCATCTGAGAACAATGTTACGCCCCATTCGAAATCATCAAACCCTACAGAACCAGtaata > 1:155236/1‑141 (MQ=255) aaGAAACTACCGAATTCACCGTATCGTGCTGTTGTTTCATCAAAGCGCATTTCGTATACAATTTTTTTGAATTGTAATATGTCATCTGAGAACAATGTTACGCCCCATTCGAAATCATCAAACCCTACAGAACCAGtaata > 1:155235/1‑141 (MQ=255) cGTATCGTGCTGTTGTTTCATCAAAGCGCATTTCGTATACAATTTTTTTGAATTGTAATACGTCATCTGAGAACAATGTTACGCCCCATTCGAAATCATCAAACCCTACAGAACCAGTAATAAATTGTTTGATTTTGCCAg < 1:198244/141‑1 (MQ=255) cGTATCGTGCTGTTGTTTCATCAAAGCGCATTTCGTATACAATTTTTTTGAATTGTAATACGTCATCTGAGAACAATGTTACGCCCCATTCGAAATCATCAAACCCTACAGAACCAGTAATAAATTGTTTGATTTTGCCAg < 2:165831/141‑1 (MQ=255) tagcGTGCTGTTGTTTCATCAAAGCGCATTTCGTATACAATTTTTTTGAATTGTAATACGTCATCTGAGAACAATGTTACGCCCCATTCGAAATCATCAAACCCTACCGAACCAGTAATAAATTGTTTGATTTTGCCAGCa < 1:187981/138‑1 (MQ=255) atcaAAGCGCATTTCGTATACAATTTTTTTGAATTGTAATACGTCATCTGAGAACAATGTTACGCCCCATTCGAAATCATCAAACCCTACAGAACCAGTAATAAATTGTTTGATTTTGCCAGCATATTTTCTACCAATCAt < 1:202536/141‑1 (MQ=255) aaaGCGCATTTCGTATACAATTTTTTTGAATTGTAATACGTCATCTGAGAACAATGTTACGCCCCATTCGAAATCATCAAACCCTACAGAACCAGTAATAAATTGTTTGATTTTGCCAGCATATTTTCTACCAATCATAcc < 1:90125/141‑1 (MQ=255) ttttttGAATTGTAATATGTCATCTGAGAACAATGTTACGCCCCATTCGAAATCATCAAACCCTACAGAACCAGTAATAAATTGTTTGATTTTGCCAGCATATTTTCTACCAATCATACCATGGTCATACATTAATTTTTg < 2:155235/141‑1 (MQ=255) ttttttGAATTGTAATATGTCATCTGAGAACAATGTTACGCCCCATTCGAAATCATCAAACCCTACAGAACCAGTAATAAATTGTTTGATTTTGCCAGCATATTTTCTACCAATCATACCATGGTCATACATTAATTTTTg < 2:155236/141‑1 (MQ=255) ttGAATTGTAATACGTCATCTGAGGACAATGTTACTCCCCATTCGAAATCATCAAACCCTACAGAACCAGTAATATATTGTTTGATTTTGCCAGCATATTTTCTACCAATCATACCATGGTCATACATTAATTTTTGGCGt < 2:160836/141‑1 (MQ=255) ttGAATTGTAATACGTCATCTGAGAACAATGTTACGCCCCATTCGAAATCATCAAACCCTACAGAACCAGTAATAAATTGTTTGATTTTGCCAGCATATTTTCTACCAATCATACCATGGTCATACATTAATTTTTGGCGt > 1:145488/1‑141 (MQ=255) | AGAAATCGCAAAGAATTGATCGAATTCGTTTGTGTTAATAATATGTCCTACAAAGAAACTACCGAATTCACCGTATCGTGCTGTTGTTTCATCAAAGCGCATTTCGTATACAATTTTTTTGAATTGTAATACGTCATCTGAGAACAATGTTACGCCCCATTCGAAATCATCAAACCCTACAGAACCAGTAATAAATTGTTTGATTTTGCCAGCATATTTTCTACCAATCATACCATGGTCATACATTAATTTTTGGCGT > USA300TCH1516_ALE/643904‑644162 |
| 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: ‑ |