| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,363,715 | 0 | A | C | 23.8% | 35.9 / 10.1 | 21 | intergenic (‑82/+18) | yfaY/yfaZ | inactive PncC family protein/outer membrane protein, putative porin |
| Reads supporting (aligned to +/- strand): ref base A (9/7); new base C (3/2); total (12/9) | |||||||||||
| 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.08e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGATTTGCCCGTGTAACACTTCATCCCCGGTGGATAACATTTCCACTTTTAACATGGTTATCTCCTGCTTTTGTTATCGAATCACTATAACGCAAAGGGAGATAAGGCGTAGTAAGAAACTGAATGGCGCGACAGAATGCCGCGCCGGAAGGGATCAGAAACTGGCGTTAACCCCAACATACGGGCCGTCAGCCACGGCGTTGTCGCGGTTACCGTCTTTACCCGACAGATTCAGGTAGCGATAACCCGCCTCAATACTGACTGGACGCAT > NC_000913/2363579‑2363849 | cGATTTGCCCGTGTAACACTTCATCCCCGGTGGATAACATTTCCACTTTTAACATGGTTATCTCCTGCTTTTGTTATCGAATCACTATAACGCAAAGGGAGATAAGGCGTAGTAAGAAACTGAATGGCGCGGCATTCTGTCGCGccatt > 1:156793/1‑146 (MQ=255)cGATTTGCCCGTGTAACACTTCATCCCCGGTGGATAACATTTCCACTTTTAACATGGTTATCTCCTGCTTTTGTTATCGAATCACTATAACGCAAAGGGAGATAAGGCGTAGTAAGAAACTGAATGGCGCGGCATTCTGTCGCGccatt < 2:156793/149‑4 (MQ=255) gCCCGTGTAACACTTCATCCCCGGTGGAAAACATTTCAACTTTAAACATGGTTATCTCCTGCTTTTGTTATCGAATCACTATAACGCAAAGGGAGATAAGGCGTAGTAAGAAACTGAATGGTGCGGCATTCTGTCGCGccat > 1:16792/1‑140 (MQ=255) cccGTGTAACACTTCATCCCCGGTGGATAACATTTCCACTTTTAACATGGTTATCTCCTGCTTTTGTTATCGAATCACTATAACGCAAAGGGAGATAAGGCGTAGTAAGAAACTGAATGGCGCGACAGAATGCCgcgc < 1:39977/138‑1 (MQ=255) cccGTGTAACACTTCATCCCCGGTGGATAACATTTCCACTTTTAACATGGTTATCTCCTGCTTTTGTTATCGAATCACTATAACGCAAAGGGAGATAAGGCGTAGTAAGAAACTGAATGGCGCGACAGAATGCCgcgc > 2:39977/1‑138 (MQ=255) aTTTCCACTTTTAACATGGTTATCTCCTGCTTTTGTTATCGAATCACTATAACGCAAAGGGAGATAAGGCGTAGTAAGAAACTGAATGGCGCGGCATTCTGTCGCGcc < 1:22701/108‑1 (MQ=255) aTTTCCACTTTTAACATGGTTATCTCCTGCTTTTGTTATCGAATCACTATAACGCAAAGGGAGATAAGGCGTAGTAAGAAACTGAATGGCGCGGCATTCTGTCGCGcc > 2:22701/1‑108 (MQ=255) cACTTTTAACATGGTTATCTCCTGCTTTTGTTATCGAATCACTATAACGCAAAGGGAGATAAGGCGTAGTAAGAAACTGAATGGCGCGACAGAATGCCGCGCCGGAAGGGATCAGAAACTgg < 8:38086/122‑1 (MQ=255) cACTTTTAACATGGTTATCTCCTGCTTTTGTTATCGAATCACTATAACGCAAAGGGAGATAAGGCGTAGTAAGAAACTGAATGGCGCGACAGAATGCCGCGCCGGAAGGGATCAGAAACTgg > 7:38086/1‑122 (MQ=255) cTTTTGTTATCGAATCACTATAACGCAAAGGGAGATAAGGCGTAGTAAGAAACTGAATGGCGCGACAGAATGCCGCGCCGGAAGGGATCAGAAACTGGCGTTAAc < 1:7752/105‑1 (MQ=255) cTTTTGTTATCGAATCACTATAACGCAAAGGGAGATAAGGCGTAGTAAGAAACTGAATGGCGCGACAGAATGCCGCGCCGGAAGGGATCAGAAACTGGCGTTAAc > 2:7752/1‑105 (MQ=255) cTTTTGTTATCGAATCACTATAACGCAAAGGGAGATAAGGCGTAGTAAGAAACTGAATGGCGCGACAGAATGCCGCGCCGGAAGGGATCAGAAACTGGCGTTAACCCCAACATACGGGCCGTCAGCCACGGCGTTGTCGCGGTTACCGt > 3:131372/1‑149 (MQ=255) tATCGAATCACTATAACGCAAAGGGAGATAAGGCGTAGTAAGAAACTGAATGGCGCGACAGAATGCCGCGCCGGAAGGGATCAGAAACTGGCGTTAACCCCAACATACGGGCCGTCAGCCACGGCGTTGTCGCGGTTACCGTCTTTAcc > 8:70191/1‑149 (MQ=255) aaaGGGAGATAAGGCGTAGTAAGAAACTGAATGGCGCGACAGAATGCCGCGCCGGAAGGGATCAGAAACTGGCGTTAACCCCAACATACGGGCCGTCAGCCACGGCGTTGTCGCGGTTACCGTCTTTACCCGACAGATTCAGGTAGCGa > 6:16677/1‑149 (MQ=255) aaGGGAGATAAGGCGTAGTAAGAAACTGAATGGCGCGACAGAATGCCGCGCCGGAAGGGATCAGAAACTGGCGTTAACCCCAACATACGGGCCGTCAGCCACGGCGTTGTCGCGGTTACCGTCTTTACCCGAc > 5:137080/1‑133 (MQ=255) aaGGGAGATAAGGCGTAGTAAGAAACTGAATGGCGCGACAGAATGCCGCGCCGGAAGGGATCAGAAACTGGCGTTAACCCCAACATACGGGCCGTCAGCCACGGCGTTGTCGCGGTTACCGTCTTTACCCGAc < 6:137080/133‑1 (MQ=255) gagaTAAGGCGTAGTAAGAAACTGAATGGCGCGACAGAATGCCGCGCCGGAAGGGATCAGAAACTGGCGTTAACCCCAACATACGGGCCGTCAGCCACGGCGTTGTCGCGGTTACCGTCTTTACCCGACAGATTCAGGTAGCGATAAcc < 5:16677/149‑1 (MQ=255) cTGAATGGCGCGACAGAATGCCGCGCCGGAAGGGATCAGAAACTGGCGTTAACCCCAACATACGGGCCGTCAGCCACGGCGTTGTCGCGGTTACCGTCTTTACCCGACAGATTCAGGTAGCGATAACCCGCCTCAATACTGACTGGACg > 2:67249/1‑149 (MQ=255) aaTGGCGCGACAGAATGCCGCGCCGGAAGGGATCAGAAACTGGCGTTAACCCCAACATACGGGCCGTCAGCCACGGCGTTGTCGCGGTTACCGTCTTTACCCGACAGATTCAGGTAGCGATAACCCGCCTCAATACTGACTGGACGCAt < 7:70191/149‑1 (MQ=255) aCAGAATGCCGCGCCGGAAGGGATCAGAAACTGGCGTTAACCCCAACATACGGGCCATCAGCCACGGCGTTGTCGCGGTTACCGTCTTTACCCGACAGATTCAGGTAGCGATAACCCGCCTCAATACTGACTGGAc > 2:138957/1‑136 (MQ=255) aCAGAATGCCGCGCCGGAAGGGATCAGAAACTGGCGTTAACCCCAACATACGGGCCATCAGCCACGGCGTTGTCGCGGTTACCGTCTTTACCCGACAGATTCAGGTAGCGATAACCCGCCTCAATACTGACTGGAc < 1:138957/136‑1 (MQ=255) | CGATTTGCCCGTGTAACACTTCATCCCCGGTGGATAACATTTCCACTTTTAACATGGTTATCTCCTGCTTTTGTTATCGAATCACTATAACGCAAAGGGAGATAAGGCGTAGTAAGAAACTGAATGGCGCGACAGAATGCCGCGCCGGAAGGGATCAGAAACTGGCGTTAACCCCAACATACGGGCCGTCAGCCACGGCGTTGTCGCGGTTACCGTCTTTACCCGACAGATTCAGGTAGCGATAACCCGCCTCAATACTGACTGGACGCAT > NC_000913/2363579‑2363849 |
| 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: ‑ |