| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | ACB122 | 1,974,160 | 0 | T | G | 52.2% | 12.2 / 22.3 | 23 | intergenic (+42/‑258) | rpsA/ihfB | 30S ribosomal protein S1/Integration host factor subunit beta |
| Reads supporting (aligned to +/- strand): ref base T (5/6); new base G (12/0); total (17/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.58e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTA > ACB122/1974026‑1974289 | gAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCt > 3:129578/1‑135 (MQ=255) agagCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCt > 7:436428/1‑135 (MQ=255) agagCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCGGTTGTTTg > 1:161073/1‑133 (MQ=255) cttaaaaaCGCTCCGTAAGCTGCTGCCGACACCACCTTGGCTGCGCTGCTGCGAGATGCAATGGCCATGCACAACTGAGTTGTGTTTGATGGTTAAAAAGGGCCACCCTCGGGTCGCCCTTTTTTTTGCTCcgcg < 4:124088/132‑1 (MQ=255) aaaacCCCTCCGGACGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGAcc < 8:156209/135‑1 (MQ=255) aaaCGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCGTTTTTTTGCTCCGCGCCCAg > 7:184570/1‑135 (MQ=255) aaaCGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCGTTTTTTTGCTCCGCGACCac > 5:95850/1‑134 (MQ=255) gAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTgcggc < 8:436428/135‑1 (MQ=255) gcggctgcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTGCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGa < 4:242475/132‑1 (MQ=255) gcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCGCTTTTTTTTGCTCCGCGCCCAGGTGGTGCGGCGGATGGCAc > 2:365317/1‑135 (MQ=255) gcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCGCGGTTTTTTTCTACGCGCCCCGGTGGTGGGGGGGGTGGccc > 3:52563/1‑133 (MQ=255) gcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCGTTTTTTTGGTCCGCGACCAGGTGGTGGGGCGGATGGccc > 2:98898/1‑133 (MQ=255) caccaTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGA‑CCTCGGGGCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAac < 4:313224/135‑1 (MQ=255) caTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGTTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAacac < 2:405018/135‑1 (MQ=255) ctgctgCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCCCTGCTGGACCACGCCCGCGTCa > 6:372414/1‑135 (MQ=255) tgctgcGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCGGTTTTTTTGCTCCGCGACCAGGTGGTGGGGCGGCTGGCGCGGCTGGTACACGCCCGCGTCAg > 7:480484/1‑135 (MQ=255) tgctgcGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGCCCCGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTGAg > 2:71575/1‑135 (MQ=255) cgcgAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCGCGTTTTTTTGCTCCGCGCCCAGGTGGTGCGGCGGATGGGCCTGCTGGAACACGCCCGCTTCCGGCTgg > 2:316186/1‑135 (MQ=255) tGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCGTTTTTTTGCTCCGCGCCCAGGTGGTGCGGCGGATGGATCTGCTTGAACACGCCACCTTAAGGCTGGCGAAGcccct > 8:515796/1‑133 (MQ=255) tGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCGATTTTTTGCACCGCCACCACGTGCAGCCGGCGATGGCACTGCTGGAACACGCCAGCTGAAGGCTGGCGAAGCCCga > 8:523697/1‑134 (MQ=255) ccAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGCCCAGGTGGTGCGGCGGCTGGCACTGCTGGAACACCCCACCGTCAGGCTGGCGCAGCCCCtt > 3:239265/1‑133 (MQ=255) cAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCGCGGTTTTTTGCTCCGCGCCCAGGGGGTGCGGCGGCTGGGCCGGCTGGAACACGCCAGCTTCAGGCGGGCGGAGCCCGGTGTAgggg > 4:515871/1‑133 (MQ=255) tCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCGCGTTTTTTTTGTCCGCGCCCCGGGGGTGGGGGGGCTGGCACGGGGGGAACACCCCCGCGTCAGGCGGGCGCAGGCCGGTGTTGGCGCGGGttgcccc > 6:312831/1‑130 (MQ=255) ccTCGGGTCGCCCGTTTTTTTGCTCCGCGCCCAGGTGGTGCGGCGGATGGCACTGCGGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGCAGCAGGCGCCCCTGGCCTTTc > 8:378670/1‑134 (MQ=255) ggctgtttttttgGCCCCGCGACCGGGTGGGGGGGGGGTGGCCAGTGGGGCACCCCGCCAGGTTCCGGCTGGCCAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTa < 5:372414/130‑1 (MQ=19) | GAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTA > ACB122/1974026‑1974289 |
| 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: ‑ |