| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 1,073,307 | Δ1 bp | coding (705/1149 nt) | rutA ← | pyrimidine oxygenase, FMN‑dependent |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,073,305 | 0 | A | . | 95.7% | 197.5 / 2.5 | 47 | coding (707/1149 nt) | rutA | pyrimidine oxygenase, FMN‑dependent |
| Reads supporting (aligned to +/- strand): ref base A (0/2); new base . (25/20); total (25/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.14e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01 | |||||||||||
TCGTCGGTTTCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAG > NC_000913/1073206‑1073400 | tcgtcgGTTTCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGa > 2:357573/1‑100 (MQ=255) gtcgGTTTCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGaag > 2:478382/1‑99 (MQ=255) gtcgGTTTCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGaag > 1:1099953/1‑99 (MQ=255) tCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTggg < 2:1430932/100‑1 (MQ=255) cGGCAATCACCATAAACAACCCATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTa > 2:964338/1‑100 (MQ=255) cGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTa > 2:538068/1‑100 (MQ=255) cGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTa < 1:501662/100‑1 (MQ=255) aTCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTAc < 1:43943/100‑1 (MQ=255) tCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACg < 2:1304501/100‑1 (MQ=255) cACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACg < 2:4385/99‑1 (MQ=255) cACCATAAACAACACATAAGAGCCAACGTCGCGTCCGCTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGc > 2:844007/1‑100 (MQ=255) aCCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGcc > 1:1288495/1‑100 (MQ=255) ccATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGcc < 2:483820/100‑1 (MQ=255) aTAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCttt > 2:1174550/1‑100 (MQ=255) cacaTAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAAc < 2:103825/100‑1 (MQ=255) caTAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAg > 2:267527/1‑100 (MQ=255) aTAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGa > 2:254372/1‑100 (MQ=255) tAAGAGCCAATGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGaa > 2:1399004/1‑100 (MQ=255) aaGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAg > 1:146289/1‑100 (MQ=255) aaGAGCCAACGGCGCGTCCGGTTTGCGGTACGGCCTGTTTCATCCGCGCAGCGGTCCGGGCG‑AAGCCGTGAGTGTATTTACGCCTTTGCCGAAACAg < 2:1396082/97‑1 (MQ=255) agagCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGt > 2:1312859/1‑100 (MQ=255) gggCCAACGTCGCGTCCGGTTGGCTCTGCGGCCTGTTGCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTGCGCCTTTGCCGAAACa < 2:1283543/93‑1 (MQ=255) aaCGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAAt < 2:517370/100‑1 (MQ=255) cGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCg < 2:534728/100‑1 (MQ=255) tCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGc < 1:282589/100‑1 (MQ=255) tCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGc < 1:424567/100‑1 (MQ=255) cgcgTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCa > 1:266347/1‑100 (MQ=255) cgTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGTCGAAACAGAAGTTGAAATCGGCATa < 1:537956/100‑1 (MQ=255) gTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATAc > 1:570257/1‑100 (MQ=255) cGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTTAAATCGGCATACCgg > 1:1553939/1‑100 (MQ=255) ggTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCggg > 1:511946/1‑100 (MQ=255) tGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGga > 1:293820/1‑100 (MQ=255) gCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGgag > 2:591474/1‑100 (MQ=255) gCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGgag > 1:524152/1‑100 (MQ=255) gCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGgag < 1:1206993/100‑1 (MQ=255) ctctGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGgaga > 2:37439/1‑100 (MQ=255) tGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACg < 2:21877/100‑1 (MQ=255) cGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGcc < 2:1028837/100‑1 (MQ=255) tttCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAg < 1:80309/100‑1 (MQ=255) tCCGGGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAAGCGGCATACCGGGCGGAGAACGCCATGCCAGCGTc > 2:1174700/1‑99 (MQ=255) ccGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGc > 2:543115/1‑100 (MQ=255) cgcAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGATCGCCATGCCAGCGTCGCttt < 2:609627/100‑1 (MQ=255) aGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGcc > 1:619200/1‑100 (MQ=255) ggTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCgg > 1:1510761/1‑100 (MQ=255) gTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGgc < 1:541183/100‑1 (MQ=255) gTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGgc > 1:1483064/1‑100 (MQ=255) gTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGgc > 1:1153070/1‑100 (MQ=255) ggggCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCa > 2:378569/1‑100 (MQ=255) gggCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAg < 2:619198/100‑1 (MQ=255) gggCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAg < 1:1399012/100‑1 (MQ=255) | TCGTCGGTTTCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAG > NC_000913/1073206‑1073400 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |