Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,718,336 | A→C | *237E (TAA→GAA) | yiaF ← | DUF3053 domain‑containing protein YiaF |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,718,336 | 0 | A | C | 91.7% | 62.6 / ‑2.0 | 24 | *237E (TAA→GAA) | yiaF | DUF3053 domain‑containing protein YiaF |
Reads supporting (aligned to +/- strand): ref base A (1/1); new base C (16/6); total (17/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.07e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01 |
GGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCACCTGAGTACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTGTTGCGTGAAAGTTTGTGCCGCCGGAATCAACGGTTGCAATGCATCAGCTGGCGTCGTT > NC_000913/3718278‑3718574 | ggACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGATGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCACCTGAGTACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTGTTGCGTGAAAGt > 1:346694/1‑251 (MQ=255) aagaaaatgccgatttgcggccgcaTCGCTCGGGCTTTTACTCTTTCTTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCACCTGAGTACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTGTTGCGTGAAAGTTTGTGCCGCCgg > 1:163809/25‑251 (MQ=255) gaaaatgccgatttgaggccgcaTCGCTCGGGCTTTTACTCTTTCTTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTACACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCACCTGAGTACCTTGCTGAGCACTGTAGTCTCCAACCATCACCAGTTGTTGCGTGAAAGTTTGTGCCGCCGGa > 2:99389/23‑250 (MQ=255) aaatgccgatttgcggccgcaTCGCTCGGGCTTTTACTCTTTCTTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCACCTGAGTACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTGTTGCGTGAAAGTTTGTGCCGCCGGAATc > 1:111251/21‑251 (MQ=255) atgccgatttgcggccgcaTCGCTCGGGCTTTTACTCTTTCTTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCACCTGAGTACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTGTTGCGTGAAAGTTTGTGCCGCCGGAATCa > 1:2844/19‑250 (MQ=255) cgatttgcggccgcaTCGCTCGGGCTTTTACTCTTTCTTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCACCTGAGTACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTGTTGCGTAAAAGTTTGTGCCGCCAGAATCAACGGt < 2:405208/237‑1 (MQ=255) gTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCTCCTGAGTTCCCTGCTGAGCAATGTAGTCTCCAACCATCCCCAGTTGTTGCGTGAAAGTTTGTGCCGCCGGAATCAACGGtt < 2:346694/251‑1 (MQ=255) atttgcggccgcaTCGCTCGGGCTTTTACTCTTTCTTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCACCTGAGTACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTGTTGCGTGAAAGTTTGTGCCGCCGGAATCAACGGTTg > 2:403172/13‑251 (MQ=255) ttgcggccgcaTCGCTCGGGCTTTTACTCTTTCTTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCACCGGAGTACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTGTTGCGTGAAAGTTTGTGCCGCCGGAATCAACGGTTGc < 2:265236/240‑1 (MQ=255) cggccgcaTCGCTCGGGCTTTTACTCTTTCTTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCACCTGAGTACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTGTTGCGTGAAAGTTTGTGCCGCCGGAATCAACGGTTGCAAt < 1:200580/243‑1 (MQ=255) cggccgcaTCGCTCGGGCTTTTACTCTTTCTTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCACCTGAGTACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTGTTGCGTGAAAGTTTGTGCCGCCGGAATCAACGGTTGCAAt > 2:200580/8‑250 (MQ=255) ggccgcaTCGCTCGGGCTTTTACTCTTTCTTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAg < 1:422765/124‑1 (MQ=255) ggccgcaTCGCTCGGGCTTTTACTCTTTCTTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAg > 2:422765/7‑130 (MQ=255) ggccgcaTCGCTCGGGCTTTTACTCTTTCTTGGGTTGCAGTAAATGCTGCAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCACTCAGATTGTTATATTCACTTGCCTGATGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTAACCTGAGTACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTGTTGCGTGAAAGCTTGTGCCGCCGGAATCAACGGTTCCAATGc > 1:557360/7‑251 (MQ=255) ggccgcaTCGCTCGGGCTTTTACTCTTTCTTGGGTTGCAGCAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCACTCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCACCTGAGTACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTGTTGCGTGAAAGTTTGTGCCGCCGGAATCAACGGTTGCAATGc > 1:524353/7‑251 (MQ=255) gccgcaTCGCTCGGGCTTTTACTCTTTCTTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCACCTGAGTACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTGTTGCGTGAAAGTTTGTGCCGCCGGAATCAACGGTTGCAATGc > 1:425176/6‑250 (MQ=255) cgcaTCGCTCGGGCTTTTACTCTTTCTTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCACCTGAGTACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTGTTGCGTGAAAGTTTGTGCCGCCGGAATCAACGGTTGCAATGCATc > 1:473838/4‑251 (MQ=255) gCTCGGGCTTTTACTCTTTCTTGCGTTGCAGTAACTGCTGTAGTACCGGCCTGATTAAACGCCTGATGCTGTGCCGGTCATGGCGCAATCAGTTTGTTATATGCACTTGCCTGCTGGGCAGTCGGGAACTGGATGCCATTTGCGACAATGCTCACATGAGTAACTTGCTGAGCCATGTAGTCTACAGCCATCACCAGTTGATGCTTGTAAGTTTGTGTCGCCGGCAGCAACGGTTGAAATTCAGCAGCttg > 2:218658/1‑249 (MQ=255) cTCGGGCTTTTACTCTTTCTTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCACCTGAGTACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTGTTGCGTGAAAGTTTGTGCCGCCGGAATCAACGGTTGCAATGCATCAGCTGGc > 2:460008/1‑251 (MQ=255) gggCTTTTACTTTTTCTTGGGTTGCAGTCACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGACCTGGATGCCATTTGCGCCAAAGCTCACCTGAGTACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTATTGCGTGAAAGTTTGTGCCGCCGGAATCAACGGTTGCAATGCATCTGCTGGcgtc > 2:320819/1‑251 (MQ=255) gggCTTTTACTCTTTCTTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCACCTGAGTACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTGTTGCGTGAAAGGTTGTGCCGCCGGAATCAACGGTTGCAATGCATCAGCTGGcgtc < 1:402856/251‑1 (MQ=255) gggCTTTTACTCTTTCTTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCACCTGAGAACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTGTTGCGTGAAAGTTTGTGCCGCCGGAATCAACGGTTGCAATGCATCAGCTGGcgtc < 1:421397/251‑1 (MQ=255) gCTTTTACTCTTTCTTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCACCTGAGTACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTGTTGCGTGACAGTTTGTGCCGCCGGAATCAACGGTTGCAATGCATCAGCTGGcgtcgt > 2:55478/1‑251 (MQ=255) cTTTTACTCTTTCTTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGCTTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCACCTGAGTACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTGTTGCGTGAAAGTTTGTGCCGCCGGAATCAACGGTTGCAATGCATCAGCTGGCGTCGtt > 2:425633/1‑251 (MQ=255) | GGACTAAGCCGCGACTGCGTAGAGTAAAGCCCGATAATCGCTCGGGCTTTTACTCTTTATTGGGTTGCAGTAACTGCTGTAGTCCAGGCCTGATTAAACGCCTGATGCTGTGCCGGTAATGGCGCAATCAGTTTGTTATATTCACTTGCCTGCTGTGAAGTCGGGAACTGGATGCCATTTGCGACAAAGCTCACCTGAGTACCTTGCTGAGCAATGTAGTCTCCAACCATCACCAGTTGTTGCGTGAAAGTTTGTGCCGCCGGAATCAACGGTTGCAATGCATCAGCTGGCGTCGTT > NC_000913/3718278‑3718574 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |