| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | AE015451 | 4,741,228 | 1 | . | C | 70.6% | 38.7 / 37.3 | 34 | intergenic (+465/+288) | gltA/yeiW | citrate synthase/putative oxidoreductase |
| Reads supporting (aligned to +/- strand): ref base . (6/4); new base C (14/10); total (20/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.26e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACCTCCGCCAACTGCAATTCGTGAA‑CCAAGCAAAAGGCCCCGAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCC > AE015451/4741104‑4741339 | ttACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgag < 1:4022754/125‑1 (MQ=255) aCGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagc > 1:4022887/1‑123 (MQ=255) ggTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagcagct < 1:4022756/125‑6 (MQ=255) ggTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagcagct > 1:4022940/1‑120 (MQ=255) ggTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagcagc > 1:4022915/1‑120 (MQ=255) gTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagcagctt > 2:4022983/1‑119 (MQ=255) gTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagcagctt > 1:4022827/1‑119 (MQ=255) gTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagcagctt < 2:4022758/125‑7 (MQ=255) aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGgcgc < 1:4022809/125‑1 (MQ=255) aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGcc < 1:4022811/125‑1 (MQ=255) cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCt < 1:4022812/125‑1 (MQ=255) tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCtt > 2:4022947/1‑125 (MQ=255) gggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGcc < 1:4022814/122‑1 (MQ=255) ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTg > 1:4022984/1‑124 (MQ=255) ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGa > 1:4022853/1‑125 (MQ=255) ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGa > 1:4022954/1‑125 (MQ=255) tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc < 1:4022815/125‑1 (MQ=255) ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc > 1:4022955/1‑124 (MQ=255) ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc > 1:4022978/1‑124 (MQ=255) ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa < 1:4022816/125‑1 (MQ=255) ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc > 2:4023027/1‑124 (MQ=255) ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACATCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa < 1:4022817/125‑1 (MQ=255) aGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg > 1:4023097/1‑125 (MQ=255) ggTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACgg > 2:4022925/1‑125 (MQ=255) tCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc > 1:4023041/1‑125 (MQ=255) ccAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc > 2:4022888/1‑124 (MQ=255) cAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcac < 2:4022819/125‑1 (MQ=255) cAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcac > 1:4022914/1‑125 (MQ=255) aGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcaca > 2:4023020/1‑125 (MQ=255) aGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcac > 1:4022991/1‑124 (MQ=255) gCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac > 2:4023105/1‑125 (MQ=255) gCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac > 2:4022957/1‑125 (MQ=255) gCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac > 1:4023257/1‑125 (MQ=255) gCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac > 1:4023070/1‑125 (MQ=255) tctGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcg < 2:4022825/125‑1 (MQ=255) tGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgcg < 1:4022834/125‑1 (MQ=255) gTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGATACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGAt < 1:4022850/125‑1 (MQ=255) cACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg < 2:4022851/125‑1 (MQ=255) aCCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAgg > 2:4023061/1‑125 (MQ=255) ccTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAATCCAGGCAC‑AGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc < 1:4022857/125‑1 (MQ=255) cTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc < 1:4022832/125‑1 (MQ=255) cTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc > 1:4023015/1‑125 (MQ=255) | TTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACCTCCGCCAACTGCAATTCGTGAA‑CCAAGCAAAAGGCCCCGAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCC > AE015451/4741104‑4741339 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |