| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,470,108 | 0 | A | C | 20.0% | 48.7 / 10.4 | 25 | intergenic (‑255/+37) | chiA/tufA | periplasmic endochitinase/translation elongation factor EF‑Tu 1 |
| Reads supporting (aligned to +/- strand): ref base A (10/10); new base C (3/2); total (13/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AAGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGG > NC_000913/3469966‑3470244 | aaGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGcc < 5:45282/149‑1 (MQ=255) tGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCTTCGGCGCCCtttt > 2:52353/1‑149 (MQ=255) aCAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCTTCGGCGCCCttttt < 1:52353/148‑2 (MQ=255) aCAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCtttt < 8:64405/147‑1 (MQ=255) gAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATCAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTttaaaataa > 8:95044/1‑144 (MQ=255) aGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCaaaa > 5:81067/1‑149 (MQ=255) aTTAATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGc < 7:95044/149‑1 (MQ=255) aaTGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAg < 6:81067/149‑1 (MQ=255) gTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGcaa > 6:139130/1‑149 (MQ=255) tAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGcaac > 6:41361/1‑148 (MQ=255) ttttAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCaaaa > 1:41565/1‑97 (MQ=255) ttttAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCaaaa < 2:41565/97‑1 (MQ=255) tttATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCAc < 5:139130/149‑1 (MQ=255) cGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGAtt < 5:41361/148‑1 (MQ=255) aCAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGgccgcc < 2:147041/131‑1 (MQ=255) aCAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGgccgcc > 1:147041/1‑131 (MQ=255) aaaTGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTg > 2:28560/1‑149 (MQ=255) gCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGAt > 2:138773/1‑148 (MQ=255) tgaaaaGGGCGCTTCGGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGAt > 3:32919/3‑147 (MQ=255) aaaaaGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACg > 4:100234/1‑87 (MQ=255) aaaaaGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACg < 3:100234/87‑1 (MQ=255) aaaaGGGCGCTTCGGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAgg < 4:32919/149‑1 (MQ=255) aaaaGGGCGCTACGGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAgg > 4:18041/1‑149 (MQ=255) aaaaGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTa > 5:62580/1‑52 (MQ=255) aaaaGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTa < 6:62580/52‑1 (MQ=255) | AAGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGG > NC_000913/3469966‑3470244 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |