Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 772,568 | C→A | Q371K (CAG→AAG) | cydA → | cytochrome bd‑I ubiquinol oxidase subunit I |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 772,568 | 0 | C | A | 97.0% | 111.3 / ‑2.8 | 33 | Q371K (CAG→AAG) | cydA | cytochrome bd‑I ubiquinol oxidase subunit I |
Reads supporting (aligned to +/- strand): ref base C (0/1); new base A (15/17); total (15/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.03e-01 |
TCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCAT > NC_000913/772435‑772704 | tctGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTc > 1:33276/1‑140 (MQ=255) tGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCaa > 1:242750/1‑140 (MQ=255) aacaacTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGGCTGAAGCGAAGATTCAACAGGCa > 1:615135/1‑140 (MQ=255) acaacTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCaa < 1:215734/140‑1 (MQ=255) acaacTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCaa < 2:587856/140‑1 (MQ=255) tGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCaaa < 2:242750/140‑1 (MQ=255) tCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTc > 1:110512/1‑140 (MQ=255) cGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCgtg < 1:149512/140‑1 (MQ=255) gACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCgtgt > 2:511417/1‑140 (MQ=255) gACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCgtgt > 2:370813/1‑140 (MQ=255) ccAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAg < 2:552013/140‑1 (MQ=255) ccAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAg < 1:238843/140‑1 (MQ=255) ccAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAg < 1:511417/140‑1 (MQ=255) cAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAgc > 2:145258/1‑140 (MQ=255) aGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAgcg < 2:364518/140‑1 (MQ=255) aCCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACttt < 2:110512/140‑1 (MQ=255) gTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGt > 1:21631/1‑140 (MQ=255) agaaagaCCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTa < 1:439155/132‑1 (MQ=255) agaaagaCCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTa > 2:439155/1‑132 (MQ=255) aCCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCgtgt > 1:655004/1‑140 (MQ=255) aCCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCgtgt < 1:219677/140‑1 (MQ=255) gTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCgt > 1:372336/1‑132 (MQ=255) gTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCgt < 2:372336/132‑1 (MQ=255) gTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTc > 1:415134/1‑140 (MQ=255) tACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCt < 2:415134/140‑1 (MQ=255) cGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGc > 1:236976/1‑140 (MQ=255) ctgctgAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGc > 1:461753/1‑140 (MQ=255) tgAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAatc < 2:461753/140‑1 (MQ=255) gCTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATcgcg < 2:414156/140‑1 (MQ=255) cTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCgcgc < 2:33276/140‑1 (MQ=255) cTATACGCCAAACGTGGCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCgcgc < 2:21631/140‑1 (MQ=255) ggCTGATGCGACTGAAGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAgtg > 2:68031/1‑140 (MQ=255) aGCGAAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCa < 2:317525/140‑1 (MQ=255) gcgaAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCAt < 1:318676/140‑1 (MQ=255) | TCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGCGACTGAAGCGCAGATTCAACAGGCAACCAAAGACTCCATCCCGCGTGTAGCGCCGCTGTACTTTGCGTTCCGTATCATGGTGGCGTGTGGCTTCCTGCTTCTGGCAATCATCGCGCTCTCTTTCTGGAGTGTCATCCGCAACCGCAT > NC_000913/772435‑772704 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |