| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,029,274 | 0 | C | T | 25.0% | 132.9 / 12.6 | 60 | K31K (AAG→AAA) | yedI | conserved inner membrane protein YedI |
| Reads supporting (aligned to +/- strand): ref base C (14/31); new base T (15/0); total (29/31) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.46e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCG > NC_000913/2029225‑2029322 | cTTGTTGCGCATTGAGGGATAAGTCATCCCCTAATACACCGGCGGTTttc > 1:2763179/1‑50 (MQ=255)cTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTttc < 1:5149232/50‑1 (MQ=255)cTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTttc > 1:7160836/1‑50 (MQ=255)cTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTttc > 1:8443816/1‑50 (MQ=255) tgttgCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTtttttt > 1:5024903/1‑50 (MQ=255) tgttgCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTttctt > 1:3200389/1‑50 (MQ=255) gttgCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTttctt < 1:323663/49‑1 (MQ=255) gttgCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTttct < 1:8518311/48‑1 (MQ=255) gttgCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTttct < 1:6392020/48‑1 (MQ=255) ttgCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTttcttc < 1:3636905/49‑1 (MQ=255) ttgCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTttcttc < 1:1571003/49‑1 (MQ=255) ttgCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTttcttc < 1:2138220/49‑1 (MQ=255) ttgCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTcg > 1:8591275/1‑50 (MQ=255) tgCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTtttttctc > 1:4580421/1‑48 (MQ=255) gcgcATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTcgc < 1:3357711/49‑1 (MQ=255) cgcATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTcgccc > 1:5846180/1‑49 (MQ=255) cgcATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTcgcc < 1:1278622/49‑1 (MQ=255) gcATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTctcccc > 1:8124057/1‑48 (MQ=255) aTTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTcgccgc > 1:1141624/1‑48 (MQ=255) aTTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCa < 1:2250661/50‑1 (MQ=255) aTTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCa < 1:1500869/50‑1 (MQ=255) tGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTTCCCCCCAgg > 1:3428899/1‑49 (MQ=255) tGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTTGCCGCCCGt > 1:8099199/1‑50 (MQ=255) tGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGTCGCCAGt > 1:2257558/1‑50 (MQ=255) gAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAGt < 1:221936/49‑1 (MQ=255) aGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAGttt < 1:1077541/50‑1 (MQ=255) aGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAGtt < 1:1610761/49‑1 (MQ=255) aGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAGtt < 1:460269/49‑1 (MQ=255) gCGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTCGCCGCCCGTTTg > 1:2302580/1‑50 (MQ=39) cGATAAGTCATCCCCTAATACACCGGCGGTTTTTTTCCCCGCCAGTTTGc > 1:6717995/1‑50 (MQ=255) aTAAGTCATCCCCTAATACACCGGCGGTTTTTTTCGCCGCCCGTTTTccc > 1:6401613/1‑50 (MQ=37) aTAAGTCATCCCCTAATACACCGGCGGTTTTTTTCGCCGCCAGTTTGccc > 1:3188185/1‑50 (MQ=255) tAAGTCATCCCCTAATACACCGGCGGTTTTTTTCGCCCCCCGTTTTcccc > 1:6399918/1‑49 (MQ=255) aaGTCATCCCCTAATACACCGGCGGATTTCTTCGCCGCCCGTTTTcccct > 1:2503763/1‑48 (MQ=37) aGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAGTTTGCCCAt < 1:2975981/49‑1 (MQ=255) aGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAGTTTGCCCATc < 1:3408209/50‑1 (MQ=255) aGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAGTTTGCCCATc < 1:1093730/50‑1 (MQ=255) gTCATCCCCTAATACACCGGCGGTTTTTTTCGCCGCCAGTTTGCCCCTCa > 1:576705/1‑50 (MQ=255) gTCATCCCCTAATACACCGGCGGTTTTTTTCGCCCCCCGCTTGCCCATCa > 1:2078706/1‑50 (MQ=255) gTCATCCCCTAATACACCGGCGGTTTTCTTTGCCCCCCGGTTGCCCCTCa > 1:3767712/1‑50 (MQ=255) gTCATCCCCTAATACACCGGCGGTTTTCTTTGCCCCCAGTTTTCCCCTCa > 1:164687/1‑50 (MQ=255) aTCCCCTAATACACCGGCGGTTTTTTTCGCCGCCAGTTTGCCCAACACg > 1:4967846/1‑49 (MQ=38) aTCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAGTTTGCCCATCACgg < 1:2720481/50‑1 (MQ=255) tCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAGTTTGCCCATCACGga < 1:3904484/50‑1 (MQ=255) cccTAATACACCGGCGGTTTTTTTCTCCGCCAGTTTGCCCCTCAcc > 1:1342903/1‑45 (MQ=37) ccTAATACACCGGCGGTTTTTTTCGCCCCCAGGTTTCCCCTCACGgaga > 1:3645463/1‑49 (MQ=255) ccTAATACACCGGCGGTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGat < 1:6368107/50‑1 (MQ=255) cTAATACACCGGCGGTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGata < 1:1441119/50‑1 (MQ=255) aaTACACCGGCGGTTTTTTTCGCCGCCCATTTGCCCCTCACGGAGATATc > 1:6775869/1‑50 (MQ=255) aTACACCGGCGGTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATc < 1:6729097/49‑1 (MQ=255) aTACACCGGCGGTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCg < 1:1073427/50‑1 (MQ=255) aTACACCGGCGGTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCg < 1:4569605/50‑1 (MQ=255) tACACCGGCGGTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGt < 1:6336378/50‑1 (MQ=255) tACACCGGCGGTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGt < 1:7643283/50‑1 (MQ=255) cacCGGCGGTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTcc < 1:319427/50‑1 (MQ=255) cacCGGCGGTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTcc < 1:2627328/50‑1 (MQ=255) acCGGCGGTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTCCa < 1:7358714/50‑1 (MQ=255) acCGGCGGTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTCCa < 1:896223/50‑1 (MQ=255) acCGGCGGTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTCCa < 1:2714613/50‑1 (MQ=255) ccGGCGGTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTCCag > 1:3407323/1‑50 (MQ=255) cggcggTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTCCagt < 1:5443477/50‑1 (MQ=255) cggcggTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTCCagt < 1:3549487/50‑1 (MQ=255) ggcggTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTCCagta < 1:7419780/50‑1 (MQ=255) ggcggTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTCCag > 1:5204734/1‑48 (MQ=255) ggcggTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTCCag > 1:7255097/1‑48 (MQ=255) gcggTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTCCagtag < 1:1907030/50‑1 (MQ=255) gcggTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTCCagtag < 1:6303950/50‑1 (MQ=255) cggTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTCCagtagt < 1:1420087/50‑1 (MQ=255) ggTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTCCAGTAgtg > 1:1552339/1‑50 (MQ=255) ggTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTCCAGTAgtg > 1:6520531/1‑50 (MQ=255) ttttCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTCCAGTAgtgcg > 1:6294745/1‑48 (MQ=255) tttCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTCCAGTAGTGTgg > 1:4967598/1‑50 (MQ=255) tttCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTCCAGTAGTGTgg > 1:2501820/1‑50 (MQ=255) cttcGCCGCCAGTTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCg < 1:4991663/49‑1 (MQ=255) cttcGCCGCCAGTTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCg < 1:1897378/49‑1 (MQ=255) | CTTGTTGCGCATTGAGCGATAAGTCATCCCCTAATACACCGGCGGTTTTCTTCGCCGCCAGTTTGCCCATCACGGAGATATCGTCCAGTAGTGTGGCG > NC_000913/2029225‑2029322 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |