| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 4,165,663 | A→G | 12.1% | R161R (CGT→CGC) | PP_3666 ← | MHS family metabolite MFS transporter |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 4,165,663 | 0 | A | G | 12.1% | 77.3 / 3.6 | 33 | R161R (CGT→CGC) | PP_3666 | MHS family metabolite MFS transporter |
| Reads supporting (aligned to +/- strand): ref base A (16/13); new base G (2/2); total (18/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01 | |||||||||||
GAACGGCGCCGCCAGCAGGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGC > NC_002947/4165516‑4165805 | gAACGGCGCCGCCAGCAGGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCAcg > 1:687076/1‑150 (MQ=255) aCGGCGCCGCCAGCAGGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCGcgcc > 2:398652/1‑150 (MQ=255) cagGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGgcgc < 2:208066/150‑1 (MQ=255) agGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCGCGCCGCCTGGgcggcgtg < 1:398652/150‑3 (MQ=255) ggCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCa > 2:515801/1‑109 (MQ=255) ggCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCa < 1:515801/109‑1 (MQ=255) ggCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCGCGGCGTa > 1:281567/1‑150 (MQ=255) ttCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGc > 1:280477/1‑132 (MQ=255) ttCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGc < 2:280477/132‑1 (MQ=255) ttCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTc > 2:340038/1‑150 (MQ=255) tCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGTCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAAc > 1:498233/1‑150 (MQ=255) cATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCGCGCCGCCTGGgcggcgtgg > 1:247726/1‑96 (MQ=255) cATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCGCGCCGCCTGGgcggcgtgg < 2:247726/99‑4 (MQ=255) aTGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAAccc < 1:168241/150‑1 (MQ=255) tGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCACCGGCGTACTCACCCGAGGCGGAAAAcccc > 2:363066/1‑150 (MQ=255) ccgcACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATcagca > 2:139848/1‑150 (MQ=255) aCCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATcagcaacag > 2:135830/1‑150 (MQ=255) ccGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATcagcaacagc > 2:136058/1‑150 (MQ=255) cGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATcagcaacagca > 2:30814/1‑150 (MQ=255) caacaacCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGc < 1:30814/150‑1 (MQ=255) acaacCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCg > 2:81722/1‑150 (MQ=255) caacCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCgg < 2:281567/150‑1 (MQ=255) cccGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCg > 1:441355/1‑100 (MQ=255) cccGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCg < 2:441355/100‑1 (MQ=255) gccgTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGc < 2:498233/150‑1 (MQ=255) atcgTGCCGGCCAGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGc > 2:777492/3‑150 (MQ=255) cTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATc < 2:687076/150‑1 (MQ=255) gCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTgc < 1:340038/150‑1 (MQ=255) gCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACGCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTgc < 1:571095/150‑1 (MQ=255) aCGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAAc > 1:504210/1‑90 (MQ=255) aCGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAAc < 2:504210/90‑1 (MQ=255) tACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAg > 2:464306/1‑150 (MQ=255) aCAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGTTCATCGGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGc < 2:707005/150‑1 (MQ=255) | GAACGGCGCCGCCAGCAGGAACGGCAGGCGCCAGCCCCACTCATGCAGGGCTTCACTGCTGAGCAGCTCATGCAGCACCGCGACGAACGCCGCACCGAACAACAACCCCGCCGCCGTGCTGGCCGGCACGATGCTGGTGTACAAGCCACGCCGCCCAGGCGGCGCGTACTCGGCCAGGAACGCAGCGGCCCCGGCGTACTCACCCGAGGCGGAAAACCCCTGTACCAGGCGGATCAGCAACAGCAAGGCCGGTGCCCACAGCCCGATCTGCGCATAGCCCGGCAACAAGC > NC_002947/4165516‑4165805 |
| 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: ‑ |