| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,129,710 | 0 | T | A | 21.4% | 55.3 / 11.9 | 28 | pseudogene (434/1101 nt) | yghO | pseudogene, IS‑interrupted |
| Reads supporting (aligned to +/- strand): ref base T (12/10); new base A (3/3); total (15/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCTGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTGATCTTACTTGCACCTTGTGACTTCAACCACG > NC_000913/3129573‑3129842 | aGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCTGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTg > 1:102801/1‑149 (MQ=255) ggTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGACAGTGAGAGCGGTGCGCTGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGGTCAATATGCGCGGCATACCACGGTTTGCCgtgt < 5:110941/148‑1 (MQ=255) ttGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCTGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCg < 2:102801/148‑1 (MQ=255) gcACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCTGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCAAACCGTGGTATGCCgc < 1:30607/132‑2 (MQ=255) gcACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCTGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCAAACCGTGGTATGCCgc > 2:30607/1‑131 (MQ=255) cAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCTGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCa > 1:73398/1‑148 (MQ=255) gTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCTGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCCTCATCGCACAGggtggt > 4:45988/1‑140 (MQ=255) gTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCTGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGggtggt < 3:45988/140‑1 (MQ=255) cGGAGAGAAAGTGAGATCGGTTCGCTGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAAt > 6:71108/1‑144 (MQ=255) cGGAGAGAAAGTGAGATCGGTTCGCTGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAAt < 5:71108/144‑1 (MQ=255) tGAGATCGGTTCGCTGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTcc < 2:73398/149‑1 (MQ=255) cGGTTCGCTGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCAAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCtgatt > 5:96552/1‑148 (MQ=255) cTGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGgtgt < 1:123205/105‑1 (MQ=255) cTGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGgtgt > 2:123205/1‑105 (MQ=255) accacGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGa > 3:59934/1‑147 (MQ=255) cacGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAg > 7:5497/1‑148 (MQ=255) acGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGa < 6:96552/149‑1 (MQ=255) cAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGc < 8:5497/149‑1 (MQ=255) cccAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAgg < 3:98134/111‑1 (MQ=255) cccAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAgg > 4:98134/1‑111 (MQ=255) cccAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTGATCTTACTTGCACCTTGTGACt > 5:88527/1‑148 (MQ=255) ccAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTGATCTTACtt > 5:14690/1‑134 (MQ=255) ccAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTGATCTTACtt < 6:14690/134‑1 (MQ=255) ttCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTGATCTTACTTGCACCTTGTGACTTCAACCACg > 1:42114/1‑149 (MQ=255) acgGCAAACCGTGGTATGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTGATCTTACTTGc > 5:140471/2‑119 (MQ=255) acgGCAAACCGTGGTATGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTGATCTTACTTGc < 6:140471/118‑1 (MQ=255) acgGCAAACCGTGGTATGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTGATCTTACTTGCAcc < 4:139123/121‑1 (MQ=255) acgGCAAACCGTGGTATGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTGATCTTACTTGCAcc > 3:139123/2‑122 (MQ=255) | AGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCTGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTGATCTTACTTGCACCTTGTGACTTCAACCACG > NC_000913/3129573‑3129842 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |