| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913_3_hsa_tpiA | 1,544,022 | 1 bp→AA | coding (39/1389 nt) | narU ← | nitrate/nitrite transporter |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913_3_hsa_tpiA | 1,544,021 | 1 | . | A | 86.2% | 84.9 / 9.4 | 29 | coding (40/1389 nt) | narU | nitrate/nitrite transporter |
| Reads supporting (aligned to +/- strand): ref base . (2/2); new base A (14/11); total (16/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01 | |||||||||||
| * | NC_000913_3_hsa_tpiA | 1,544,022 | 0 | C | A | 100.0% | 96.4 / NA | 29 | L13F (TTG→TTT) | narU | nitrate/nitrite transporter |
| Reads supporting (aligned to +/- strand): ref base C (0/0); new base A (16/13); total (16/13) | |||||||||||
TTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATT > NC_000913_3_hsa_tpiA/1543885‑1544166 ||| ttAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAAc < 2:250442/149‑1 (MQ=255) aaCGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGa > 1:520650/1‑149 (MQ=255) gACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCata < 2:87898/149‑1 (MQ=255) aGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCt < 2:520650/149‑1 (MQ=255) ggCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTc > 1:554409/1‑149 (MQ=255) tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTTCAGTGCCATATTGTTCCTCACATGCACACATTGGTa > 2:557767/1‑149 (MQ=255) atatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaa > 1:450598/1‑149 (MQ=255) tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaa < 1:284992/149‑1 (MQ=255) cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCGCATATGCACACATTGGTAATGAAAAAAAg > 1:390128/1‑149 (MQ=255) gagaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACGATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAc > 2:501434/1‑149 (MQ=255) tttCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaaa > 1:174160/1‑149 (MQ=255) tcttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAaca > 2:390649/1‑149 (MQ=255) tcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa < 1:390649/149‑1 (MQ=255) cAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTaa > 1:267187/1‑149 (MQ=255) atatGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTAGCG‑AAAAAAGATAACAACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTTACATGCCCACATTGGTACTGAAAAAAGAACAAAACACgaaggaagt > 2:116952/1‑148 (MQ=255) ttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAGTGAAACAAGGACAAAACACGAGGTAACGGGCAAtg > 1:232131/1‑148 (MQ=255) cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGcc < 2:425616/149‑1 (MQ=255) tttATTTTCGCAGTAGGCAGAATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGTCTATTTTTCTGATTTTGCAGTGCCATACTGTTACTCACATGCACACATAGCTCATGACAAAAAGACAAAACACGAGGTAATACGCAATAGCCa > 2:493133/1‑149 (MQ=255) aTTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAacac < 2:321768/126‑1 (MQ=255) aTTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAacac > 1:321806/1‑126 (MQ=255) aTTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAacac > 1:321768/1‑126 (MQ=255) aTTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGtt > 2:496464/1‑149 (MQ=255) aTTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTAGTAATAAAAAAAAGACAAAacac < 2:321806/126‑1 (MQ=255) ggCCGGATTTTCAGGTTTCCAGTAGCGAA‑AAAAGATAACAACTATTTTTCTCAGTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg < 2:232131/149‑1 (MQ=255) tttCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAAtt < 1:496464/149‑1 (MQ=255) gtCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATATTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACTCGCGGTAAGGCGCAATAGCCAATTATTAGAATTAAGGATGAATTg < 1:493133/148‑1 (MQ=255) tGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGt < 1:387007/149‑1 (MQ=255) cAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAAt < 1:407881/77‑1 (MQ=255) cAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAAt > 2:407881/1‑77 (MQ=255) gcgcgaaaaaaGATAACGACTATTTTTCGCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTACTGAAAAAAAGACAAAACACGAGGTACGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGa > 2:23495/8‑149 (MQ=255) gcgaaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACATAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTACGGATGAATTGGGTGAAGTGCTGAtt > 2:435423/6‑149 (MQ=255) aaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa > 1:579775/3‑77 (MQ=255) aaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa < 2:579775/75‑1 (MQ=255) ||| TTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATT > NC_000913_3_hsa_tpiA/1543885‑1544166 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |