Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913_3_pae_tpiA | 4,037,227 | 1 | . | T | 25.0% | 7.7 / 14.1 | 16 | intergenic (+10/‑33) | ileT/alaT | tRNA‑Ile/tRNA‑Ala |
Reads supporting (aligned to +/- strand): ref base . (6/6); new base T (0/4); total (6/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.34e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.23e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACG‑GCAAATTTGAAGAG > NC_000913_3_pae_tpiA/4037079‑4037241 | agaagCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACg < 1:143943/149‑1 (MQ=34) cGTTCTTTGAAGTGCTCACACAGATTGTCGGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACG‑GCaaa < 2:130672/149‑1 (MQ=21) tCTTTGAAGTGCTCACACAGATTGTCTTATGAAAATGAGCAGTAAAGCCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACG‑GCAAAttt < 1:243187/149‑1 (MQ=17) cTTTGAGGTGCTAACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCCTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGC‑CGTGCAAATTTg < 1:990274/149‑1 (MQ=11) cTTTGAAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGC‑CGTGCAAATTTg < 1:617716/149‑1 (MQ=12) cTTTGAAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGC‑CGTGCAAATTTg < 2:14662/149‑1 (MQ=255) cTTTGAAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGATTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAGGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCCCG‑GCAAATTTg > 2:291221/1‑149 (MQ=12) cTTTGAAGTGCTCACACAGATTGCCTGATTAAAATCAGCAGTAAAACCTTGACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACG‑GCAAATTTg < 2:372901/149‑1 (MQ=11) tttGAAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACG‑GCAAATTTGa > 1:621929/1‑149 (MQ=32) tttGAAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCAACTCAGGCCTACCAAATTTGCACG‑GCAAATTTGa > 1:907299/1‑149 (MQ=17) tttGAAGTGCTCACACAGATAGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGCGGTAGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACG‑GCAAATATGa > 2:995452/1‑149 (MQ=255) tttaaGTGATCACACAGTTTGTCTGATGAAAATGAACAGTAAACCCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCCAGTCCACTCAGGCCTACCAAATTTGCACG‑GCAACATTGaa < 2:13998/146‑1 (MQ=11) ttGAAGTGCTCACACAGATTGTCTGATGAAAATTAGCAGAAAAACCTCTACAGGCTTGTAGCTCCGGAGGTTAGAGCGCACCCCTGATACGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACG‑GCCGATTGGaa > 2:835097/1‑149 (MQ=11) ttGAAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACG‑GCAAATTTGaa < 1:18523/149‑1 (MQ=32) tGAAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGC‑CGTGCAAATTTggta < 2:942052/149‑4 (MQ=255) gAAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACCGGCTTGTAGCTCAGGTGATTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCCCG‑GCACAttgcaggt > 2:726747/1‑143 (MQ=11) aaGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTTAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACG‑GCAAATTTGAagag < 1:373728/149‑1 (MQ=11) | AGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAATGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACG‑GCAAATTTGAAGAG > NC_000913_3_pae_tpiA/4037079‑4037241 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |