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 | 98.0% | 195.7 / ‑3.9 | 51 | E13K (GAG→AAG) | cspC | stress protein, member of the CspA family |
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (22/28); total (22/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.88e-01 |
ATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAT > NC_000913/1907352‑1907444 | aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGt > 1:1977663/1‑50 (MQ=255) tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt > 1:8740/1‑50 (MQ=255) tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt < 1:2606723/46‑1 (MQ=255) tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGaa < 1:1536179/49‑1 (MQ=255) ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCa > 1:1440202/1‑50 (MQ=255) ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCAc < 1:2971133/48‑1 (MQ=255) ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt > 1:2703121/1‑49 (MQ=255) ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt < 1:506797/49‑1 (MQ=255) ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACt < 1:215366/49‑1 (MQ=255) aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa > 1:66476/1‑50 (MQ=255) aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa > 1:2161081/1‑50 (MQ=255) aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa < 1:3614868/50‑1 (MQ=255) aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc < 1:1067705/49‑1 (MQ=255) aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAAcc < 1:2167537/49‑1 (MQ=255) aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt < 1:270511/50‑1 (MQ=255) aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt > 1:3656879/1‑50 (MQ=255) ccGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa < 1:2769586/50‑1 (MQ=255) ccGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa > 1:2852887/1‑50 (MQ=255) ccGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa > 1:3052800/1‑50 (MQ=255) ccGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa > 1:1521743/1‑50 (MQ=255) cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc > 1:25667/1‑50 (MQ=255) cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc > 1:2627735/1‑50 (MQ=255) ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc < 1:626654/50‑1 (MQ=255) ggAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAcc < 1:1274282/50‑1 (MQ=255) gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt < 1:3519735/50‑1 (MQ=255) aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt > 1:167890/1‑50 (MQ=255) gTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt < 1:3371732/50‑1 (MQ=255) aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAAt < 1:3036752/50‑1 (MQ=255) tGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATc < 1:3259142/50‑1 (MQ=255) aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt > 1:941270/1‑50 (MQ=255) aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt < 1:2653183/50‑1 (MQ=255) aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt < 1:1185959/49‑1 (MQ=255) aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt > 1:2460378/1‑50 (MQ=255) aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt > 1:1525338/1‑50 (MQ=255) gCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTg < 1:1374577/50‑1 (MQ=255) ccAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc < 1:3932144/50‑1 (MQ=255) cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGcc < 1:12652/50‑1 (MQ=255) aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt < 1:299061/50‑1 (MQ=255) aCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt < 1:2558996/50‑1 (MQ=255) ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt > 1:3161783/1‑49 (MQ=255) ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg < 1:3022770/50‑1 (MQ=255) ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg < 1:2286023/50‑1 (MQ=255) ccTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg < 1:1985952/50‑1 (MQ=255) cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg > 1:280625/1‑49 (MQ=255) tttAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa > 1:76057/1‑49 (MQ=255) ttAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa > 1:2803567/1‑49 (MQ=255) tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa < 1:2305855/49‑1 (MQ=255) tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa > 1:2349149/1‑49 (MQ=255) tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaa < 1:2566297/48‑1 (MQ=255) gACTTGTTGAACCACTTACCCTGACCTTTAATCTTTGCCATTTGAAAAAt < 1:2901923/50‑1 (MQ=38) gACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt > 1:1209025/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: ‑ |