| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | ACB122 | 5,148,290 | 0 | A | C | 57.1% | 9.8 / 20.6 | 21 | intergenic (‑110/+146) | bioP/azoR | Biotin transporter/FMN‑dependent NADH‑azoreductase |
| Reads supporting (aligned to +/- strand): ref base A (1/8); new base C (12/0); total (13/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.42e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.98e-03 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCGCGACGATTATCAGATAGCCCATGGTTGCCTCTCGTGAATCAAGGCCGCGACCTTAGCGGGTTGCAGCGAATGAGTTCAACCATGAGGAATGTAGGGGCTGTACCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCTGCATATCATCCGGCGAGCGAATAGGCCGGCACTGACCACCAAAAAACCCCAGGCATAAAAAAGCCCGGCCATACAAGCCGGGCAGGTACACCCAAAGGAGCAAAAACAGGTGTCAGGTTTGGGAA > ACB122/5148156‑5148424 | gcgcGACGATTATCAGATAGCCCATGGTGGCCTCTCGTGAATCAAGGCCGCGACCTTAGCGGGTTGCAGCGAATGAGTTAAACCATGAGGAATGTAGGGGCTGTACCGGCCTCTTCGCGGGCACGCCCGCTCCcac < 2:208962/136‑1 (MQ=255) gaTTATCAGATAGCCCATGGTTGCCTCTCGTGAATCAAGGCCGCGACCTTAGCGGGTTGCAGCGAATGAGTTCAACCATGAGGAATGTAGGGGCTGTACCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGAcc < 2:19555/136‑1 (MQ=255) gCCCATGGTTGCCTCTCGTGAATCAAGGCCGCGACCTTAGCGGGTTGCAGCGAATGAGTTCAACCATGAGGAATGTAGGGGCTGTACCGGCCTCTTCGCGGGCCCGCCCGCTCCCCCAGGGACCGGCCTatcatc > 6:232526/1‑135 (MQ=255) gccAAGGTTGCCTCTCGTGAATCAAGGCCGCGACCTTAGCGGGTTGCAGCGAATGAGTTCAACCATGAGGAATGTAGGGGCTGTACCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCTGCATATCATCCg < 1:128784/135‑1 (MQ=255) ccATGGTTGCCTCTCGTGAATCAAGGCCGCGACCTTAGCGGGTTGCAGCGAATGAGTTCAACCATGAGGAATGTAGGGGCTGTACCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCTGCATATCATCCgg < 2:217751/136‑1 (MQ=255) cATGGTTGCCTCTCGTGAATCAAGGCCGCGACCTTAGCGGGTTGCAGCGAATGAGTTCAACCATGAGGAATGTAGGGGCTGTACCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCTGCATATCATCCgg < 6:222150/135‑1 (MQ=255) ttGCCTCTCGTGAATCAAGGCCGCGACCTTAGCGGGTTGCAGCGAATGAGTTCAACCATGAGGAATGTAGGGGCTGTACCGGCCTCTTCGCGGGGCCGCCCGCGCCCCCACGGCCCAGCCTATCATCCGgcgcgcg > 1:22627/1‑136 (MQ=255) ctcGTGAATCAAGGCCGCGACCTTAGCGGGTTGCAGCGAATGAGTTCAACCATGAGGAATGTAGGGGCTGTACCGGCCTCTTCGCGGGCACGCCCGCTCCCCCAGGGCCCTGCATATCATACGGCGAGCGAATAgg > 8:51425/1‑136 (MQ=255) ctcGTGAATCAAGGCCGCGACCTTAGCGGGTTGCAGCGAATGAGTTCAACCATGAGGAATGTAGGGGCTGTACCGGCCTCGTCGCGGGCACGCCCGCGCCCCCAGGGGCCGGGCTCTCATCCcgcgcgcgcggtgg > 7:225508/1‑130 (MQ=255) gAATCAAGGCCGCGACCTTAGCGGGTTGCAGCGAATGAGTTCAACCATGAGGAATGTAGGGGCTGTACCGGCCTCTTAGGGGGCCCCCCCGCGGCCCCACGGGGCGGGCTATCATCCCGCGAGGGATAATGCCGGc > 1:141676/1‑136 (MQ=255) tCAAGGCCGCGACCTTAGCGGGTTGCAGCGAATGAGTTCAACCATGAGGAATGTAGGGGCTGTACCGGCCTCTTCTCCGGCCCGCCCGCACCCCCAGGGACCGGCAGATAATCCCGCGCGCGAAAATGCCCGccag > 7:105758/1‑133 (MQ=255) ccGCGACCTTAGCGGGTTGCAGCGAATGAGTTCAACCATGAGGAATGTAGGGGCTGTACCGGCCTCTTCGCGGGGCCGCCCCCTCCCCCCGGGTCCTGCCTATACTCCGGCGCGGGCATAGGCCGGCACGGccca > 5:29643/1‑135 (MQ=255) cAGCGAATGAGTTCAACCATGAGGAATGTAGGGGCTGTACCGGCCTCTTCGCGGGCACGCCCGCTCCCCCAGGGACCTGCATATCATCCGGCGAGCGAATAGGCCGGCACTGACCACCAAAAAACCCCAGGCATa > 8:51482/1‑135 (MQ=255) gctttCAACCATTAGGATTGTAGTGGCTGTACCGGCCTCTTCGCGGGCACGCCCTCTCCCACAGGGACCTGCATATCATCCGGCGAGCGAATAGGCCGGCACTGACCACCAAAAAACCCCAGGCATAAAAAAGcc < 6:61039/132‑1 (MQ=39) gAGTTCTAACATGAGGAATGTTGCGGCTGAACAGGCATGTTCGCCGGCACGCCCGCTCCCACAGGGACCTGCATATGATCCGGCGAGCGAATAGGCCGGCACTGACCACCAAAAAACCCCAGGCATAAAAAAGcc < 6:57378/135‑1 (MQ=255) aaTGTAGGGGCTGTACCGGCCTCTTAGCGGGCACGCCCGCTCCCACACGGGCCGGCATATCAACCGGCGAGCGAATAGGCCGGCCCTGCCCACCAAAAAACCCCAGGGATAAAAAAGCCCGGCCCTACAAGCCgg > 6:102450/1‑135 (MQ=38) gTACCGGCCTCTTCGCGGGCACGCCCGCTCCCCCAGGGGCCGGCATATCATCCGGCGCGCGAATAGGCCGGCCCTGCCCACCAAAAAACCCCAGGCATAAAAAAGCCCGGCCATACAAGCCGGGCAGGTACAccc > 3:66533/1‑135 (MQ=39) ccGGCCTCTTCGCGGGCCCGCCCGCTCCCCCAGGGACCTGCATATCATCCGGCGCGCGAATAGGCCGGCCCGGCCCACCAAAAAACCCCAGGCAGAAAAAAGCCCGGCCCTACAAGCCGGGGAGGTACACCCAAAg > 7:156234/1‑136 (MQ=255) cGGCCTCTTCGCGGGGCCGCCCGCTCCCCCCGGGGCCTGGCTTTAATCCGGCGCGCGAATAGGCCGGCCCTGCCCCCCAAAAAACCCCCGGCATAAAAAAGCCCGGCCAGACAAGCCGGGCAGGTACACCCAAAgg > 8:221374/1‑136 (MQ=255) gCACGCCCGCTCCCACAGGGACCTGCATATCATCCGGCGAGCGAATAGGCCGGCACTGACCACCAAAAAACCCCAGGCATAAAAAAGCCCGGCCATACAAGCCGGGCAGGTACACCCAAAGGAGCAAAAACAGgt < 5:102450/135‑1 (MQ=255) gCTCCCCCAGGGACCTGCATATCATCCGGCGAGCGAATAGGCCGGCACTGACCACCAAAAAACCCCAGGCATAAAAAAGCCCGGCCATACAAGCCGGGCAGGTACACCCAAAGGAGCAAAAACAGGTGTCAGGtt > 6:85771/1‑135 (MQ=255) ccACAGGGTCCTGCATATCATCCGGCGAGCGAATAGGCCGGCACTGACCACCAAAAAACCCCAGGCATAAAAAAGCCCGGCCATACAAGCCGGGCAGGTACACCCAAAGGAGCAAAAACAGGTGTCAGGTTTgg < 8:105758/134‑1 (MQ=255) cacaGGGACCTGCATATCATCCGGCGAGCGAATAGGCCGGCACTGCCCCCCAAAAAACCCCAGGCAGAAAAAAGCCCGGCCGTACAAGCCGGGCAGGTACCCCCAAAGGAGCAAAAGCAGGTGTCAGGTTTGGGaa > 1:18669/1‑136 (MQ=255) cacaGGGACCTGCATATCATCCGGCGAGCGAATAGGCCGGCACTGACCACCAAAAAACCCCAGGCATAAAAAAGCCCGGCCATACAAGCCGGGCAGGTACACCCAAAGGAGCAAAAACAGGTGTCAGGTTTGGGaa > 3:101169/1‑136 (MQ=255) | GCGCGACGATTATCAGATAGCCCATGGTTGCCTCTCGTGAATCAAGGCCGCGACCTTAGCGGGTTGCAGCGAATGAGTTCAACCATGAGGAATGTAGGGGCTGTACCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCTGCATATCATCCGGCGAGCGAATAGGCCGGCACTGACCACCAAAAAACCCCAGGCATAAAAAAGCCCGGCCATACAAGCCGGGCAGGTACACCCAAAGGAGCAAAAACAGGTGTCAGGTTTGGGAA > ACB122/5148156‑5148424 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |