Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913_3_pae_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_pae_tpiA | 1,544,021 | 1 | . | A | 85.7% | 74.2 / 8.6 | 28 | coding (40/1389 nt) | narU | nitrate/nitrite transporter |
Reads supporting (aligned to +/- strand): ref base . (4/0); new base A (15/9); total (19/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.73e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01 | |||||||||||
* | NC_000913_3_pae_tpiA | 1,544,022 | 0 | C | A | 100.0% | 84.0 / NA | 28 | L13F (TTG→TTT) | narU | nitrate/nitrite transporter |
Reads supporting (aligned to +/- strand): ref base C (0/0); new base A (19/9); total (19/9) |
TATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGA > NC_000913_3_pae_tpiA/1543877‑1544164 ||| tATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑a‑aaa > 1:570523/1‑149 (MQ=255) aTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATAGCCAGAGATTTCTTCGAGCAATATGATTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATaa > 2:202167/1‑149 (MQ=255) cGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACt > 2:372005/1‑149 (MQ=255) cgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTtctc > 1:165401/1‑149 (MQ=255) ccAGACAGAGAAGGCAAGAAGTAGACAACTGAAAAATATCCTAAGTTCGCGTCGATTAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTGCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGcc < 2:60198/149‑1 (MQ=255) gACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCata < 2:571529/149‑1 (MQ=255) gACACAGAAGGCAAGAAGTAGACAAATGACTTATATCCAGAGATTTCTTTGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTCTCCAGTGCCata < 2:213824/149‑1 (MQ=255) aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCAGTTTGCAGTGCCATATTGtt > 2:460799/1‑149 (MQ=255) aaGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATACCTACTATTTTTCTCATTTTGCAGTGCCATATTGTTct > 2:115588/1‑148 (MQ=255) ggCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTc > 1:167265/1‑149 (MQ=255) gCACGATGTGGACATCTGACTGAAATCCTGAGATTTCGTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCATTCTCGAA‑AAAAGATAACGACTATTTTGCTCATTTTGCAGTGCCATATTGTTCCTca < 2:405355/149‑1 (MQ=255) aagTAGACAACTGATGGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGc < 1:505189/149‑1 (MQ=255) aagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCAGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGc > 1:224218/1‑149 (MQ=255) agTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGca < 1:380918/149‑1 (MQ=255) tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCAAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTa > 1:471940/1‑149 (MQ=255) ctgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACCTGCACACATTGGTAATg > 2:416794/1‑149 (MQ=255) atatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACCCATTGGTAATGaaa > 1:133274/1‑149 (MQ=255) ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg > 1:565352/1‑137 (MQ=255) ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg < 2:565352/137‑1 (MQ=255) cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAg > 2:512023/1‑149 (MQ=255) tttCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaaa > 1:55436/1‑149 (MQ=255) ttcGAGCAATATGTTTTCCTTTATTGTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATAGGTAATGAAAAAAAGACAGAACACg > 1:520094/1‑149 (MQ=255) cAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTaa > 1:171004/1‑149 (MQ=255) ttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTAGCCCGAGCACATGGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatt > 1:277370/1‑149 (MQ=255) cAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGaa > 2:395372/1‑149 (MQ=255) aGGCCGGATTTTCTGGTTCCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTTTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa < 1:115588/149‑1 (MQ=255) ggCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg > 1:465346/1‑149 (MQ=255) ggATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATAGTGTTCCTCACATGCACACATTGGTAATGAAAAAAGGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATg > 1:398959/1‑149 (MQ=255) tttCCAGTCGCGAA‑AAAAAATAAAGAATATTTTTCTCATTTTGAAGTGCCATATTGTCCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGaa < 2:570523/149‑1 (MQ=255) tcgcgaaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGa > 2:322072/8‑149 (MQ=255) ||| TATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGA > NC_000913_3_pae_tpiA/1543877‑1544164 |
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: ‑ |