| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 4,501,608 | A→G | 100% | K59K (AAA→AAG) | intB → | predicted integrase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 4,501,608 | 0 | A | G | 84.9% | 83.9 / 7.8 | 33 | K59K (AAA→AAG) | intB | predicted integrase |
| Reads supporting (aligned to +/- strand): ref base A (5/0); new base G (11/17); total (16/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.84e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
CGTCGTGACGAGGCGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATT > W3110S.gb/4501555‑4501665 | cgtcgtGACGAGGCGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGt > 1:706204/1‑71 (MQ=255)cgtcgtGACGAGGCGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGt > 1:378045/1‑71 (MQ=255)cgtcgtGACGAGGCGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGt > 1:430616/1‑71 (MQ=255)cgtcgtGACGAGGCGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGt > 1:181914/1‑71 (MQ=255)cgtcgtGACGAGGCGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGt > 1:1244409/1‑71 (MQ=255)cgtcgtGACGAGGCGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGt > 1:745916/1‑71 (MQ=255)cgtcgtGACGAGGCGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGt > 1:574993/1‑71 (MQ=255) gcgTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGTTGAACAGAGTAAAg < 1:81384/71‑1 (MQ=255) gcgTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGTTGAACAGAGTAAAg < 1:1521571/71‑1 (MQ=255) gcgTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGTTGAACAGAGTAAAg < 1:889485/71‑1 (MQ=255) gcgTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGTTGAACAGAGTAAAg < 1:160545/71‑1 (MQ=255) gcgTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGTTGAACAGAGTAAAg < 1:475941/71‑1 (MQ=255) gcgTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGTTGAACAGAGTAAAg < 1:1849521/71‑1 (MQ=255) gcgTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGTTGAACAGAGTAAAg < 1:1952940/71‑1 (MQ=255) gcgTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAGAAGAAAAATGACAAGGTTGAACAGAGTAAAg < 1:1896880/71‑1 (MQ=255) cgTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGTTGAACAGAGTAAAg < 1:629850/70‑1 (MQ=255) cgTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGTTGAACAGAGTAAAg < 1:686014/70‑1 (MQ=255) cgTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGTTGAACAGAGTAAAg < 1:732718/70‑1 (MQ=255) aGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACg < 1:844377/70‑1 (MQ=255) aGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACg < 1:1235627/70‑1 (MQ=255) gCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGTTGAACAGAGTAAAGCAGGAACCTTTAAAGAAGTc > 1:1303144/1‑71 (MQ=255) gCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTc > 1:1189881/1‑71 (MQ=255) gCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTc > 1:1189258/1‑71 (MQ=255) gCGTCGATCCGGGAGACAAGAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTAAAAGAAGTc > 1:916044/1‑71 (MQ=255) tCGATCCGGGAGACAAGAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTcgcg < 1:1776527/71‑1 (MQ=255) tCGATCCGGGAGACAAGAAGAAAAATGATAAGGTGGAACAGAGTAAAGCACGAACCTTTAAAGAAGTcgcg < 1:696409/71‑1 (MQ=255) cGATCCGGGAGACAAGAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTcgcg < 1:123014/70‑1 (MQ=255) ccGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGAtt > 1:257155/1‑69 (MQ=255) ccGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGAtt > 1:528440/1‑69 (MQ=255) ccGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGAtt > 1:510523/1‑69 (MQ=255) ccGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGAtt > 1:1053201/1‑69 (MQ=255) ccGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGAtt > 1:1099478/1‑69 (MQ=255) agaCAAGAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACg < 1:1578494/42‑1 (MQ=255) | CGTCGTGACGAGGCGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATT > W3110S.gb/4501555‑4501665 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |