| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913_3_vch_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_vch_tpiA | 1,544,021 | 1 | . | A | 91.7% | 70.6 / 0.3 | 24 | coding (40/1389 nt) | narU | nitrate/nitrite transporter |
| Reads supporting (aligned to +/- strand): ref base . (1/1); new base A (11/11); total (12/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01 | |||||||||||
| * | NC_000913_3_vch_tpiA | 1,544,022 | 0 | C | A | 100.0% | 74.9 / NA | 24 | L13F (TTG→TTT) | narU | nitrate/nitrite transporter |
| Reads supporting (aligned to +/- strand): ref base C (0/0); new base A (12/12); total (12/12) | |||||||||||
GTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA > NC_000913_3_vch_tpiA/1543891‑1544170 |||| gTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTAt > 1:537076/1‑149 (MQ=255) aCTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTAtttt > 2:91233/1‑149 (MQ=255) tCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCGTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGc > 1:112126/1‑149 (MQ=255) cgaGACAGAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGcc < 2:133884/147‑1 (MQ=255) aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTATTCCTcaca > 1:293322/1‑149 (MQ=255) gTAGACAACTGACCGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcac < 1:450557/149‑1 (MQ=255) caGAAAACTGACTGATATCAAGAGAGTTCTTCGAGCAATATGATTTACTTTATTGTACCAGAAGGGCGGTTTTTC‑GGTTTAAAAAAAAA‑AAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcaca < 2:107367/148‑1 (MQ=255) gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacaca > 2:210066/1‑149 (MQ=255) aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg < 1:498493/149‑1 (MQ=255) ctgactgaTATCCAGCGATTTCTTCGAGCAATCTGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGt < 1:75112/149‑1 (MQ=255) tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGa < 1:91233/149‑1 (MQ=255) ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa > 1:302669/1‑149 (MQ=255) acAGAGATTTCTTCGAGCAATATGTGTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa < 2:202306/148‑1 (MQ=255) tttCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaaa > 2:307632/1‑149 (MQ=255) tttCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaaa < 2:302669/149‑1 (MQ=255) ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACg > 1:530380/1‑149 (MQ=255) cGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTATTTTTCTCATTTTg > 1:145986/1‑90 (MQ=255) cGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTATTTTTCTCATTTTg < 2:145986/90‑1 (MQ=255) gAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAgg > 1:400366/1‑149 (MQ=255) ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGtt < 1:382360/94‑1 (MQ=255) ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGtt > 2:382360/1‑94 (MQ=255) tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCa < 1:210066/149‑1 (MQ=255) cGGATTTTCTGGTTTCCAGTCGCGAA‑‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAACACAAAAAAAAACAC‑AGATCAAGGCGCAATAGCCCGTTATTAGAATTAAAGag > 2:263206/1‑148 (MQ=255) ctCCAGTCGCGAA‑‑AAAAGATAACGACTCTTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAg < 1:291464/148‑1 (MQ=255) tcgcgaaaaaaGATAACGACTATTTTTCTTATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAAAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGa < 2:279815/142‑1 (MQ=255) tcgcgaaaaaaGAAAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGa < 1:263734/142‑1 (MQ=255) gaaaaaaGATAACGACTATTTTTTTCACTCGAAACGACCATATTGTTCCTCACATGCACGCAATGGTAATGAAAAAACGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTTTTATCCTTACGGGTGAATTGGTGGAGGTGCTGAttca > 2:123893/4‑147 (MQ=255) aaaaaGATAACCACTATTTTTCTCATGTTGCAATGCCATATTGTTCCTCACAAGCACACATTTGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa < 2:503824/148‑1 (MQ=255) |||| GTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGT‑AAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA > NC_000913_3_vch_tpiA/1543891‑1544170 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |