| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 566,597 | 0 | T | C | 52.4% | 2.8 / 17.2 | 22 | intergenic (‑619/‑87) | intD/renD | putative integrase/protein RenD |
| Reads supporting (aligned to +/- strand): ref base T (5/5); new base C (4/7); total (9/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.70e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑T‑GACCTGTATGACAGCGCGAATTAACCGTG‑GTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTA > NC_000913/566363‑566773 | gtttgatatggtttACGAGTATTGCCGGACGCGAGGCCTGTATCCGGATGAGGCGTCCTATCCGTGGAAATCAAACGAGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCC‑ATATGACTGCGAGAATTAACCGTG‑GTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATgg > 1:231847/14‑239 (MQ=255) ggtttACGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCC‑ATATGACTGCGAGAATTAACCGTG‑GTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTg > 1:44783/5‑251 (MQ=255) gcgaggCCTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCC‑ATATGACTGCGAGAATTAACCGTG‑GTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGATCTCTAAATCATGCACAGGCTCTGGCGAAGATCg < 1:202944/247‑1 (MQ=255) gTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATAGAACGCGCATTAATGGTTGGTTCCCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGATGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑T‑GACCTGTATGACAGCGCGAATTAACCGTG‑GTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGCAGATCGc > 1:173603/1‑250 (MQ=255) cTGTATCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCC‑ATATGACTGCGAGAATTAACCGTG‑GTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTg < 1:84083/223‑1 (MQ=255) tCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGAGGCGGAAGTACGCCGTAAGGCCGCAGAGGAGC‑TTGTCC‑ATATGACTGCGAGAATGAACCGTG‑GTTAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGc < 1:107010/186‑1 (MQ=255) ggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCGT‑GTCC‑ATATGACTGCGAGAATTAACCGTG‑GTGAGGCGAGCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCggn < 1:29953/154‑2 (MQ=255) gttACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑T‑GACCTGTATGACAGCGCGAATTAACCGTG‑GTGAGACGATACCTGAACCAGGAAAACAACTTCCTGTTATGGGCGGTan < 1:69281/155‑2 (MQ=255) ccAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTTCCGATGAAC‑T‑GACCTGTATGACAGCGCGAATTAACCGTG‑GTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTggg < 1:134333/251‑1 (MQ=255) aGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑T‑GACCTGTATGACAGCGCGAATTAACCGTG‑GTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCGCTGGCGAAGATCGCAGAAAGTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATGACCTggg < 1:45379/249‑1 (MQ=255) ggCCAATGCGCTTACTTATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCC‑ATATGACTGCGAGAATTAACCGTG‑GTGAGGCGATCCCTGANCCANTAAANCANCGTCCTGTCATGGGCNGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAANNCNNNGNTGAGTTCg < 2:180356/185‑1 (MQ=255) ggCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCC‑ATATGACTGCGAGAATTAACCGTG‑GTGAGGCGATACCTGAACCAGTAAAACACCTTCCTGTCATGGGCGGAAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCg > 1:180356/1‑185 (MQ=255) ggCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAA‑T‑GACCTGTATGACAGCGCGAATTAACCGCG‑GTGAGAAGATCCCTGAACCAGTAAAACACCTTCCTGTTATGGGCGGTAGCCCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCCGACATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAag > 1:131426/1‑249 (MQ=255) aGGCTGCCGATGAAC‑T‑GACCTGTATGACAGCGCGAATTAACCGTG‑GTGAGACGATACCGGATCCAGTAAAACAACTTCCTGTGATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGTACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGgct < 1:1154/251‑2 (MQ=255) gCTGCCGATGAAC‑T‑GACCTGTATGACAGCGCGAATTAACCGTG‑GTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCa > 1:135971/1‑185 (MQ=255) cGATGAAC‑T‑GACCTGTATGACAGCGCGAATTAACCGTG‑GTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCg > 1:72354/1‑250 (MQ=255) gatgagcttatccatattactgcgagAATTAACCGGG‑GTGAGGCGATCCCGGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAAGGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTTCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTa < 1:78509/226‑1 (MQ=255) gTATGACAGCGCGAATTAACCGTG‑GTGAGACGATACCTGAACCAGTAAAACAACTTCCAGTTATGGGCGGTAGAACTCTAAATCGTGcac > 1:182511/1‑90 (MQ=255) gNATGNCAGCNCGNATTAACC‑TGTGTGAGNCGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGANNTNNNANTCGTGcac < 2:182511/90‑1 (MQ=255) tATGACAGCGCGAATTAACCGTG‑GTGAGACGATACCTGAACCAGTAAAACAACTGCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGc < 1:126827/93‑1 (MQ=255) ttAACCGTG‑GTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAANNTCGCAg > 2:19767/1‑96 (MQ=255) ttAACCGTG‑GTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTg < 1:19767/200‑1 (MQ=255) | AAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑T‑GACCTGTATGACAGCGCGAATTAACCGTG‑GTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTA > NC_000913/566363‑566773 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |