| 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 | 72.4% | 74.1 / 26.3 | 29 | coding (40/1389 nt) | narU | nitrate/nitrite transporter |
| Reads supporting (aligned to +/- strand): ref base . (5/3); new base A (11/10); total (16/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.97e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAT > NC_000913_3_hsa_tpiA/1543872‑1544165 | gATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTcgcg > 3:161614/1‑150 (MQ=255) aCGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACt < 1:111321/150‑1 (MQ=255) aCTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTAttttt < 1:332695/150‑1 (MQ=255) tAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCAttt < 4:90874/150‑1 (MQ=255) cacaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTg < 3:164019/150‑1 (MQ=255) caGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGtt < 1:142200/150‑1 (MQ=255) aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc < 2:361425/150‑1 (MQ=255) agaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATg < 1:431659/150‑1 (MQ=255) aGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacaca > 1:289776/1‑150 (MQ=255) gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt > 2:321448/1‑150 (MQ=255) cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa > 1:102602/1‑150 (MQ=255) agaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaa > 1:176808/1‑150 (MQ=255) tttCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAaca > 2:10217/1‑150 (MQ=255) tcttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAacac > 1:208843/1‑150 (MQ=255) ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa > 3:242158/1‑150 (MQ=255) ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGCTTGCCAGTCGCG‑C‑AAA‑AGATAACGACTATTTTTCTCAATTTGCAGCGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAg > 3:223449/1‑150 (MQ=255) tcGAGCAATATGTTTTCCTTTATTTTCCCGGAAGGCCGGTTTCTATGGTTAGCAATCACA‑A‑AAAGCGATAACGACTATTCTTCTGATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAg < 1:376414/150‑1 (MQ=255) cAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAg < 1:143898/150‑1 (MQ=255) ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAAt > 1:289455/1‑150 (MQ=255) tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg > 4:155211/1‑150 (MQ=255) ttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGtt < 4:306927/150‑1 (MQ=255) aTTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGtta > 2:433872/1‑150 (MQ=255) ttttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttat < 4:99782/150‑1 (MQ=255) tCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAg > 4:219991/1‑150 (MQ=255) ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGaa < 3:155211/150‑1 (MQ=255) ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGaa > 1:403131/1‑150 (MQ=255) ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAATACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAAt > 2:215969/1‑150 (MQ=255) gAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTa > 2:466668/1‑150 (MQ=255) cGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATg > 4:288275/1‑150 (MQ=255) gTCGCG‑C‑AAA‑AGATACCGACTATTTCTCTTATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTAATAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAt < 4:279583/150‑1 (MQ=255) | GATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAT > NC_000913_3_hsa_tpiA/1543872‑1544165 |
| 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: ‑ |