| 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 | 76.9% | 66.0 / 15.4 | 26 | coding (40/1389 nt) | narU | nitrate/nitrite transporter |
| Reads supporting (aligned to +/- strand): ref base . (6/0); new base A (9/11); total (15/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.37e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| * | NC_000913_3_vch_tpiA | 1,544,022 | 0 | C | A | 100.0% | 80.5 / NA | 26 | L13F (TTG→TTT) | narU | nitrate/nitrite transporter |
| Reads supporting (aligned to +/- strand): ref base C (0/0); new base A (15/11); total (15/11) | |||||||||||
TTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑C‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGA > NC_000913_3_vch_tpiA/1543879‑1544168 ||| ttGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑A‑a‑aaaa > 1:590890/1‑148 (MQ=255)ttGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCATTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑A‑a‑aaaa > 1:480961/1‑148 (MQ=255) cGGTAACTGCGCTAAATAGCATCCAGCCACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AAAGATAACGACt < 1:589888/149‑1 (MQ=255) tAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTCCGAGCAATATGTTTTCCTTTATTTTCCCGGAAAGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AAAGATAACGACTAtt > 1:125785/1‑149 (MQ=255) tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCCGTCG‑CG‑AAA‑AAAGATAACGCCTATTTTtc > 2:130205/1‑149 (MQ=255) gcgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AAAGATAACGACTATTTTtct < 1:31327/149‑1 (MQ=255) tCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGc > 2:15554/1‑149 (MQ=255) aaacaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCatat < 1:207637/147‑1 (MQ=255) acaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCCA‑A‑A‑AAAGATAACGACTCTTTTTCTCCTTTTGAACGGCCatatct > 2:178841/1‑147 (MQ=255) caGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AAAGATCGCTACTGTTTTTCTCATTTTGCAGTGCCATATTGt < 2:317407/149‑1 (MQ=255) gAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc < 2:239312/149‑1 (MQ=255) aaGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑AAA‑AAAGATAACGCCTAATCTTATCATTTTGCAGAGCCATATTGTTcc > 1:293806/1‑149 (MQ=255) gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca < 2:524944/149‑1 (MQ=255) aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca > 1:513758/1‑149 (MQ=255) aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca < 1:70251/149‑1 (MQ=255) tAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCCTATTGTTCCTCACATGcaca > 1:318148/1‑149 (MQ=255) aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGTCGGATTTTCTGGTTTCCCGTCG‑CG‑A‑AGAAAGATAACGACTATTTTTCTCATTTTGCAATGCTCTATCGCTCGACACTAGCACACCTTg > 1:345088/1‑149 (MQ=255) tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCATTGCCATATTGTTCCTCACATGCACACATTGGTa > 2:575889/1‑149 (MQ=255) tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc > 2:659157/1‑125 (MQ=255) tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc < 1:659157/125‑1 (MQ=255) atCCAGAGCTTTCTTCGAGCAATCTGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AAAGATAAAGACTGTTTTTCTCATTTTTCAGTGCCATATTGTTCATCCCAAGAGAACATTGGTAATGtaaca > 2:309820/1‑147 (MQ=255) tcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa > 2:129959/1‑149 (MQ=255) tcGAGCAATATGTTTTCCTTTATTTTACAAAAAGAAAGAATTTTCTGGTTTCCAGTGG‑CGAA‑A‑AAAGATAACGACTATTTTCCTCATTTTGCAGTGCCTTATTGTTCCCCTCTTGCACTGAAAAAAAAAAAAAAAAAGACAAAACACGa < 2:318148/149‑1 (MQ=255) aGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGt > 1:172709/1‑149 (MQ=255) ttCTGGTTTCCAGTCG‑CGAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTg < 2:222835/149‑1 (MQ=255) ttCCAGCCG‑CGAA‑A‑AAAGATAACGACTATTTTGATCATGTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAg < 2:542982/149‑1 (MQ=255) tcgcgaaa‑aaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGa > 1:546375/8‑149 (MQ=255) gcgaaa‑aaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAtt > 2:556144/6‑149 (MQ=255) gaaa‑aaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGa < 1:129959/146‑1 (MQ=255) ||| TTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑C‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGA > NC_000913_3_vch_tpiA/1543879‑1544168 |
| 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: ‑ |