Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913_3_pae_tpiA | 275,065 | G→A | T18I (ACC→ATC) | insH1 ← | IS5 transposase and trans‑activator |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913_3_pae_tpiA | 275,065 | 0 | G | A | 94.7% | 46.6 / ‑6.8 | 19 | T18I (ACC→ATC) | insH1 | IS5 transposase and trans‑activator |
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (7/11); total (8/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.21e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTATTA > NC_000913_3_pae_tpiA/275002‑275209 | ggggAGAAGAAAATCTGTTTTTGGGTCTGACGGCGCTGACTGGTGAAGTCACTGGCGGCGAAGATAAGTAGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAAGATGATCTCATACCAGGGACTTGGTCGCACCTTCCg < 2:43463/148‑1 (MQ=11) cGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCg < 2:836627/149‑1 (MQ=11) gggACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCgg > 2:770260/1‑149 (MQ=14) gggACAAGAAAATCGCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCgg < 2:951982/149‑1 (MQ=14) ggACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGa < 2:804255/149‑1 (MQ=17) gAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGtt > 1:396448/1‑149 (MQ=37) aaaaTCGCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTa < 2:623768/149‑1 (MQ=255) aTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGc > 1:154844/1‑142 (MQ=21) aTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGc < 2:154844/142‑1 (MQ=21) aCGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTg > 1:855030/1‑149 (MQ=255) ggCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTa < 1:15706/149‑1 (MQ=255) gcgcTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTAt < 1:892589/149‑1 (MQ=255) aCTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGatt < 1:826665/149‑1 (MQ=255) aCTGCTGAATTCACCGTGGGCGAAGATAAGTTGAAGACTCCTGATGAACCCTGTTCTATGGCTCCAGATGACAAAGATGATCGCAGATCCGGGACTTGTTCGCATCTTCCGGGGGCGTTATGAGCTGGCGGCCCTTTTTGTGTCTGatt < 2:517354/149‑1 (MQ=17) tgctgAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGAGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTTattct > 2:625502/1‑147 (MQ=255) aaTTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGTCTCGATCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATc > 1:549591/1‑149 (MQ=255) tttaCTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATccc < 1:306800/146‑1 (MQ=255) ttgaCTGTCGGCGAAGGTAAGTTGTTGACTCATGATGAACTCTGTTCTATGGAGACAGATGATAAACATGAACTCATATCAGGTACATGTTAGCACCTTCCGGAAGCGTTATGTAGTGGAGGACCGTTTTGTCTCTGATTATTAATcct > 2:278899/4‑148 (MQ=12) gAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATAATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTATTa > 2:266986/1‑149 (MQ=255) | CGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTATTA > NC_000913_3_pae_tpiA/275002‑275209 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |