Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,907,399 | C→T | E13K (GAG→AAG) | cspC ← | stress protein, member of the CspA family |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,907,399 | 0 | C | T | 84.1% | 594.1 / 94.8 | 214 | E13K (GAG→AAG) | cspC | stress protein, member of the CspA family |
Reads supporting (aligned to +/- strand): ref base C (18/16); new base T (86/94); total (104/110) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.09e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01 |
Only 100 of 216 total aligned reads displayed.
CATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCC > NC_000913/1907351‑1907447 | aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGt < 2:5718543/50‑1 (MQ=255) aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGt > 2:3491803/1‑50 (MQ=255) aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGt < 1:526153/50‑1 (MQ=255) tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACttgttg > 1:4680830/1‑49 (MQ=255) tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACttgttg > 2:3413513/1‑49 (MQ=255) tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt < 2:4567089/46‑1 (MQ=255) tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGaa < 2:389273/49‑1 (MQ=255) tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc > 2:756393/1‑50 (MQ=255) tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc > 1:4866345/1‑50 (MQ=255) tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc > 2:3035425/1‑50 (MQ=255) tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc > 1:2919723/1‑50 (MQ=255) ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAcc < 2:1969767/49‑1 (MQ=255) ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAcc < 2:2562585/49‑1 (MQ=255) ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCa > 1:2566055/1‑50 (MQ=255) ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCa > 1:728405/1‑50 (MQ=255) ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGGTGAACCa < 1:2396900/50‑1 (MQ=25) gcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt > 1:4727654/1‑50 (MQ=255) gcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt > 2:114558/1‑50 (MQ=255) ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAGCCAc < 2:4998572/48‑1 (MQ=37) ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt < 2:1457705/50‑1 (MQ=255) ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt < 1:5056317/50‑1 (MQ=255) ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt < 1:182191/50‑1 (MQ=255) ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt < 1:3714348/49‑1 (MQ=255) ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt < 2:719601/49‑1 (MQ=255) ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACtt < 2:3776372/50‑1 (MQ=255) ccATCAGCCAGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt > 1:5069752/1‑49 (MQ=38) cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTa < 1:389387/50‑1 (MQ=255) cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTa < 1:4310086/50‑1 (MQ=255) cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTa < 1:3940637/50‑1 (MQ=255) cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTa < 1:3294154/50‑1 (MQ=255) cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACtt > 2:2469410/1‑49 (MQ=255) aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa < 1:1405255/50‑1 (MQ=255) aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa < 1:52662/50‑1 (MQ=255) aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTaa > 2:4239656/1‑50 (MQ=255) aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt < 1:422914/50‑1 (MQ=255) aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt > 2:3225016/1‑50 (MQ=255) aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt < 1:3573756/50‑1 (MQ=255) gCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt > 2:2576492/1‑49 (MQ=255) cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc < 1:1268865/50‑1 (MQ=255) cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc > 2:1776791/1‑50 (MQ=255) cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc < 1:2356451/50‑1 (MQ=255) cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc > 1:3499213/1‑50 (MQ=255) cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc > 1:3784949/1‑50 (MQ=255) ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc > 1:3200024/1‑50 (MQ=255) ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc > 1:5484293/1‑50 (MQ=255) gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt > 1:2327191/1‑50 (MQ=255) gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt < 1:3835778/50‑1 (MQ=255) aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt < 2:5252440/50‑1 (MQ=255) gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt > 1:1610032/1‑49 (MQ=255) gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt > 1:4373540/1‑49 (MQ=255) tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt > 2:3920614/1‑49 (MQ=255) tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt > 2:5485586/1‑49 (MQ=255) tAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa < 2:472039/50‑1 (MQ=255) aaTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa > 2:1399875/1‑50 (MQ=255) aaTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTaa > 1:2434550/1‑50 (MQ=255) aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa < 1:609434/49‑1 (MQ=255) aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAAt > 2:3638320/1‑50 (MQ=255) tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc > 2:1301571/1‑50 (MQ=255) tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc > 1:4784753/1‑50 (MQ=255) tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc < 2:631631/50‑1 (MQ=255) tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc > 1:5459331/1‑50 (MQ=255) gAAGCCAAAACCTTTGGACTTGTTGAACCACTTAACCTGACCTTTAATCt > 2:4327760/1‑50 (MQ=25) aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt > 1:5435138/1‑50 (MQ=255) aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt < 2:1623337/50‑1 (MQ=255) aGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCttt > 1:1077764/1‑50 (MQ=255) gCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTg < 1:2206617/50‑1 (MQ=255) gCCAAAACCTTTAGACTTGTTAATCCACTTAACCTGACCTTTAATCTTTg < 1:1903241/50‑1 (MQ=25) ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc < 1:5798383/50‑1 (MQ=255) ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc > 2:1727837/1‑50 (MQ=255) ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc < 2:3207099/50‑1 (MQ=255) ccAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGc > 1:5364842/1‑50 (MQ=255) cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGcc < 2:2708364/50‑1 (MQ=255) cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc > 2:1005789/1‑49 (MQ=255) aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa < 2:1564293/50‑1 (MQ=255) aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa > 1:5729024/1‑50 (MQ=255) aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt > 2:5220156/1‑50 (MQ=255) aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt > 2:4493840/1‑50 (MQ=255) aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt > 2:4861696/1‑50 (MQ=255) aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt > 2:3620118/1‑50 (MQ=255) aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt < 2:1342641/50‑1 (MQ=255) aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt < 1:3632634/50‑1 (MQ=255) aCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt > 2:5287760/1‑50 (MQ=255) ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt > 2:1684801/1‑49 (MQ=255) ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg < 1:3980452/50‑1 (MQ=255) cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa < 1:2732280/50‑1 (MQ=255) cTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa < 1:173308/50‑1 (MQ=255) tttAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa > 1:64550/1‑49 (MQ=255) ttAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa > 2:2853062/1‑49 (MQ=255) ttAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa > 2:2260868/1‑49 (MQ=255) tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa > 2:4382785/1‑49 (MQ=255) tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa > 1:191652/1‑49 (MQ=255) tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa > 2:5200288/1‑49 (MQ=255) tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa < 1:3724790/48‑1 (MQ=255) tAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa < 2:2039338/49‑1 (MQ=255) aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa > 2:868283/1‑50 (MQ=255) aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa > 1:1870253/1‑50 (MQ=255) gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa < 2:3092790/49‑1 (MQ=255) gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt < 2:3269078/50‑1 (MQ=255) gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt > 1:843134/1‑50 (MQ=255) gACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt < 1:93587/50‑1 (MQ=255) | CATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCC > NC_000913/1907351‑1907447 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |