| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913_3_pae_pgi | 1,736,097 | 0 | T | A | 28.6% | 22.5 / 15.6 | 21 | intergenic (+138/+24) | sodB/ydhP | superoxide dismutase, Fe/putative MFS transporter, inner membrane protein |
| Reads supporting (aligned to +/- strand): ref base T (9/6); new base A (3/3); 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 = 3.86e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATC > NC_000913_3_pae_pgi/1735949‑1736230 | tcGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACCTTGTCAGCAACTGTCACGCAGAAGGtt > 1:6203/1‑150 (MQ=255) cGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTa > 2:80759/1‑150 (MQ=255) cGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTa > 2:80745/1‑150 (MQ=255) tGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATcc < 2:140950/150‑1 (MQ=255) tGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGCCAAGATAAATCGCACGGTGTCAGCAACTGTAACGCAGAAGGTTATcc > 2:145100/1‑150 (MQ=255) cATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCtt > 2:85334/1‑150 (MQ=255) ataaCTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGc < 1:105111/150‑1 (MQ=255) tGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAAc > 1:7855/1‑150 (MQ=255) gCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGttt > 1:46059/1‑150 (MQ=255) tACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGc < 1:80670/100‑1 (MQ=255) tACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGc > 2:80670/1‑100 (MQ=255) gATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAGTCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTATCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTg < 1:104359/150‑1 (MQ=255) tCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCa < 1:85334/150‑1 (MQ=255) tCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCAcc > 4:110574/1‑150 (MQ=255) cACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCTGc < 3:25061/150‑1 (MQ=255) aCGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTg > 2:19516/1‑92 (MQ=255) aCGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTg < 1:19516/92‑1 (MQ=255) acaaaaaCGCAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATc > 2:118000/6‑150 (MQ=255) aaCGCAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTgg < 3:57081/76‑1 (MQ=38) aaCGCAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTgg > 4:57081/1‑76 (MQ=38) cGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGcaca > 2:197650/1‑133 (MQ=255) cGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGcaca < 1:197650/133‑1 (MQ=255) | TCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCGACAATCGCTCCCATCACCGGCACAAAGCTGTATC > NC_000913_3_pae_pgi/1735949‑1736230 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |