| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | USA300TCH1516_ALE | 108,148 | A→G | 11.8% | K1037R (AAA→AGA) | recD → | RecBCD enzyme subunit RecD |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | USA300TCH1516_ALE | 108,148 | 0 | A | G | 11.8% | 90.7 / 5.0 | 34 | K1037R (AAA→AGA) | recD | RecBCD enzyme subunit RecD |
| Reads supporting (aligned to +/- strand): ref base A (9/21); new base G (2/2); total (11/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.80e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01 | |||||||||||
TATTTTGTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACTGGATATTAAATGCAATAAGTGG > USA300TCH1516_ALE/108009‑108288 | tATTTTGTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGaaa < 2:154309/141‑1 (MQ=255) ttttGTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAAcc < 1:56091/141‑1 (MQ=255) ttttGTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAAcc < 1:443731/141‑1 (MQ=255) ttttGTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTGGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATATAGATGAAAcc < 1:60695/141‑1 (MQ=255) gTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCAttt > 1:293801/1‑141 (MQ=255) gTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCAttt > 1:432536/1‑141 (MQ=255) taataCCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAGACCATTTTATGAGACGAtttt > 1:71913/1‑141 (MQ=255) taataCCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAGACCATTTTATGAGACGAtttt > 1:302309/1‑141 (MQ=255) ccccAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTaaa < 2:359129/140‑1 (MQ=255) aCCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTaaa < 2:84766/141‑1 (MQ=255) aCCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTaaa < 2:269728/141‑1 (MQ=255) ccAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTaaaga > 2:189597/1‑141 (MQ=255) aGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTaaagaaag > 1:227057/1‑141 (MQ=255) tGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAgaaagaaa < 1:143909/141‑1 (MQ=255) ggTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAAt < 1:442946/141‑1 (MQ=255) gTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATg < 1:216517/141‑1 (MQ=255) gTGCTATGAAAAGGTCATGAGAAAAACCCAAATTGTTAAAAGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAACGAATACAGCTGAAACCATCTTATGAGACGATTTTTAAAGCAAGAAATg < 2:298691/141‑1 (MQ=255) gTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATAc < 2:372951/141‑1 (MQ=255) gTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATAc < 1:460406/141‑1 (MQ=255) ccAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGtatata < 2:432536/141‑1 (MQ=255) ccAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAAAAAAAAAAAAAAAAGtatata < 2:293801/141‑1 (MQ=255) aCTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAgg > 1:221888/1‑141 (MQ=255) cTTGTTAAACGTTGCAGTGACAAGAGGAAAGAGAGAGGTGTGTGGGATGGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGa > 2:20480/1‑141 (MQ=255) ttGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAAAAAAAAAAAAAGTATATAGAGGAAGTTATACAtt < 2:364321/141‑1 (MQ=255) ttGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAAAAACAAAAAAAAAAGTATATAGAGGAAGTTATACAtt < 2:418026/141‑1 (MQ=255) tGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACAttt < 1:98540/141‑1 (MQ=255) gACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAg < 2:221888/141‑1 (MQ=255) aCAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAgg < 2:317835/141‑1 (MQ=255) aCAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAgg < 2:348154/141‑1 (MQ=255) tGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTaac < 1:298101/141‑1 (MQ=255) gACATGCAAAGAATACAGATGAGACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACt < 2:302309/141‑1 (MQ=255) gACATGCAAAGAATACAGATGAGACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAAAAAAAAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACt < 2:71913/141‑1 (MQ=255) cAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACTGGatat > 1:161637/1‑141 (MQ=255) aaGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACTGGATATTa > 2:178362/1‑141 (MQ=255) aCAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACTGGATATTAAATGCa > 1:378502/1‑141 (MQ=255) aaCCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACTGGATATTAAATGCAATAAGTgg < 1:177476/141‑1 (MQ=255) | TATTTTGTAATAGGTACTGATAATACCCAAGATGGTGCTGTGAACTGGTCATGCGAAAAACCAAACTTGTTAAACGTTGCAGTGACAAGAGCTAAGAAAGAGTTTTATGTAATTGGCGACATGCAAAGAATACAGATGAAACCATTTTATGAGACGATTTTTAAAGAAAGAAATGTAAAATAACATACAAAAAAGTATATAGAGGAAGTTATACATTTTAAAAGGAGCAAAATTGAATAATGGAGAAATTTAACAACTGGATATTAAATGCAATAAGTGG > USA300TCH1516_ALE/108009‑108288 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |