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 | 97.5% | 618.3 / 3.3 | 160 | E13K (GAG→AAG) | cspC | stress protein, member of the CspA family |
Reads supporting (aligned to +/- strand): ref base C (3/1); new base T (82/74); total (85/75) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.23e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01 |
Only 100 of 164 total aligned reads displayed.
ATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATT > NC_000913/1907352‑1907445 | aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGt < 1:5209960/50‑1 (MQ=255) aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGt < 1:4170583/50‑1 (MQ=255) aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGt < 1:3421949/50‑1 (MQ=255) cTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt > 2:5038494/1‑49 (MQ=255) tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACttgttg > 2:4949528/1‑49 (MQ=255) tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACttgttg > 2:2262656/1‑49 (MQ=255) ttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGaa > 1:198452/1‑50 (MQ=255) tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACttgttg < 2:3096324/47‑1 (MQ=255) tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt < 2:1614450/46‑1 (MQ=255) tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGa < 2:5329250/48‑1 (MQ=255) tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGa < 1:2704388/48‑1 (MQ=255) tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGa < 2:1715653/48‑1 (MQ=255) tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc > 2:4463978/1‑50 (MQ=255) tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc > 2:4335750/1‑50 (MQ=255) gctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc > 1:4618913/1‑49 (MQ=255) tgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc > 2:616892/1‑50 (MQ=255) tgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc > 1:5357409/1‑50 (MQ=255) tgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc > 1:3365831/1‑50 (MQ=255) gcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt > 1:2729549/1‑50 (MQ=255) ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc < 1:2831582/48‑1 (MQ=255) ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt < 1:4462329/50‑1 (MQ=255) ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt < 1:5293757/49‑1 (MQ=255) cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTa < 1:3197695/50‑1 (MQ=255) aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa > 1:4091427/1‑50 (MQ=255) aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa > 1:1691451/1‑50 (MQ=255) aGCCGGAGTAATGAAGCCACAACCTTTAGACTTGTTGAACCACTTAACCt > 2:1186650/1‑50 (MQ=37) aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt < 2:2062563/50‑1 (MQ=255) aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt > 1:706196/1‑50 (MQ=255) aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt > 1:1929132/1‑50 (MQ=255) gCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTg < 2:1212808/50‑1 (MQ=255) ccGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa > 2:2708231/1‑50 (MQ=255) cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc < 2:1760151/50‑1 (MQ=255) ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc > 1:1024699/1‑50 (MQ=255) ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc > 2:381071/1‑50 (MQ=255) ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc < 1:209997/50‑1 (MQ=255) gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt > 1:2907163/1‑50 (MQ=255) gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt > 1:4553450/1‑50 (MQ=255) aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt > 1:3579521/1‑50 (MQ=255) aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt > 1:3957936/1‑50 (MQ=255) aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt < 1:4811346/50‑1 (MQ=255) aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt > 2:4756700/1‑50 (MQ=255) gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt < 1:2720744/50‑1 (MQ=255) gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt < 1:5518318/50‑1 (MQ=255) gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt < 2:5229378/50‑1 (MQ=255) tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt > 1:4371589/1‑49 (MQ=255) tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt > 1:2356016/1‑49 (MQ=255) tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTa < 2:4005360/50‑1 (MQ=255) aaTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa > 1:1087657/1‑50 (MQ=255) aaTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa > 1:5032911/1‑50 (MQ=255) aaTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa > 1:4986272/1‑50 (MQ=255) aaTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa > 1:1916397/1‑50 (MQ=255) aaTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa > 1:4606742/1‑50 (MQ=255) aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa < 2:1361302/49‑1 (MQ=255) aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAAt > 2:4630912/1‑50 (MQ=255) aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAAt < 1:3069632/50‑1 (MQ=255) aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAAt > 2:3591034/1‑50 (MQ=255) tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc > 2:1583094/1‑50 (MQ=255) tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc > 1:4635416/1‑50 (MQ=255) gAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt > 1:1767113/1‑50 (MQ=255) gAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt > 2:297853/1‑50 (MQ=255) aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt > 2:3852957/1‑50 (MQ=255) aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt < 1:3638453/50‑1 (MQ=255) aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt > 2:1264292/1‑50 (MQ=255) aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt > 2:2486563/1‑50 (MQ=255) aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt < 1:991141/50‑1 (MQ=255) aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt > 2:2858014/1‑50 (MQ=255) gCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTg < 2:219010/50‑1 (MQ=255) gCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTg < 1:250520/50‑1 (MQ=255) ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc > 1:76521/1‑50 (MQ=255) cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc > 1:5406311/1‑49 (MQ=255) aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa > 1:2119326/1‑50 (MQ=255) aaaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa < 2:745811/50‑1 (MQ=255) aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt > 2:195687/1‑50 (MQ=255) aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt > 1:846691/1‑50 (MQ=255) aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt > 1:1844087/1‑50 (MQ=255) aCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt < 1:4484784/50‑1 (MQ=255) aCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt < 1:1534623/50‑1 (MQ=255) ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt > 2:4263132/1‑49 (MQ=255) ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt > 2:4054704/1‑49 (MQ=255) ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg < 1:2829395/50‑1 (MQ=255) ccTTTAGACTTGTTGAACCACTTAACCTGACCTTAAATCTTTGCCATTTg < 1:3098813/50‑1 (MQ=38) cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg > 1:5566972/1‑49 (MQ=255) cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg > 2:1512092/1‑49 (MQ=255) tttAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa > 2:2478122/1‑49 (MQ=255) tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGac < 2:1047918/48‑2 (MQ=255) tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa > 1:2339327/1‑49 (MQ=255) tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa < 1:610302/48‑1 (MQ=255) aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa > 1:2439846/1‑50 (MQ=255) aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa > 1:1158572/1‑50 (MQ=255) aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa > 1:1960303/1‑50 (MQ=255) aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa > 1:1691921/1‑50 (MQ=255) aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaa < 1:5187786/49‑1 (MQ=255) aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaa < 1:3018589/49‑1 (MQ=255) gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa < 2:2570490/49‑1 (MQ=255) gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaaa < 2:1565994/49‑1 (MQ=255) gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt < 1:1349436/50‑1 (MQ=255) gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt < 2:2607275/50‑1 (MQ=255) gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt > 2:913525/1‑50 (MQ=255) aCTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAtt > 2:1454880/1‑50 (MQ=255) aCTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAtt < 1:2941081/50‑1 (MQ=255) | ATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATT > NC_000913/1907352‑1907445 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |