| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 710,914 | 0 | A | T | 35.0% | 16.6 / 16.7 | 20 | intergenic (‑189/+21) | uof/fldA | ryhB‑regulated fur leader peptide/flavodoxin 1 |
| Reads supporting (aligned to +/- strand): ref base A (6/7); new base T (3/4); total (9/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCAAATATTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGA > NC_000913/710766‑711043 | gCAAATATTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACa < 4:120096/149‑1 (MQ=255) aataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACTCATGTCTAcgc < 4:65241/148‑1 (MQ=255) taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcgccg > 3:62750/1‑148 (MQ=255) aaaaaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGca > 6:100105/1‑149 (MQ=255) aaaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACTCATGTCTACGCCGtat > 3:140035/1‑148 (MQ=255) aaaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGt < 6:131101/137‑1 (MQ=255) gTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACTCATGTCTACGCCGtattaata < 4:140035/149‑6 (MQ=255) aCATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTga < 3:75570/147‑1 (MQ=255) aCATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTga < 5:100105/147‑1 (MQ=255) aaGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATg > 4:31268/1‑149 (MQ=255) cACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACTCATGTCTAcgc > 3:124064/1‑108 (MQ=255) cACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACTCATGTCTAcgc < 4:124064/108‑1 (MQ=255) aaCACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCa < 3:31268/148‑1 (MQ=255) ttGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTaa < 4:62750/149‑1 (MQ=255) aacaGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTa < 1:133562/149‑1 (MQ=255) aGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTacac > 3:118929/1‑149 (MQ=255) gTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTcc < 4:118929/148‑1 (MQ=255) aTTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTtcgt > 2:70259/1‑149 (MQ=255) cTATTAAAACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCAATGAGAATTTCGTAGAGATGCAACTCTTCTGAAATCTGATTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGa > 8:63524/1‑149 (MQ=255) tgcGGCGTAGACTCATGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGa > 5:39600/3‑121 (MQ=255) tgcGGCGTAGACTCATGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGa < 6:39600/119‑1 (MQ=255) aTGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACg > 5:46190/1‑111 (MQ=255) aTGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACg < 6:46190/111‑1 (MQ=255) | GCAAATATTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGA > NC_000913/710766‑711043 |
| 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: ‑ |