Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | USA300TCH1516_ALE | 2,072,028 | Δ1 bp | 5.8% | intergenic (+523/+50) | USA300HOU_RS10445 → / ← USA300HOU_RS10450 | hypothetical protein/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 2,072,028 | 0 | T | . | 5.8% | 233.8 / 7.4 | 69 | intergenic (+523/+50) | USA300HOU_RS10445/USA300HOU_RS10450 | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (37/27); new base . (2/2); total (40/29) | |||||||||||
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.93e-01 |
AGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAG > USA300TCH1516_ALE/2071980‑2072118 | aGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTa < 1:1384085/91‑1 (MQ=31) gTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTaa < 1:2010221/91‑1 (MQ=32) tATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTaaa > 1:2077910/1‑91 (MQ=33) tATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTaaa > 1:602115/1‑91 (MQ=33) aTTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAAt > 1:2298306/1‑91 (MQ=34) ttGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATg > 1:1412371/1‑91 (MQ=35) tGTATAAATAGAGAACAGCAGTAAGAA‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGt < 1:3030107/91‑1 (MQ=255) gTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGtt > 1:2181391/1‑91 (MQ=38) gTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGtt > 1:1293517/1‑91 (MQ=38) tataAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGtta > 1:4128324/1‑91 (MQ=39) ataAATAGAGAACAGCAGTAAGATAATTTTCAATTA‑GAAAA‑TATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGttat < 1:1269775/91‑1 (MQ=11) ataAATAGAGAACAGCAGTAAGATAATTTTC‑AATTAGAAAA‑TATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGttat < 1:1269383/91‑1 (MQ=11) ataAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGttat < 1:1777900/91‑1 (MQ=40) ataAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGttat < 1:1778918/91‑1 (MQ=40) aaaTAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGttatta > 1:3522207/1‑91 (MQ=255) aaaTAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGttatta > 1:3992393/1‑91 (MQ=255) aaTAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAg < 1:2132211/91‑1 (MQ=255) aTAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGt < 1:1762033/91‑1 (MQ=255) aTAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGt < 1:3768336/91‑1 (MQ=255) aTAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGt < 1:3400672/91‑1 (MQ=255) tAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGtt < 1:222444/91‑1 (MQ=255) tAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGtt > 1:685406/1‑91 (MQ=255) tAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGtt < 1:1387606/91‑1 (MQ=255) gagaACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTAt > 1:361944/1‑91 (MQ=255) agaACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTAtt > 1:3456973/1‑91 (MQ=255) agaACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTAtt > 1:2995190/1‑91 (MQ=255) agaACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTAtt < 1:2725666/91‑1 (MQ=255) gaACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTAttt > 1:4084601/1‑91 (MQ=255) aaCAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTAtttt < 1:2748919/91‑1 (MQ=255) aCAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCAGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTAttttt > 1:1562629/1‑91 (MQ=255) cagcagTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTa > 1:2126812/1‑91 (MQ=255) agcagTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTAc > 1:2775860/1‑91 (MQ=255) agcagTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTAc > 1:2775970/1‑91 (MQ=255) gcagTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACt > 1:501990/1‑91 (MQ=255) gcagTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACt > 1:4148700/1‑91 (MQ=255) gcagTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACt < 1:4108091/91‑1 (MQ=255) cagTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACtt > 1:403928/1‑91 (MQ=255) cagTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACtt < 1:2016248/91‑1 (MQ=255) cagCAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTACTTATTTTTACtt < 1:3206976/91‑1 (MQ=255) agTAAGATAATTTTC‑AATTAGAAAA‑TATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCAAAATGTTATTAGTTATTTTTACttt > 1:2810516/1‑91 (MQ=17) agTAAGATAATTTTC‑AATT‑AGAAAATATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACttt > 1:3804176/1‑91 (MQ=18) agTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACttt > 1:3835218/1‑91 (MQ=255) agTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACttt > 1:597378/1‑91 (MQ=255) agTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACttt < 1:1064889/91‑1 (MQ=255) aaGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCaa > 1:656141/1‑91 (MQ=255) aaGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCaa > 1:1692801/1‑91 (MQ=255) aaGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCaa > 1:2135407/1‑91 (MQ=255) aGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAat < 1:3794742/91‑1 (MQ=255) gAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAata > 1:3081622/1‑91 (MQ=255) at‑atTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAatat < 1:1713566/91‑1 (MQ=255) atTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCa < 1:2232177/91‑1 (MQ=255) ttttCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAg > 1:1458642/1‑91 (MQ=255) ttttCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAg < 1:446107/91‑1 (MQ=255) ttttCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAg > 1:1458509/1‑91 (MQ=255) tttCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGa > 1:2754804/1‑91 (MQ=255) cTAATT‑GAAAATTATGTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAaat > 1:3532451/1‑91 (MQ=255) cTAATT‑GAAAATTATGTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGCTATTAGTTATTTTTACTTTCAATATCAGAaat > 1:3530433/1‑91 (MQ=255) cTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGGTATTAGTTATTTTTACTTTCAATATCAGAaat > 1:3532214/1‑91 (MQ=255) cTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGCTATTAGTTATTTTTACTTTCAATATCAGAaat > 1:3531628/1‑91 (MQ=255) tAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAaata < 1:2658946/91‑1 (MQ=255) tAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAaata < 1:4038014/91‑1 (MQ=255) aaTT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAaataa > 1:2390545/1‑91 (MQ=255) aTT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAaataat > 1:173530/1‑91 (MQ=255) t‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAAttt > 1:1061275/1‑91 (MQ=255) agAAAA‑TATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTc > 1:2802219/2‑91 (MQ=255) gAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTc < 1:4109009/91‑1 (MQ=255) aaaTTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCtt < 1:2285817/91‑1 (MQ=255) aaaTTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCtt < 1:775877/91‑1 (MQ=255) aaTTAGCTTACTGCTGTTTTTTAGGCATTTATGTCCCAGGCTCTAAATGTTATGAGTTATTTTTACTTTCAATATCAGAAATAATTTCttt < 1:3862702/91‑1 (MQ=37) ttATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAg < 1:980491/91‑1 (MQ=255) ttATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAg > 1:3360151/1‑91 (MQ=255) ttATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAg < 1:2614543/91‑1 (MQ=255) ttATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAg < 1:3496108/91‑1 (MQ=255) ttATCTTACTGCTGTTTTTTAGGGAGTTATGTCCCAGCCTCTAGATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAg < 1:3224036/91‑1 (MQ=255) | AGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAG > USA300TCH1516_ALE/2071980‑2072118 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |