Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,401,844 | A→G | I471T (ATC→ACC) | abgB ← | p‑aminobenzoyl‑glutamate hydrolase, B subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,401,844 | 0 | A | G | 94.1% | 161.0 / 0.9 | 51 | I471T (ATC→ACC) | abgB | p‑aminobenzoyl‑glutamate hydrolase, B subunit |
Reads supporting (aligned to +/- strand): ref base A (3/0); new base G (25/23); total (28/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.42e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01 |
ACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGG > NC_000913/1401752‑1401937 | acgGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCa > 2:754811/1‑100 (MQ=255) gTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTg < 2:738802/100‑1 (MQ=255) gTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTg > 2:509694/1‑100 (MQ=255) ggATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATAc > 1:902961/1‑100 (MQ=255) aCTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCgtg < 1:642636/100‑1 (MQ=255) aCTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCgag > 1:1110524/1‑98 (MQ=255) aCTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCg < 2:1110514/98‑1 (MQ=255) cATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTcc > 1:833977/1‑100 (MQ=255) cATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTc > 1:994792/1‑99 (MQ=255) aTGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCg < 2:1211484/100‑1 (MQ=255) tGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTcc < 2:778915/98‑1 (MQ=255) gggCATTCCTCGGTTGTTGTGTTTGCGTTGTTATTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGtt > 1:1235673/1‑100 (MQ=255) ggCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTGTTCGGAGTAGGGCAGTGATACGGTTGCGTGTCCGTTa > 2:876522/1‑100 (MQ=255) cATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTAGCGTGTCCGTTACt > 2:887192/1‑100 (MQ=255) ttCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTg < 2:571944/100‑1 (MQ=255) ttCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTg > 1:420651/1‑100 (MQ=255) ttCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTg > 1:1231549/1‑100 (MQ=255) ttCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTg < 1:865968/100‑1 (MQ=255) ccTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCt < 2:100062/100‑1 (MQ=255) ccTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCt > 2:84811/1‑100 (MQ=255) cTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTg < 2:420647/100‑1 (MQ=255) cTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTg < 1:104495/100‑1 (MQ=255) cGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGAt < 1:731980/100‑1 (MQ=255) cGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGAt < 2:283721/100‑1 (MQ=255) ggTTGTTGTGTTTGCGATGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATg > 1:25943/1‑100 (MQ=255) tgtgTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCttg < 2:902954/100‑1 (MQ=255) tgtTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCttgtt > 2:456765/1‑100 (MQ=255) tgtTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCttgt > 2:896374/1‑99 (MQ=255) gtTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCttgttg < 2:701826/100‑1 (MQ=255) ttgttgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCt < 1:76540/100‑1 (MQ=255) gttgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTg < 1:962696/100‑1 (MQ=255) gttgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTg > 2:638958/1‑100 (MQ=255) ttgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGc < 2:1231538/100‑1 (MQ=255) tgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCa < 1:458812/100‑1 (MQ=255) tgttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCa > 1:63284/1‑100 (MQ=255) tgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATa > 2:1295364/1‑100 (MQ=255) gttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAg > 1:620779/1‑100 (MQ=255) gttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAg > 2:534701/1‑100 (MQ=255) ttttAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCt > 2:937718/1‑100 (MQ=255) ttttAAAGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCt > 1:798641/1‑100 (MQ=255) tttAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTg > 2:543860/1‑100 (MQ=255) ttAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGa > 1:139016/1‑100 (MQ=255) aaaGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAAt > 2:957672/1‑100 (MQ=255) aGGTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCa > 1:1040997/1‑100 (MQ=255) gTGACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAg > 1:59411/1‑100 (MQ=255) gACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCa > 2:536622/1‑96 (MQ=255) gACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGa < 2:708836/100‑1 (MQ=255) gACGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGa < 1:1102199/100‑1 (MQ=255) cGGTGTCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAag < 1:1144485/100‑1 (MQ=255) gtgtCACGTTTTTCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAagag < 1:299555/100‑1 (MQ=255) ttCGGGGTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTgg < 1:896381/100‑1 (MQ=255) | ACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGG > NC_000913/1401752‑1401937 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |