Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,341,504 | 0 | T | C | 67.7% | 34.4 / 27.1 | 31 | F80S (TTC→TCC) | kdsD | D‑arabinose 5‑phosphate isomerase KdsD |
Reads supporting (aligned to +/- strand): ref base T (5/5); new base C (8/13); total (13/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.01e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGGGGATGGGAAAATCGGGGCATATTGGGCGAAAAATGGCGGCAACGTTTGCCAGCACCGGTACACCTTCATTTTTCGTCCATCCTGGTGAAGCCGCGCATGGTGATTTAGGCATGGTTACCCCACAGGATGTGGTGATTGCTATCTCTAACTCT > NC_000913/3341428‑3341583 | aTGGGGATGGGAAAATCGGGGCATATTGGGCGAAAAATGGCGGCAACGTTTGCCAGCACCGGTACACCTTCATTTTTCGt < 2:580058/80‑1 (MQ=255) tGGGGATGGGAAAATCGGGGCATATTGGGCGAAAAATGGCGGCAACGTTTGCCAGCACCGGTACACCTTCATTTTCCGTcc < 1:580059/81‑1 (MQ=255) aTGGGAAAATCGGGGCATATTGGGCGAAAAATGGCGGCAACGTTTGCCAGCACCGGTACACCTTCATTTTCCGTCCATCCt < 2:580060/81‑1 (MQ=255) aTGGGAAAATCGGGGCATATTGGGCGAAAAATGGCGGCAACGTTTGCCAGCACCGGTACACCTTCATTTTCCGTCCATCCt < 2:580061/81‑1 (MQ=255) aTGGGAAAATCGGGGCATATTGGGCGAAAAATGGCGGCAACGTTTGCCAGCACCGGTACACCTTCATTTTCCGTCCATCCt > 1:580109/1‑81 (MQ=255) tGGGAAAATCGGGGCATATTGGGCGAAAAATGGCGGCAACGTTTGCCAGCACCGGTACACCTTCATTTTCCGTCCATCCTg < 1:580062/81‑1 (MQ=255) ggAAAATCGGGGCATATTGGGCGAAAAATGGCGGCAACGTTTGCCAGCACCGGTACACCTTCATTTTCCGTCCATCCTGGt < 2:580063/81‑1 (MQ=255) gAAAATCGGGGCATATTGGGCGAAAAATGGCGGCAACGTTTGCCAGCACCGGTACACCTTCATTTTCCGTCCATCCTGGTg < 1:580064/81‑1 (MQ=255) aaaaTCGGGGCATATTGGGCGAAAAATGGCGGCAACGTTTGCCAGCACCGGTACACCTTCATTTTTCGTCCATCCTGGTg > 1:580104/1‑80 (MQ=255) aaaaTCGGGGCATATTGGGCGAAAAATGGCGGCAACGTTTGCCAGCACCGGTACACCTTCATTTTTCGTCCATCCTGGTg > 2:580112/1‑80 (MQ=255) aaaTCGGGGCATATTGGGCGAAAAATGGCGGCAACGTTTGCCAGCACCGGTACACCTTCATTTTTCGTCCATCCTGGTGaa < 1:580065/81‑1 (MQ=255) cATATTGGGCGAAAAATGGCGGCAACGTTTGCCAGCACCGGTACACCTTCATTTTTCGTCCATCCTGGTGAAGCCGCGCAt < 2:580066/81‑1 (MQ=255) cATATTGGGCGAAAAATGGCGGCAACGTTTGCCAGCACCGGTACACCTTCATTTTCCGTCCATCCTGGTGAAGCCGCGCAt > 2:580108/1‑81 (MQ=255) gggCGAAAAATGGCGGCAACGTTTGCCAGCACCGGTACACCTTCATTTTTCGTCCATCCTGGTGAAGCCGCGCATGGTGAt > 2:580105/1‑81 (MQ=255) gggCGAAAAATGGCGGCAACGTTTGCCAGCACCGGTACACCTTCATTTTCCGTCCATCCTGGTGAAGCCGCGCATGGTGa < 2:580067/80‑1 (MQ=255) gggCGAAAAATGGCGGCAACGTTTGCCAGCACCGGTACACCTTCATTTTCCGTCCATCCTGGTGAAGCCGCGCATGGTGAt > 1:580107/1‑81 (MQ=255) gCGAAAAATGGCGGCAACGTTTGCCAGCACCGGTACACCTTCATTTTCCGTCCATCCTGGTGAAGCCGCGCATGGTGAttt < 2:580068/81‑1 (MQ=255) gcggcAACGTTTGCCAGCACCGGTACACCTTCATTTTTCGTCCATCCTGGTGAAGCCGCGCATGGTGATTTAGGCATGGt > 1:580121/1‑80 (MQ=255) gcggcAACGTTTGCCAGCACCGGTACACCTTCATTTTCCGTCCATCCTGGTGAAGCCGCGCATGGTGATTTAGGCATGGt > 1:580117/1‑80 (MQ=255) gTTTGCCAGCACCGGTACACCTTCATTTTTCGTCCATCCTGGTGAAGCCGCGCATGGTGATTTAGGCATGGTTACCCca > 1:580115/1‑79 (MQ=255) tttGCCAGCACCGGTACACCTTCATTTTCCGTCCATCCTGGTGAAGCCGTGCATGGTGATTTAGGCATGGTTACCCCACAg > 2:580131/1‑81 (MQ=255) aGCACCGGTACACCTTCATTTTCCGTCCATCCTGGTGAAGCCGCGCATGGTGATTTAGGCATGGTTACCCCACAGGAtgtg > 1:580124/1‑81 (MQ=255) cACCGGTACACCTTCATTTTCCGTCCATCCTGGTGAAGCCGCGCATGGTGATTTAGGCATGGTTACCCCACAGGATgtggt < 1:580069/81‑1 (MQ=255) aCCGGTACACCTTCATTTTCCGTCCATCCTGGTGAAGCCGCGCATGGTGATTTAGGCATGGTTACCCCACAGGATgtggt > 1:580113/1‑80 (MQ=255) ggTACACCTTCATTTTCCGTCCATCCTGGTGAAGCCGCGCATGGTGATTTAGGCATGGTTACCCCACAGGATGTGGTGAtt < 2:580070/81‑1 (MQ=255) cacCTTCATTTTTCGTCCATCCTGGTGAAGCCGCGCATGGTGATTTAGGCATGGTTACCCCACAGGATGTGGTGATTGCt < 1:580072/80‑1 (MQ=255) cacCTTCATTTTCCGTCCATCCTGGTGAAGCCGCGCATGGTGATTTAGGCATGGTTACCCCACAGGATGTGGTGATTGCt < 1:580073/80‑1 (MQ=255) cacCTTCATTTTCCGTCCATCCTGGTGAAGCCGCGCATGGTGATTTAGGCATGGTTACCCCACAGGATGTGGTGATTGCt < 1:580071/80‑1 (MQ=255) cTTCATTTTCCGTCCATCCTGGTGATGCCGCGCATGGTGATTTAGGCATGGTTACCCCACAGGATGTGGTGATTGCTAtct < 1:580074/81‑1 (MQ=255) cATTTTCCGTCCATCCTGGTGAAGCCGCGCATGGTGATTTAGGCATGGTTACCCCACAGGATGTGGTGATTGCTATctct > 2:580136/1‑80 (MQ=255) aTTTTTCGTCCATCCTGGTGAAGCCGCGCATGGTGATTTAGGCATGGTTACCCCACAGGATGTGGTGATTGCTATCTCTaa < 2:580075/81‑1 (MQ=255) tttCGTCCATCCTGGTGAAGCCGCGCATGGTGATTTAGGCATGGTTACCCCACAGGATGTGGTGATTGCTATCTCTAAct < 1:580077/80‑1 (MQ=255) tCGTCCATCCTGGTGAAGCCGCGCATGGTGATTTAGGCATGGTTACCCCACAGGATGTGGTGATTGCTATCTCTAActct < 2:580078/80‑1 (MQ=255) | ATGGGGATGGGAAAATCGGGGCATATTGGGCGAAAAATGGCGGCAACGTTTGCCAGCACCGGTACACCTTCATTTTTCGTCCATCCTGGTGAAGCCGCGCATGGTGATTTAGGCATGGTTACCCCACAGGATGTGGTGATTGCTATCTCTAACTCT > NC_000913/3341428‑3341583 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |