| Read alignment evidence... | ||||||||||
|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | change | freq | score | reads | annotation | genes | product | ||
| * | NC_000913 | 2,785,196 | 1 | .→A | 100.0% | 2.4 | 4 | intergenic (‑185/+566) | ypjC/ileY | pseudogene/tRNA‑Ile |
| Rejected: E-value exceeds prediction threshold. | ||||||||||
| Reads supporting (aligned to +/- strand): new base (1/1): ref base (1/1): total (2/2) | ||||||||||
| Fisher's exact test strand distribution p-value = 1.00e+00 | ||||||||||
| Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.07e-01 | ||||||||||
ATTAATGTCTACATGGCTATTTTTAATGTTATTACTGTTTGTCACTATAAAAAAACGCTCATTTGAGACAATTACTGACATTAACTGCTTCACTTGCTACGCATGGAACTTTTAATTAAATTAGCACAGGAATGTTAAATTTAATAAACAAAAGGTTATTTCGCTGTATGAT‑AAAAAAAACCGTTATAATTTATTAGTGAAAATCGTTTTTCAAGTGTTAGAAATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTACTCAC > NC_000913/2785025‑2785419 | aTTAATGTCTACATGGCTATTTTTAATGTTATTACTGTTTGTCACTATAAAAAAACGCTCATTTGAGACAATTACTGACATTAACTGCTTCACTTGCTACACATGGAACTTTTAATTAAATTAGCACAGGAATGTTAAATTTAATAAACAAAAGGTTATTTCGCTGTATGAT‑AAAAAAAACCGTTATAATTTATTAGTGAAAATCGTTTTTCAAGTGTTAGAAATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGcggcgg > 1:15779/1‑301 (MQ=255) ttGAGACAATTACTGACATTAACTGCTTCACTTGCTACACATGGAACTTTTAATTAAATTAGCACAGGAATGTTAAATTTAATAAACAAAAGGTTATTTCGCTGTATGAT‑AAAAAAAACCGTTATAATTTATTAGTGAAAATCGTTTTTCAAGTGTTAGAAATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGc < 2:15779/301‑1 (MQ=255) gACATTAACTGCTTCACTTGCTACGCATGGAACTTTTAATTAAATTAGCACAGGAATGTTAAATTTAATAAACAAAAGGTTATTTCGCTGTATGATAAAAAAAAACCGTTATAATTTATTAGTGAAAATCGTTTTTCAAGTGTTAGAAATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCtacata > 1:303231/1‑301 (MQ=255) cTACGCATGGAACTTTTAATTAAATTAGCACAGGAATGTTAAATTTAATAAACAAAAGGTTATTTCGCTGTATGATAAAAAAAAACCGTTATAATTTATTAGTGAAAATCGTTTTTCAAGTGTTAGAAATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTACTCAc < 2:303236/300‑1 (MQ=255) | ATTAATGTCTACATGGCTATTTTTAATGTTATTACTGTTTGTCACTATAAAAAAACGCTCATTTGAGACAATTACTGACATTAACTGCTTCACTTGCTACGCATGGAACTTTTAATTAAATTAGCACAGGAATGTTAAATTTAATAAACAAAAGGTTATTTCGCTGTATGAT‑AAAAAAAACCGTTATAATTTATTAGTGAAAATCGTTTTTCAAGTGTTAGAAATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTACTCAC > NC_000913/2785025‑2785419 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
| Reads not counted as support for junction |
| read_name Not counted due to insufficient overlap past the breakpoint. |
| read_name Not counted due to not crossing MOB target site duplication. |