Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913_3_bme_pgi | 4,296,382 | +GC | intergenic (+587/+55) | gltP → / ← yjcO | glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913_3_bme_pgi | 4,296,381 | 1 | . | C | 93.9% | 79.9 / 0.4 | 33 | intergenic (+586/+56) | gltP/yjcO | glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein |
Reads supporting (aligned to +/- strand): ref base . (0/2); new base C (13/18); total (13/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.08e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01 | |||||||||||
* | NC_000913_3_bme_pgi | 4,296,381 | 2 | . | G | 93.8% | 82.1 / 0.2 | 33 | intergenic (+586/+56) | gltP/yjcO | glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein |
Reads supporting (aligned to +/- strand): ref base . (0/2); new base G (13/17); total (13/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.02e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01 |
GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTG > NC_000913_3_bme_pgi/4296295‑4296511 || gCGACGCTGAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTAATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCTACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc < 1:107073/149‑1 (MQ=1) gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACCTCAATGCCTGATGCGACGCTTGCCCCGTCTTATCAATCCTATCTTAACCg > 1:29719/1‑149 (MQ=11) aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGt < 1:231120/149‑1 (MQ=12) acgcacggcgGCTTATCAGGGCTACGGCAGGCAGCGCAATAGCCTGATTTAGGGTGAATTTGTAGGTCGGATAAGGCGTTTCCG‑‑CCGCATACGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg < 2:294909/142‑1 (MQ=11) aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa > 2:697886/1‑149 (MQ=255) tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACCCTTGCCGCGTCATATCAGACCTATCTTAACCGTTAGTTAattatt > 1:256840/1‑149 (MQ=11) tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc > 1:71305/1‑149 (MQ=255) gCCTACGCCAGACAGCGCAAATGCCGGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt < 1:214410/149‑1 (MQ=17) cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCCTca > 1:312961/1‑149 (MQ=255) cTACGCCAGACAGCGCAATAGCCTAATTGAGCGTGAATTTGTAAGGCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca < 2:310547/149‑1 (MQ=17) cAGACGGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg < 1:270562/149‑1 (MQ=255) aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg < 1:529452/149‑1 (MQ=255) aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCTGGCCGGTGTGAACCTTTGGTTAATTTTTCAAACTCTTCGCACCCGGTa > 1:220789/1‑149 (MQ=14) gCCTAATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCGGCGTCTTAGCAGGCTCATCCTAACCCTTGGTAAATTTTTCTAACTCTTCACACCCGGTATCAAACCCTTc < 1:663920/149‑1 (MQ=17) gCCCGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc < 1:37073/149‑1 (MQ=255) ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc > 2:64273/1‑149 (MQ=255) tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCAt > 2:515777/1‑149 (MQ=255) gagttAGTGTGATTTTGTAGTGCGGATAAGGCGTGT‑‑ACGCCGCATCAGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTCACCGTTGGTTAATTTTTCAAACGCTTCACACTCGGTATCAAACCCTTCCATa < 2:736930/146‑1 (MQ=18) gATTGAGCGGTATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa < 1:466556/149‑1 (MQ=21) ttttAGGGGGATTTGTTAGGGCGGTTAATGCGGTC‑‑ACTCCGGATCCGACATCAATGCCTGATGGGACGCTTGCCGCGTCGTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc < 1:750390/148‑1 (MQ=255) tttAGCGTGATTTTGGAGGCCGGATAAGGCGTTTCCG‑‑CCGCCTACGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACa < 1:206013/149‑1 (MQ=21) gATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa > 1:626622/1‑108 (MQ=255) gATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa < 2:626622/108‑1 (MQ=255) gATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt < 1:172752/149‑1 (MQ=255) ttttGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTc < 1:329149/149‑1 (MQ=255) tAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCTACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAACCCCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCa < 2:203802/149‑1 (MQ=25) tAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCa < 1:490438/149‑1 (MQ=255) aGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGGCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAg < 1:768942/149‑1 (MQ=255) ggTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTTCATTCAGCTCAGGTTTAACCaat > 2:284438/1‑147 (MQ=25) gccGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAGTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCCGta > 2:723763/3‑148 (MQ=25) gTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTg > 1:279355/1‑149 (MQ=255) aaGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc > 1:477785/1‑88 (MQ=17) aaGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc < 2:477785/88‑1 (MQ=17) || GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTG > NC_000913_3_bme_pgi/4296295‑4296511 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |