| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | AM260480 | 153,689 | +G | intergenic (‑67/+93) | h16_B0133 ← / ← h16_B0134 | conserved hypothetical protein/outer membrane protein (porin) |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | AM260480 | 153,687 | 1 | . | G | 93.9% | 116.3 / 0.5 | 33 | intergenic (‑65/+95) | h16_B0133/h16_B0134 | conserved hypothetical protein/outer membrane protein (porin) |
| Reads supporting (aligned to +/- strand): ref base . (1/1); new base G (13/18); total (14/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.23e-01 | |||||||||||
GTCGGTGGTGTTGGCCGTGTTCTCGGGCGGGGCGTTGTGGATCGTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCA‑GG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGCCCGCTCGGCGCGGGCATGGCGTCTACGAAAAGTCAGAACTGGTGACGCAGGCCGGCCA > AM260480/153560‑153807 | gTCGGTGGTGTTGGCCGTGTTCTCGGGCGGGGCGTTGTGGATCGTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCA‑gggctt > 2:760796/1‑130 (MQ=255) tCGGTGGTGTTGGCCGTGTTCTCGGGCGGGGCGTTGTGGATCGTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCA‑gggcttc < 1:382297/134‑6 (MQ=255) gtggtgTTGGCCGTGTTCTCGGGCGGGGCGTTGTGGATCGTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCat > 1:611080/1‑133 (MQ=255) ggtgTTGGCCGTGTTCTCGGGCGGGGCGTTGTGGATCGTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGGACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCatgaa > 1:952386/1‑134 (MQ=255) gCCGTGTTCTCGGGCGGGGCGTTGTGGATCGTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAg < 1:832145/134‑1 (MQ=255) gtgGATCGTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCt < 1:239382/133‑1 (MQ=255) gTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCaa < 2:666206/134‑1 (MQ=255) tGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAg > 2:343836/1‑103 (MQ=255) tGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAg < 1:343836/103‑1 (MQ=255) tGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCa > 1:1159753/1‑132 (MQ=255) gTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCA‑GGACTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAg > 1:472822/1‑133 (MQ=255) cACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAGCGTctg < 2:978169/105‑1 (MQ=255) cACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAGCGTctg > 1:978169/1‑105 (MQ=255) ccACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAGCGTCTTCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGc < 2:194666/133‑1 (MQ=255) ccACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGc < 1:33893/133‑1 (MQ=255) cTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCAt > 1:927118/1‑107 (MQ=255) cTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCAt < 2:927118/107‑1 (MQ=255) cTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGcc < 2:1159753/134‑1 (MQ=255) cggcATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGCCCGCTCGGCGCggg < 1:800225/134‑1 (MQ=255) cgcgAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGCCCGCTCGGCGCGGGCATGGCGTCt < 2:952386/134‑1 (MQ=255) cAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGCCCGCTCGGCGCGGGCATGGCGTCTACGaa > 1:304991/1‑133 (MQ=255) aGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCatgaatg < 1:946034/54‑1 (MQ=255) aGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCatgaatg > 2:946034/1‑54 (MQ=255) acCGCACCGTTGGCGTGCGGGCTGTCCA‑GGACTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGCCCGCTCGGCGCGGGCATGGCGTCTACGAAAAGTCAGAACTGGt < 2:472822/134‑1 (MQ=255) cGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGCCCGCTCGGCGCGGGCAt < 2:185879/106‑1 (MQ=255) cGCACCGTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGCCCGCTCGGCGCGGGCAt > 1:185879/1‑106 (MQ=255) gTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGAt < 1:566711/74‑1 (MQ=255) gTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGAt > 2:566711/1‑74 (MQ=255) gTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGcc > 1:1038520/1‑84 (MQ=255) gTTGGCGTGCGGGCTGTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGcc < 2:1038520/84‑1 (MQ=255) tGCGGGCTGTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGCCCGCTCGGCGCGGGCATGGCGTCTACGAAAAGTCAGAACTGGTGACGCAGGCCGGCCa < 1:760796/134‑1 (MQ=255) ggCTGTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCa < 2:564866/45‑1 (MQ=38) ggCTGTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCa > 1:564866/1‑45 (MQ=38) gTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCAt < 2:928103/42‑1 (MQ=38) gTCCAGGG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCAt > 1:928103/1‑42 (MQ=38) | GTCGGTGGTGTTGGCCGTGTTCTCGGGCGGGGCGTTGTGGATCGTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCA‑GG‑CTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGCCCGCTCGGCGCGGGCATGGCGTCTACGAAAAGTCAGAACTGGTGACGCAGGCCGGCCA > AM260480/153560‑153807 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |