Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,503,834 | A→G | 9.5% | intergenic (+39/‑59) | xanP → / → tusD | xanthine permease/sulfur transfer protein complex subunit TusD |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,503,834 | 0 | A | G | 9.5% | 118.9 / 3.4 | 42 | intergenic (+39/‑59) | xanP/tusD | xanthine permease/sulfur transfer protein complex subunit TusD |
Reads supporting (aligned to +/- strand): ref base A (20/18); new base G (2/2); total (22/20) | |||||||||||
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.26e-01 |
GGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGC > NC_002947/4503699‑4503975 | ggcgATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgggcgcggg > 2:118672/1‑150 (MQ=255) gATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCtt > 1:1472328/1‑150 (MQ=255) cAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGc > 1:1490952/1‑150 (MQ=255) tgtgCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCttt < 2:473984/150‑1 (MQ=255) tgtgCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCttt > 1:67803/1‑150 (MQ=255) cTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCa < 2:67803/150‑1 (MQ=255) tgtgGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTc < 1:118672/150‑1 (MQ=255) tACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATgg > 2:888487/1‑136 (MQ=255) tACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATgg < 1:888487/136‑1 (MQ=255) gggggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAg > 2:1126268/1‑150 (MQ=255) cTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgggcgcgg > 2:842276/1‑91 (MQ=255) cTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgggcgcgg < 1:842276/91‑1 (MQ=255) tGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTAtt > 1:255396/1‑149 (MQ=255) aaCTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGa < 2:107766/137‑1 (MQ=255) aaCTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGa > 1:107766/1‑137 (MQ=255) aCTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGTGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATcgcg < 2:879106/150‑1 (MQ=255) cTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCgg > 2:727193/1‑150 (MQ=255) tGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGgt < 2:1439429/150‑1 (MQ=255) gTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATc < 1:1346074/143‑1 (MQ=255) gTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATc > 2:1346074/1‑143 (MQ=255) gTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCGTGAGACTGAATGAGAGTAGCGCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGgtgt > 2:17381/1‑150 (MQ=255) ttCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTTCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGtt > 1:1541401/1‑150 (MQ=255) ttGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCtttttt > 1:1187588/1‑87 (MQ=255) ttGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCtttttt < 2:1187588/87‑1 (MQ=255) gCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTccc < 1:1172956/150‑1 (MQ=255) cGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCAt < 2:1476341/93‑1 (MQ=255) cGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCAt > 1:1476341/1‑93 (MQ=255) gAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGc < 1:657736/150‑1 (MQ=255) ggAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATgcgc < 2:1472328/150‑1 (MQ=255) ttGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGTGCAAgcgggcgcgg > 2:1600589/1‑56 (MQ=21) ttGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGTGCAAgcgggcgcgg < 1:1600589/56‑1 (MQ=21) ttGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGTGCAAgcgggcgc > 2:1049626/1‑54 (MQ=21) ttGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGTGCAAgcgggcgc < 1:1049626/54‑1 (MQ=21) ttGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGtt > 2:1445589/1‑97 (MQ=255) ttGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGtt < 1:1445589/97‑1 (MQ=255) aaTGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCtt > 1:1029221/1‑150 (MQ=255) agTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGa < 2:1541401/150‑1 (MQ=255) agTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGa > 2:1132885/1‑150 (MQ=255) ccgcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGt < 1:727193/150‑1 (MQ=255) ccgcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGt > 1:50234/1‑150 (MQ=255) cggcatacacaAGCGTGCGCTGGCTGTGTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGc > 2:1257558/8‑150 (MQ=255) ccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGgt > 2:482207/1‑87 (MQ=255) ccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGgt < 1:482207/87‑1 (MQ=255) | GGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGC > NC_002947/4503699‑4503975 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |