| Read alignment evidence... | ||||||||||
|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | change | freq | score | reads | annotation | genes | product | ||
| * | NC_000913 | 1,212,080 | 1 | .→C | 100.0% | 2.4 | 4 | intergenic (‑77/+623) | iraM/ycgX | RpoS stabilzer during Mg starvation, anti‑RssB factor/DUF1398 family protein |
| 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 | ||||||||||
GGATAATACCAGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGC > NC_000913/1211850‑1212279 | ggATAATACCAGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGACTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑CCCCCCCCTCAACAAACTAAATTAATTAAAATGATGGCTTCTATAAAATAAAATTTAAATCAAGGAAt > 2:1402047/1‑299 (MQ=255) gATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATCTAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTtgtg > 2:346588/1‑229 (MQ=255) gATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTtgtg < 1:346587/229‑1 (MQ=255) ccGTTTTAATTATTATACAATTTATTTGTTATGTGCTATCCCACGGTATTTAATATGGTTAATTAGGAGGGTTATATCTTTAGTATGCATAGGAGTATATTA‑CCCCCCCCTCAAAAAAATAAATTAATGAAAATGATGGCTTATAGAAAATAAAATTTAAAGCAATGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTATCTATAATGTTATACATTCAAATCACTGAACATCCTGTAATTAAAACATAGCCTTTATGc < 1:1402018/301‑1 (MQ=255) | GGATAATACCAGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGC > NC_000913/1211850‑1212279 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 32 ≤ 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. |