Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 596,346 | 0 | C | T | 55.2% | 9.2 / 31.2 | 29 | Y249Y (TAC→TAT) | cusC | copper/silver export system outer membrane channel |
Reads supporting (aligned to +/- strand): ref base C (8/5); new base T (7/9); total (15/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.62e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCAGCGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGCGGGATTGTCGTCG > NC_000913/596273‑596421 | gCAGCGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGc > 2:129978/1‑81 (MQ=255) cAGCGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCt < 1:129951/81‑1 (MQ=255) cAGCGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCt > 2:129997/1‑81 (MQ=255) cGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTATGGCAAGCTGcc < 2:129952/81‑1 (MQ=255) gACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTATGGCAAGCTGCCg > 1:129976/1‑81 (MQ=255) gCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTATGGCAAGCTGCCGCAAgcg > 1:129981/1‑81 (MQ=255) aGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTATGGCAAGCTGCCGCAAGCGCAGACAGTa > 1:129988/1‑80 (MQ=255) aGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTATGGCAAGCTGCCGCAAGCGCAGACAGTa > 1:129982/1‑80 (MQ=255) aGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTa > 2:129991/1‑80 (MQ=255) ggAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAg > 1:129985/1‑81 (MQ=255) gAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTATGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGc < 1:129953/81‑1 (MQ=255) cTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTATGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGAc < 2:129954/81‑1 (MQ=255) gcgcAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGc > 2:129994/1‑81 (MQ=255) cgcAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGcc < 1:129955/81‑1 (MQ=255) cAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTATGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTg < 2:129957/81‑1 (MQ=255) cAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTATGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTg > 2:129986/1‑81 (MQ=255) cAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTATGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTg < 1:129956/81‑1 (MQ=255) cGAATAATGCATTGCAACTGTTGCTGGGAAGCTATGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCa > 2:504370/1‑79 (MQ=255) aatGCATTGCAACTGTTGCTGGGAAGCTATGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGt > 2:129995/1‑80 (MQ=255) tGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTaa < 1:129958/81‑1 (MQ=255) cAACTGTTGCTGGGAAGCTATGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACcg < 2:129959/81‑1 (MQ=255) aCTGTTGCTGGGAAGCTATGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACcggc < 2:129960/81‑1 (MQ=255) ttGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGCggg < 1:129961/80‑1 (MQ=255) tGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGCGGGAt > 2:130004/1‑81 (MQ=255) cTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGCGGGAtt < 1:129962/80‑1 (MQ=255) ggAAGCTATGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGCGGGATTgt < 2:129963/79‑1 (MQ=255) ggAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGCGGGATTgtcg > 1:129993/1‑81 (MQ=255) aaGCTATGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGCGGGATTgtcgtc < 2:129964/81‑1 (MQ=255) gCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGCGGGATTGtcgtcg > 1:129998/1‑80 (MQ=255) | GCAGCGACATCGCTAAACGTCAGGGGGAACTGGCGCAGGCGAATAATGCATTGCAACTGTTGCTGGGAAGCTACGGCAAGCTGCCGCAAGCGCAGACAGTAAACAGCGACAGCCTGCAAAGCGTTAAATTACCGGCGGGATTGTCGTCG > NC_000913/596273‑596421 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |