| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913_3_pae_tpiA | 3,959,515 | 0 | T | A | 20.9% | 88.5 / 25.8 | 43 | intergenic (+70/+17) | ilvC/ppiC | ketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C) |
| Reads supporting (aligned to +/- strand): ref base T (19/15); new base A (4/5); total (23/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.11e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.16e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACC > NC_000913_3_pae_tpiA/3959366‑3959663 | gCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGGCCCTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCt < 3:201056/150‑1 (MQ=255) cGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTgg < 4:102918/143‑1 (MQ=255) cGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTgg > 3:102918/1‑143 (MQ=255) aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCt < 2:203582/150‑1 (MQ=255) aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCt > 4:278163/1‑150 (MQ=255) ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTa < 3:278163/150‑1 (MQ=255) ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCAGGAGAAGGCCTTgagtgtt < 1:363250/150‑7 (MQ=255) gtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGA‑AAACACCCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAg > 4:363586/1‑150 (MQ=255) aggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTgg < 1:653041/130‑1 (MQ=255) aggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTgg > 2:653041/1‑130 (MQ=255) taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCg < 1:213726/150‑1 (MQ=255) cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGA‑AAACACGCAAGGCCTTC‑TCCAGGAGAAGGCCTTgagt > 1:3669/1‑115 (MQ=255) cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGA‑AAACACGCAAGGCCTTC‑TCCAGGAGAAGGCCTTgagt < 2:3669/118‑4 (MQ=255) tATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgt < 2:158096/150‑1 (MQ=255) ttGCTGTTGCGGGTTAAGTGCGCGCTGATATCCTAGCCACGACCCTCTCCAACTGGGAG‑GGGATATTAGTAGTCAGCGCTTTCGTCGTTG‑GGAGCCCTTTCGTTTAGTTTCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgtgt > 3:211155/1‑150 (MQ=255) gTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACGGGGCGCGGGAGA‑AAACACTCAACGCCTTC‑TCCTGGAGAATGCCTTGCTATTAGTGGCGGTAGAGCACCTTACTGATGTGATATCCGAAATGGGTGTGCAGCGGGCCGGTCGGc > 4:436549/1‑150 (MQ=38) gcgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTc > 3:239841/1‑150 (MQ=255) cgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGTAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTcc < 2:109323/150‑1 (MQ=255) aTGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCg < 1:28295/150‑1 (MQ=255) cccTCACCCCGACCCTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCg < 1:148707/110‑1 (MQ=255) cccTCACCCCGACCCTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCg > 2:148707/1‑110 (MQ=255) cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACaa > 1:602601/1‑150 (MQ=255) ttctcccTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaaa > 1:692794/5‑150 (MQ=255) aCCCTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAcc < 4:239841/150‑1 (MQ=255) cccTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCa > 4:456458/1‑150 (MQ=255) tcCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTa > 1:292614/1‑150 (MQ=255) caGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGaa < 1:145429/150‑1 (MQ=255) ggagagggaga‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGc > 1:733959/1‑150 (MQ=255) gaga‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTa > 2:321171/1‑135 (MQ=255) gaga‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTa < 1:321171/135‑1 (MQ=255) tCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCt > 1:232729/1‑150 (MQ=255) gcAAGGCCTTC‑TCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCt > 2:248569/2‑150 (MQ=255) aaGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgtgt > 3:445559/1‑76 (MQ=255) aaGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgtgt < 4:445559/76‑1 (MQ=255) aaGGCCTTC‑TCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCg < 1:644599/103‑1 (MQ=255) aaGGCCTTC‑TCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCg > 2:644599/1‑103 (MQ=255) aGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGCACCATCTGACCCTGGc < 3:363586/150‑1 (MQ=255) aGGCCTTC‑TCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGc < 3:459222/150‑1 (MQ=255) ggCCTTC‑TCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGcc < 1:407449/130‑1 (MQ=255) ggCCTTC‑TCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGcc > 2:407449/1‑130 (MQ=255) c‑tcCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc > 3:18664/1‑117 (MQ=255) c‑tcCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc < 4:18664/117‑1 (MQ=255) tcCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCa > 2:637442/1‑150 (MQ=255) cTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCAcc > 1:270125/1‑150 (MQ=255) tGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCgg > 1:541068/1‑71 (MQ=255) tGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCgg < 2:541068/71‑1 (MQ=255) | GCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGA‑AAACACTCAAGGCCTTC‑TCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACC > NC_000913_3_pae_tpiA/3959366‑3959663 |
| 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: ‑ |