Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | USA300TCH1516_ALE | 643,642 | T→A | 21.5% | intergenic (‑290/+259) | yjjP_1 ← / ← USA300HOU_RS03070 | Inner membrane protein YjjP/Putative heme‑dependent peroxidase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 643,642 | 0 | T | A | 21.5% | 45.1 / 11.4 | 28 | intergenic (‑290/+259) | yjjP_1/USA300HOU_RS03070 | Inner membrane protein YjjP/Putative heme‑dependent peroxidase |
Reads supporting (aligned to +/- strand): ref base T (13/9); new base A (3/3); total (16/12) | |||||||||||
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.75e-01 |
TCTCTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACATATATATAAGTTTAAACATACAGATTTC > USA300TCH1516_ALE/643507‑643782 | tctcTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGt < 1:61685/141‑1 (MQ=255) gAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCGAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGt > 1:2216/1‑141 (MQ=255) ataAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTtata > 1:169905/1‑141 (MQ=255) ataAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATAAGTAAACGATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTtata > 2:213802/1‑141 (MQ=255) gAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAAATGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGa > 2:239532/1‑141 (MQ=255) aaaaTTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGATTATATGTAGTGAta > 2:236085/1‑141 (MQ=255) tGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTaa < 1:212729/141‑1 (MQ=255) aaaaaaaTTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAg < 1:67050/141‑1 (MQ=255) tttcaaatatttcaaataGGTTAAAATAATATTAGTTAACCATTACGAAGCTTGTCACCCTTCCCGACTGTATAATTTTTGTAATGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCg < 1:116104/141‑1 (MQ=255) tcaaatatttcaaataGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTTTGTAATGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCggg < 1:107723/141‑1 (MQ=255) atatttcaaataGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTaca < 2:169905/141‑1 (MQ=255) atttcaaataGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTAGAGTaaaat > 2:84627/1‑137 (MQ=255) tttcaaataGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACtt > 2:105870/1‑141 (MQ=255) aaaTAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTTTGTAATGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAg < 2:48191/141‑1 (MQ=255) tAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAAGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGCGTACACTTTGCTATGAGCAAAGTGTACtttc > 2:71748/1‑140 (MQ=255) aGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTTTGTAATGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGAATACTCGGGTACACTTTGATATGAGCAAAGTGTACTTTGt > 2:223308/1‑141 (MQ=255) cATTACAAAAATTTTATAGAGTAGTGATTGTATAATTTCTATTGAGTTTAATGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGat < 2:242674/141‑1 (MQ=255) cATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTTTGTAATGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTTTTATTGat > 2:126882/1‑141 (MQ=255) cATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGat < 2:122812/141‑1 (MQ=255) aCAAAAATTGTATAGAGTATCTGCTGTATAATTTCGATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGataata > 1:183379/1‑141 (MQ=255) aCAAAAATGGTATAGAGTAGCGACAGTATAATTTTTGTAATGGTTAACGTCTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGAAAGACTACTCGCGTACACTTTGCTAAGAGCAAAGTCTCCTTTGTTATTGataata > 2:48237/1‑141 (MQ=255) aaaaatTTTTATGGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATAc < 2:34112/140‑1 (MQ=255) tataGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGcaca > 1:172518/1‑141 (MQ=255) ataGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACat > 1:139177/1‑141 (MQ=255) ataATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACAAATATATAAGTTTaacc > 2:127323/1‑139 (MQ=255) aTTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACATATATATAAGTTTAAACATa > 1:233119/1‑141 (MQ=255) tttCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACATATATATAAGTTTAAACATAc < 1:13579/141‑1 (MQ=255) cTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACATATATATAAGTTTAAACATACAGa < 1:207919/141‑1 (MQ=255) tGAGGTTAACGTTTATATGTAGTGATAGTAATTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACATATATATAAGTTTAAACATACAGATTTc < 1:93349/141‑1 (MQ=255) | TCTCTTGAAGTATAAAAAGGAAAAATTATTGAAAAAAATTCGACAAAAATATGTTTCAAATATTTCAAATAGGTTAAAATAATATTAGTTAACCATTACAAAAATTGTATAGAGTAGCGACTGTATAATTTCTATTGAGGTTAACGTTTATATGTAGTGATAGTAGTTAAAGTTCTCCCAAGGAAGACTACTCGGGTACACTTTGCTATGAGCAAAGTGTACTTTGTTATTGATAATACATTAGCACATATATATAAGTTTAAACATACAGATTTC > USA300TCH1516_ALE/643507‑643782 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |