Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,479,751 | 1 | . | C | 21.1% | 131.9 / 49.0 | 52 | coding (857/1806 nt) | kefB | potassium:proton antiporter |
Reads supporting (aligned to +/- strand): ref base . (20/21); new base C (6/5); total (26/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01 |
ATACAATCGCGCCAGCAGATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATGGCGATATTCA > NC_000913/3479660‑3479849 | aTACAATCGCGCCAGCAGATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAgg < 1:1260797/100‑1 (MQ=255) tACAATCGCGCCAGCAGATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGt < 2:942718/100‑1 (MQ=255) aTCGCGCCAGCAGATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAg > 2:518174/1‑100 (MQ=255) tCGCGCCAGCAGATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAg < 2:1032005/99‑1 (MQ=255) cgcCAGCAGATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGAc < 2:438407/100‑1 (MQ=255) cagcagATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACa < 2:790410/98‑1 (MQ=255) agcagATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCCGAGGTTGAGTGACATg > 2:220572/1‑100 (MQ=255) cagATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCg < 2:1127589/100‑1 (MQ=255) gATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGAc > 1:643721/1‑100 (MQ=255) aTACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCCGAGGTTGAGTGACATGCCGAc > 1:77623/1‑100 (MQ=255) cAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCCGAGGTTGAGTGACATGCCGACaga < 2:1081602/100‑1 (MQ=255) cAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACagag > 1:585240/1‑100 (MQ=255) aGCACTAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAgaga < 2:173563/100‑1 (MQ=255) agaATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCCGAGGTTGAGTGACATGCCGACAGAGATAAAg < 2:786308/100‑1 (MQ=255) ttttCACCGCCACCAGTACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACaa > 2:436138/1‑100 (MQ=255) ttttCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACaa > 1:1254577/1‑100 (MQ=255) ttttCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACaa > 2:968098/1‑100 (MQ=255) tCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAAcc > 1:1084709/1‑100 (MQ=255) aCCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGa > 2:680748/1‑100 (MQ=255) ccaccaGCACAACCACGCTTATCACTACCCACAACATATCGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGagca < 2:1001036/99‑1 (MQ=255) ccaccaGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGagca > 2:765871/1‑100 (MQ=255) ccaccaGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGagcag < 1:1064356/100‑1 (MQ=255) ccaGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCaa > 2:764636/1‑100 (MQ=255) cacaACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCt > 1:1204449/1‑100 (MQ=255) cacaACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCt < 2:50581/100‑1 (MQ=255) cacaACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCt < 2:50569/100‑1 (MQ=255) acaACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCtt > 1:887012/1‑100 (MQ=255) caACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCttt < 1:372637/100‑1 (MQ=255) aaCCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTg < 2:1254566/100‑1 (MQ=255) aCCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTg < 1:220574/100‑1 (MQ=255) aCCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGa > 2:368141/1‑100 (MQ=255) cACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGaa < 2:1030028/100‑1 (MQ=255) cACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAg < 2:779418/100‑1 (MQ=255) aCGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAgg < 2:443284/100‑1 (MQ=255) cTTATCACTACCCACAACAGATGGGTATAAAGCACCCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGa < 2:768996/100‑1 (MQ=255) cTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGAt < 2:160135/100‑1 (MQ=255) aTCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATc < 1:1203723/98‑1 (MQ=255) aCCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCg > 2:791948/1‑100 (MQ=255) cccACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGt > 1:622139/1‑100 (MQ=255) caacaGATGGGTATCAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTcc < 2:887005/100‑1 (MQ=255) caacaGATGGGTATAAAGCA‑CCCCGATGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTcc < 2:286670/100‑1 (MQ=255) caGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGt < 2:101977/100‑1 (MQ=255) caGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGt > 2:332146/1‑100 (MQ=255) caGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGt > 2:531977/1‑100 (MQ=255) gATGGGTATAAAGCACCCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGtt > 1:173544/1‑100 (MQ=255) gATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTc < 2:456028/100‑1 (MQ=255) tGGGTATAAAGCACCCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCa > 1:755618/1‑100 (MQ=255) ggTATAAAGCACCCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATg > 2:841477/1‑100 (MQ=255) ggTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATg > 2:309021/1‑99 (MQ=255) ataAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCGATGAAGGGATCGATAGCCGTCTCCAGTTCCTGGCg > 2:189182/1‑99 (MQ=255) taAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATGGCGa > 1:171697/1‑99 (MQ=255) cA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATGGCGATATTCa > 1:455229/1‑100 (MQ=255) | ATACAATCGCGCCAGCAGATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCA‑CCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATGGCGATATTCA > NC_000913/3479660‑3479849 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |