| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | USA300TCH1516_ALE | 673,010 | T→C | 12.1% | G78G (GGT→GGC) | btuF → | Vitamin B12‑binding protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | USA300TCH1516_ALE | 673,010 | 0 | T | C | 12.1% | 79.0 / 2.9 | 33 | G78G (GGT→GGC) | btuF | Vitamin B12‑binding protein |
| Reads supporting (aligned to +/- strand): ref base T (16/13); new base C (2/2); total (18/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 = 9.28e-01 | |||||||||||
GGCTCGTCTAAATCGCCGTACCATAGAATTGTTTCGTTAATGCCTAGTAATACTGAAATTTTATATGAATTAGGATTAGGTAAATACATAGT‑TGGTGTTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATGAGTCGCAAAAGGCAACTGCTAATAAAGTATTGTCATCATTAGAGAAACAAGGCATCAAAGT > USA300TCH1516_ALE/672873‑673147 | ggCTCGTCTAAATCGCCGTACCATAGAATTGTTTCGTTAATGCCTAGTAATACTGAAATTTTATATGAATTAGGATTAGGTAAATACATAGT‑TGGTGTTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGTAAg < 1:2842/141‑1 (MQ=255) cGTCTAAATCGCCGTACCATAGAATTGTTTCGTTAATGCCAAGTAATACTGAAATTTTATATGAATTAGGATTAGGTAAATACATAGT‑CGGTGTTTCAACGGTTGATGATGATCCAAACGATGTGAAAAAGGGTAAGAAAc > 2:248324/1‑141 (MQ=255) tCTAAATAGCCAAACAATAAAATTGTTTCGTTAATGTCTAGTAATACTGAAATTTTATATGAATTAGGATTAGGTAAATACACAGT‑TGGTGTTTCAACGGTTGATGATTATCAAAAAAAAAAAAAAAAGGGTAAGAAACaa < 2:118012/141‑1 (MQ=255) cGCCGTACCATAGAATTGTTTCGTTAATGCCTAGTAATACTGAAATTTTATATGAATTAGGATTAGGTAAATACATAGT‑TGGTGTTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATg < 2:238254/141‑1 (MQ=255) cGCCGGACCATAGAATTGTTTAGTAAATGGCTATTAATACTCAAATTTTAGATGAATTAGGATTAGGTAAATACATAGT‑TGGTGTTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATg < 1:150083/141‑1 (MQ=255) actaaaGAATTGTTTCGTTAATGCCTAGTAATACTGAAATTTTATATGAATTAGGATTAGGTAAATACATAGT‑TGGTGTTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGCAAGAAACAATTTGATGCTTTGa > 2:102642/6‑141 (MQ=255) aCCATAGAATTGTTTCGTTAATGCCTAGTAATACTGAAATTTTATATGAATTAGGATTAGGTAAATACATAGT‑TGGTGTTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGCAAGAAACAATTTGATGCTTTGa > 2:225439/1‑141 (MQ=255) aCCATAGAATTGTTTCGTTAATGCCTAGTAATACTGAAATTTTATATGAATTAGGATTAGGTAAATACATAGT‑TGGTGTTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGAGTAAGAAAGAATTTGATGCTGTGa > 1:221740/1‑141 (MQ=255) gTAATACTGAAATTTTATATGAATTAGGATTAGGTAATTACATTGT‑TGGTGTTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGAAACTTTTCAAgg < 1:278826/141‑1 (MQ=255) gTAATACTGAAATTTTATATGAATTAGGATTAGGTAAATACATAGT‑TGGTGTTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAgg < 1:254923/141‑1 (MQ=255) aaTATTGAAATTTTATATGAATTAGGATTTGGTAAATACAGTGG‑TGGTGTTTGAACGGTTGATGATTATTCAAAAAAAGAGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCa < 2:203063/141‑1 (MQ=255) tataTGAATTAGGATTAGGTAAATACATAGT‑TGGTGTTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTaa < 1:44498/141‑1 (MQ=255) atatGAATTAGGATTAGGTAAATACATAGT‑TGGTGTTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCACAGCCAGATCTAAt > 1:10489/1‑141 (MQ=255) atatGAATTAGGATTAGGTAAATACATAGT‑TGGTGTTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAAt > 1:203849/1‑141 (MQ=255) gcttaggTAAATACATAGT‑TGGTGTTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATg > 2:46835/3‑141 (MQ=255) attaggTAAA‑ACCTATTCTGGTGTTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATGa > 2:183975/1‑141 (MQ=255) aaaTACATAGT‑TGGTGTTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTCAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATGAGTCGCaa > 1:388316/1‑141 (MQ=255) aaaTACATAGT‑TGGTGTTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATGAGTAGCaa > 2:137381/1‑141 (MQ=255) aTAGT‑TGGTGTTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATGAGTCGCAAAAGGCa > 1:280995/1‑141 (MQ=255) gtgtTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATGAGTCGCAAAAGGCAACTGCta > 2:222481/1‑141 (MQ=255) gtgtTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATGAGTCGCAAAAGGCAACTGCta < 1:201292/141‑1 (MQ=255) tgtTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGCAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATGAGTCGCAAAAGGCAACTGCtaa < 1:102642/141‑1 (MQ=255) tgtTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGCAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATGAGTCGCAAAAGGCAACTGCtaa < 1:225439/141‑1 (MQ=255) ttCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATGAGTCGCAAAAGGCAACTGCtaataa > 2:17548/1‑141 (MQ=255) cAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATGAGTCGCAAAAGGCAACTGCTAATAAAg > 1:100285/1‑141 (MQ=255) aCGGTTGATGATTATCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATGAGTCGCAAAAGGCAACTGCTAATAAAGTa > 2:270481/1‑141 (MQ=255) gatgatTATCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATGAGTCGCAAAAGGCAACTGCTAATAAAGTATTGtca < 2:203849/141‑1 (MQ=255) atgatTAGCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATGAGTCGCAAAAGGCAACTGCTAATAAAGTATTGtcat < 2:315267/141‑1 (MQ=255) aTCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATGAGTCGCAAAAGGCAACTGCTAATAAAGTATTGTCATCATTag > 1:18057/1‑141 (MQ=255) aTCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATGAGTCGCAAAAGGCAACTGCTAATAAAGTATTGTCATCATTaa > 2:99584/1‑140 (MQ=255) aTGTGAAAAAGGGTAAGAAAAAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATGAGTCGCAAAAGGCAACTGCTAATAAAGTATTGTCATCATTAGACAAACAAg < 1:248324/141‑1 (MQ=255) aGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATGAGTCGCAAAAGGCAACTGCTAATAAAGTATTGTCATCATGAGAGAAACAAGGCATCAAag > 1:367248/1‑141 (MQ=255) gggTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATGAGTCGCAAAAGGCAACTGCTAATAAAGTATTGTCATCATTAGAGAAACAAGGCATCAAagt < 1:222481/141‑1 (MQ=255) | GGCTCGTCTAAATCGCCGTACCATAGAATTGTTTCGTTAATGCCTAGTAATACTGAAATTTTATATGAATTAGGATTAGGTAAATACATAGT‑TGGTGTTTCAACGGTTGATGATTATCCAAAAGATGTGAAAAAGGGTAAGAAACAATTTGATGCTTTGAATCTAAATAAAGAGGAACTTTTAAAGGCAAAGCCAGATCTAATTCTTGCGCATGAGTCGCAAAAGGCAACTGCTAATAAAGTATTGTCATCATTAGAGAAACAAGGCATCAAAGT > USA300TCH1516_ALE/672873‑673147 |
| 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: ‑ |