Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,550,055 | 0 | A | C | 52.0% | 38.2 / 41.7 | 25 | V164G (GTG→GGG) | yfeX | porphyrinogen oxidase, cytoplasmic |
Reads supporting (aligned to +/- strand): ref base A (10/2); new base C (0/13); total (10/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.02e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-06 |
GACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAA‑CGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTG‑ATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCACGCTCTTCAACCCAACGGAAGCCGTGGATCTCTTCTTTCACTTCAATGC > NC_000913/2549915‑2550187 | gACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAA‑CGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCtt < 2:865454/150‑1 (MQ=255) ttCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAA‑CGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTG‑At < 2:273687/150‑1 (MQ=255) ctTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAA‑CGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTG‑ATAACCGCCACTTCGCGAcg > 2:230428/1‑151 (MQ=255) ccTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAA‑CGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTG‑ATAACCGCCACTTCGCGACGCGTCTCtt > 2:388687/1‑150 (MQ=255) cccGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAA‑CGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTG‑ATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGttt > 2:208413/1‑151 (MQ=255) gATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAA‑CGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTG‑ATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACaaa > 2:507975/1‑151 (MQ=255) acccctcccccGGTGGGGCTGCTTCAGTTTGTTTCCCCAC‑CGCTGAACAAACACATGGCTCCCCCCCCCACCCCCCCCTTTTTTG‑ATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTc < 1:849408/142‑1 (MQ=255) accGCTCACCCGGTTGAGCTCCTTCAGTTGTTGTCCCCAA‑CGCGGAAAAACCACATAGCTCCCGCCCGCACCCCCGCCTTTTTTG‑ATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTc < 2:654972/149‑1 (MQ=255) cGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAA‑CGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTG‑ATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAg > 2:669120/1‑151 (MQ=255) cggtggggctctttCAGGTTGTGTCCCCAC‑CGCTGAACAAACACAAAGCTGCCGCCCGCACCCCCGCCTTTTTTG‑ATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCACGct < 2:563683/139‑1 (MQ=255) gTTGCGCTCCTTCAGTTTGTTTCCCCAC‑CGCGGACCAACCACAACCCGCCCCCCCCCCCCCCCCCCGTCTTTA‑AAACCCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCACGctc < 1:564869/150‑1 (MQ=255) cctgctTCGGGTTGTTCCCCCAC‑CGCGGAACAAACACATAGCGCCCCCCCCCCCCCCCCCCGTTTTTG‑AAAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCACGCTCTTCAAc < 1:388685/150‑1 (MQ=255) gctTCAGGTTGTGTTCCCAA‑CGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTG‑ATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCACGCTCTTCAACCCa > 2:57987/1‑151 (MQ=255) ggTTGTGTTCCCAA‑CGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTG‑ATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCACGCTCTTCAACCCAACGGaa > 1:397864/1‑151 (MQ=255) ttttttccccAA‑CGCCGAACACACACATAGCCCCCCCCCTC‑CCCCCGCCCTTTTTGTAAAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCACGCTCTTCAACCCAACGGAAGc < 1:286879/144‑1 (MQ=255) tggtgtCCCCAC‑CGCGGCAAAAACACTTCGCCCCCCCCCCCCCCCCCGCCGTTTTTA‑AAACCCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCACGCTCTTCAACCCAACGGAAGc < 1:230427/149‑1 (MQ=255) ttttgtCCCAACCGCTGAACAACCACAACCCTGCCCCCCCCCCCCCCGCCTTTTTTA‑AAAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCACGCTCTTCAACCCAACGGAAGc < 1:59581/146‑1 (MQ=255) tCCCAA‑CGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTG‑ATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCACGCTCTTCAACCCAACGGAAGCCGTGGa > 1:370440/1‑151 (MQ=255) ccAC‑CGCGGGAAAACCACATGCCGCCCCCCCGCCCCCCCGCCGTTTTTG‑AAACCCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCACGCTCTTCAACCCAACGGAAGCCGTGGAtc < 1:384665/151‑1 (MQ=255) caaccgtgggACAACCACATCGCCCCCCCCCCCACCCCCGCGCTTTTTA‑ATACCCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCACGCTCTTCAACCCAACGGAAGCCGTGGAtc < 1:208413/143‑1 (MQ=255) accgcggaaCAAACACATAGCTCCCCCCCGCACCCCCCCCGTTTTTG‑ATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCACGCTCTTCAACCCAACGGAAGCCGTGGAtctc < 2:260001/143‑1 (MQ=255) acacATAGCTGCCGCCCGCATCCACGCCGTCTTTG‑ATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCACGCTCTTCAACCCAACGGAAGCCGTGGATCTCTTCTTTCACtt > 1:509785/1‑150 (MQ=255) cacaTAGCTGCCGCCCGCATCCACGCCGTCTTTG‑ATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCACGCTCTTCAACCCAACGGAAGCCGTGGATCTCTTCTTTCACTTCa > 1:373354/1‑151 (MQ=255) tagcgcccccccgcACCCCCGCCGTCTTTA‑ATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCACGCTCTTCAACCCAACGGAAGCCGTGGATCTCTTCTTTCACTTCAATGc < 2:540704/142‑1 (MQ=255) tacctcccgcccgcACCCCCGCCTTTTTTA‑ATACCCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCACGCTCTTCAACCCAACGGAAGCCGTGGATCTCTTCTTTCACTTCAATGc < 2:866896/145‑1 (MQ=255) | GACGTTCGTCGCCGTCGATCTCTTCGTTGGCCTCTTTGGTGCGCCCGATCACCATCTCCTGATCGTGAACGCTCATCCGGTTGAGCTGCTTCAGGTTGTGTTCCCAA‑CGCTGGACAAACACATAGCTGCCGCCCGCATCCACGCCGTCTTTG‑ATAACCGCCACTTCGCGACGCGTCTCTTCACCCGCCGGGTTTTCCGTACCGTCAACAAAGCCGCTCAGGTCACGCTCTTCAACCCAACGGAAGCCGTGGATCTCTTCTTTCACTTCAATGC > NC_000913/2549915‑2550187 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |