Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,398,903 | C→G | 15.1% | intergenic (+123/‑119) | aroQ‑III → / → PP_3004 | type II 3‑dehydroquinate dehydratase/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,398,903 | 0 | C | G | 15.1% | 80.8 / 6.4 | 33 | intergenic (+123/‑119) | aroQ‑III/PP_3004 | type II 3‑dehydroquinate dehydratase/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (19/9); new base G (2/3); total (21/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.28e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01 |
GTGCGGCGATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCG‑TGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTCAA > NC_002947/3398758‑3399048 | gtgcgGCGATGAGCATGCTCTGAGCAGGGCGCAGCGGTCGGATCTCCCAGCGCGTCTACCCGCGAAGCAGGTCGCGCGGGGCATGCGGCGACCGCGCCTGAAGTCACCGTAGGTGTCGTGGATGTCG‑TGCTGCCTACCCATCCAGCatg > 4:173468/1‑149 (MQ=255) gcgATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCG‑TGCTGCCTACCCATCCATGctggatgg < 3:63351/149‑1 (MQ=255) aGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGCTACTGGCTTTCCATCCATGctggatgggtag < 1:150311/148‑5 (MQ=38) gCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCG‑TGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACtt < 5:286867/148‑1 (MQ=255) gCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCG‑TGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACtt > 6:286867/1‑148 (MQ=255) tCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCG‑TGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAg > 6:166900/1‑148 (MQ=255) gAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGAGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCG‑TGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGAATGCATAGTACGGTCCGcc > 2:81943/1‑149 (MQ=255) cgcgAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCG‑TGCTGCCTACCCATCCAGCATGGATGGaaa > 6:70066/1‑97 (MQ=255) cgcgAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCG‑TGCTGCCTACCCATCCAGCATGGATGGaaa < 5:70066/97‑1 (MQ=255) gcgAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCG‑TGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCAcc > 5:169095/1‑149 (MQ=255) gcgcGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCG‑TGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACACGCCACCCGGTTCGTTGcc < 3:173468/149‑1 (MQ=255) gcgcGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCG‑TGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGAATGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGc < 1:81943/148‑1 (MQ=255) tGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCG‑TGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCa < 6:169095/149‑1 (MQ=255) ccGTAGGTGTGGTGGATGTCG‑TGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTg > 3:180086/1‑149 (MQ=255) gTAGGTGTGGTGGATGTCG‑TGCTGCCTACCCATCCAGCACGGATGGAAAGCCAGTAGCGACTTCCGTGAc < 4:268739/70‑1 (MQ=255) tgtgGTGGATGTCG‑TGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCAACCCGGTTCGATGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGACGTc > 4:213138/1‑149 (MQ=255) gATGTCG‑TGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGa > 5:103540/1‑149 (MQ=255) tGTCG‑TGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCGCTCATGACATGGGAGGTCGAATGGACg > 2:238158/1‑149 (MQ=255) gTCG‑TGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTgg > 2:58514/1‑133 (MQ=255) tCG‑TGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTc > 5:97569/1‑149 (MQ=255) tgctgcCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGc > 6:98427/1‑148 (MQ=255) ccTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGt > 6:70241/1‑82 (MQ=255) ccTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACaa < 4:180086/148‑1 (MQ=255) ccTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACaa < 6:97569/148‑1 (MQ=255) ccTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTAAAACCCAACGGACCGGGGCTCGCATGACCTGGGAGGTCGAATGGACGACGCTACAAc > 5:293301/1‑149 (MQ=255) tttccatccatGCTGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTa > 2:85406/4‑143 (MQ=255) tttccatccatGCTGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTa < 1:85406/140‑1 (MQ=255) tttccatccatGCTGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGa > 2:126837/4‑149 (MQ=255) tACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGAc < 2:270959/136‑1 (MQ=255) tACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGAc > 1:270959/1‑136 (MQ=255) aCCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGAt > 4:299044/1‑149 (MQ=255) tCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGt > 5:74660/1‑148 (MQ=255) aGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTcaa > 1:294274/1‑148 (MQ=255) | GTGCGGCGATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCG‑TGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTCAA > NC_002947/3398758‑3399048 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |