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 | 94.0% | 296.1 / 7.0 | 83 | E13K (GAG→AAG) | cspC | stress protein, member of the CspA family |
Reads supporting (aligned to +/- strand): ref base C (2/3); new base T (30/48); total (32/51) | |||||||||||
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.98e-01 |
ATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAT > NC_000913/1907352‑1907444 | aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGt > 1:1339143/1‑50 (MQ=255) tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt > 1:2890803/1‑50 (MQ=255) cTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt > 1:2677379/1‑49 (MQ=255) tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACttgttg > 1:2741747/1‑49 (MQ=255) ttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGaa > 1:907295/1‑50 (MQ=255) tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGaa < 1:1212751/49‑1 (MQ=255) tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGa < 1:2259196/48‑1 (MQ=255) gctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc > 1:3251355/1‑49 (MQ=255) gctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc > 1:2577490/1‑49 (MQ=255) ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAcc < 1:1484912/49‑1 (MQ=255) ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAcc < 1:4393350/49‑1 (MQ=255) ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCa > 1:287370/1‑50 (MQ=255) ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCa < 1:200963/50‑1 (MQ=255) gcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt > 1:3876705/1‑50 (MQ=255) ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc < 1:1761507/48‑1 (MQ=255) ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt < 1:792710/49‑1 (MQ=255) ccATCAGCCGGAGTAATGAAACCAAAACCTTTAGACTTGTTGAACCACtt < 1:4299107/50‑1 (MQ=25) cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTa < 1:1075399/50‑1 (MQ=255) aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa < 1:1447526/50‑1 (MQ=255) aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa < 1:3010420/50‑1 (MQ=255) tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAc < 1:1061383/50‑1 (MQ=255) tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAAc < 1:360547/50‑1 (MQ=255) cAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc > 1:2312437/1‑50 (MQ=255) cAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc > 1:1566352/1‑50 (MQ=255) aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt < 1:3338402/50‑1 (MQ=255) gCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt > 1:1340953/1‑49 (MQ=255) gCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTg < 1:828389/50‑1 (MQ=255) ccGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa > 1:3242608/1‑50 (MQ=255) cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc > 1:3341131/1‑50 (MQ=255) cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc > 1:3507779/1‑50 (MQ=255) ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc > 1:2072664/1‑50 (MQ=255) ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc > 1:3031417/1‑50 (MQ=255) ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc > 1:961561/1‑50 (MQ=255) gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt > 1:1760363/1‑50 (MQ=255) aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt < 1:4237838/50‑1 (MQ=255) aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt < 1:3942406/50‑1 (MQ=255) aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt < 1:3878337/50‑1 (MQ=255) aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt < 1:3140107/50‑1 (MQ=255) gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt < 1:450368/50‑1 (MQ=255) gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt > 1:2911666/1‑49 (MQ=255) tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt > 1:2366306/1‑49 (MQ=255) tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTa < 1:3964633/50‑1 (MQ=255) tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTa < 1:1390306/50‑1 (MQ=255) tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTa < 1:2494616/50‑1 (MQ=255) aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa < 1:971441/49‑1 (MQ=255) aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa < 1:4324172/49‑1 (MQ=255) aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAAt > 1:439508/1‑50 (MQ=255) tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc > 1:3993033/1‑50 (MQ=255) tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc < 1:1505479/50‑1 (MQ=255) tGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATc < 1:43347/50‑1 (MQ=255) tGAAACCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc < 1:1746535/50‑1 (MQ=25) gAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt > 1:3121498/1‑50 (MQ=255) aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt < 1:4367012/50‑1 (MQ=255) aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt < 1:756871/49‑1 (MQ=255) aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt < 1:3701069/49‑1 (MQ=255) aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt < 1:669328/50‑1 (MQ=255) aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt < 1:302967/50‑1 (MQ=255) aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt < 1:2713667/50‑1 (MQ=255) ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc < 1:2611294/50‑1 (MQ=255) cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGcc < 1:4184220/50‑1 (MQ=255) cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc > 1:950778/1‑49 (MQ=255) aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa < 1:642130/50‑1 (MQ=255) aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa > 1:20096/1‑50 (MQ=255) aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa < 1:3759523/50‑1 (MQ=255) aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt > 1:3499097/1‑50 (MQ=255) aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt < 1:1081943/50‑1 (MQ=255) aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt > 1:3086892/1‑50 (MQ=255) aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt > 1:1444179/1‑50 (MQ=255) aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt < 1:2704058/50‑1 (MQ=255) aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt < 1:2615977/50‑1 (MQ=255) aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt < 1:497527/50‑1 (MQ=255) aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt < 1:4033996/50‑1 (MQ=255) ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg < 1:4016154/50‑1 (MQ=255) ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg < 1:596191/50‑1 (MQ=255) ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg < 1:653257/50‑1 (MQ=255) cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg > 1:4142620/1‑49 (MQ=255) tttAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa > 1:1073520/1‑49 (MQ=255) tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa < 1:4021133/49‑1 (MQ=255) tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa > 1:300970/1‑49 (MQ=255) tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa < 1:3184869/48‑1 (MQ=255) gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa < 1:4405525/49‑1 (MQ=255) gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt < 1:918675/50‑1 (MQ=255) gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt < 1:3670126/50‑1 (MQ=255) gACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt > 1:3625107/1‑50 (MQ=255) | ATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAT > NC_000913/1907352‑1907444 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |