| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 1,107,556 | G→C | 12.1% | K236N (AAG→AAC) | hisC → | histidinol‑phosphate aminotransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 1,107,556 | 0 | G | C | 12.1% | 88.2 / 4.8 | 33 | K236N (AAG→AAC) | hisC | histidinol‑phosphate aminotransferase |
| Reads supporting (aligned to +/- strand): ref base G (14/15); new base C (2/2); total (16/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.00e-01 | |||||||||||
GGTGGTGAAACGGCGATCAGCCTGGTGGACCGCTACGATAACCTGCTGGTTACCCAGACGCTGTCGAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTCGGCGCGGCGGTGGCGTTCGAAGACCGTGAGTACTTCGAGGAAACCTGCCGCAAGGTAATAGACAGCCGCGAGGTGTTGGTCGGTCAGCTGCAGGCCAAGGGTTTT > NC_002947/1107410‑1107703 | ggtggtGAAACGGCGATCAGCCTGGTGGACCGCTACGATAACCTGCTGGTTACCCAGACGCTGTCGAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGaa < 8:41476/149‑1 (MQ=255) gtggtgAAACGGCGATCAGCCTGGTGGACCGCTACGATAACCTGCTGGTTACCCAGACGCTGTCGAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAAc < 6:279231/149‑1 (MQ=255) gtgAAACGGCGATCAGCCTGGTGGACCGCTACGATAACCTGCTGGTTACCCAGACGCTGTCGAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGc > 1:48822/1‑149 (MQ=255) aaaCGGCGATCAGCCTGGTGGACCGCTACGATAACCTGCTGGTTACCCAGACGCTGTCGAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCtt > 6:299851/1‑148 (MQ=255) ggCGATCAGCCTGGTGGACCGCTACGATAACCTGCTGGTTACCCAGACGCTGTCGAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACt < 5:182004/149‑1 (MQ=255) aGCCTGGTGGACCGCTACGATAACCTGCTGGTTACCCAGACGCTGTCGAAGTCGCGCTCGCTGGCGGGGCTACTAGTGGGCCTGGCGGTGTGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATc < 5:236133/148‑1 (MQ=255) cTACGATAACCTGCTGGTTACCCAGACGCTGTCTAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTGCAg < 5:299851/149‑1 (MQ=255) tACGATAACCTGCTGGTTACCCAGACGCTGTCGAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTg > 3:101524/1‑145 (MQ=255) tACGATAACCTGCTGGTTACCCAGACGCTGTCGAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTg < 4:101524/145‑1 (MQ=255) cGATAACCTGCTGGTTACCCAGACGCTGTCGAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAACAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTa > 8:84077/1‑149 (MQ=255) cGATAACCTGCTGGTTACCCAGACGCTGTCGAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAACAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTa > 8:79783/1‑149 (MQ=255) aaCCTGCTGGTTACCCAGACGCTGTCGAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTAtcgt > 2:55743/1‑149 (MQ=255) ctgctgGTTACCCAGACGCTGTCGAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCACCTCCTATCCACTGGACCGTGCAGCTATCGTcgg > 8:11045/1‑149 (MQ=255) ctgctgGTTACCCAGACGCTGTCGAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTcgg > 3:236407/1‑149 (MQ=255) ttACCCAGACGCTGTCGAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAACAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTcggcgcggcg < 7:84077/149‑1 (MQ=255) ttACCCAGACGCTGTCGAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAACAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTcggcgcggcg < 7:79783/149‑1 (MQ=255) ccAGACGCTGTCGAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTCGGCGCGGCGGTgg < 1:52761/149‑1 (MQ=255) aGACGCTGTCGAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTCGGCGCGGCGGTGGCg > 7:89137/1‑149 (MQ=255) aCGCTGTCGAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTCGGCGCGGCGGTGGCGt > 4:74272/1‑148 (MQ=255) tCGAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTCGGCGCGGCGGTGGCGTTCGAAg < 3:269720/148‑1 (MQ=255) aGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTCGGCGCGGCGGTGGCGTTCGAAGACCGt < 2:48822/149‑1 (MQ=255) tCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTCGGCGCGGCGGTGGCGTTCGAAGACCGt > 2:64413/1‑147 (MQ=255) gcgcTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTCGGCGCGGCGGTGGCGTTCGAAGACCGTGAGt < 1:64413/149‑1 (MQ=255) ctGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTCGGCGCGGCGGTGGCGTTCGAAGACCGTGAGTACTTCGa < 4:236407/149‑1 (MQ=255) cGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTCGGCGCGGCGGTGGCGTTCGAAGACCGTGAGTACTTCGAGGAAACCTGc < 3:74272/147‑1 (MQ=255) gAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTCGGCGCGGCGGTGGCGTTCGAAGACCGTGAGTACTTCGAgg > 3:259922/1‑138 (MQ=255) gAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTCGGCGCGGCGGTGGCGTTCGAAGACCGTGAGTACTTCGAgg < 4:259922/138‑1 (MQ=255) ggCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTCGGCGCGGCGGTGGCGTTCGAAGACCGTGAGTACTTCGAGGAAACCTGCCGCAAgg > 7:7662/1‑148 (MQ=255) cTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTc < 3:18677/57‑1 (MQ=255) cTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTc > 4:18677/1‑57 (MQ=255) aTCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTCGGCGCGGCGGTGGCGTTCGAAGACCGTGAGTACTTCGAGGAAACCTGCCGCAAGGTAATAGACAGCCGCGAg < 8:7662/149‑1 (MQ=255) gTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTCGGCGCGGCGGTGGCGTTCGAAGACCGTGAGTACTTCGAGGAAACCTGCCGCAAGGTAATAGACAGCCGCGAGGTGTTGGTCGGTCAGCTGc > 8:289471/1‑149 (MQ=255) aGAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTCGGCGCGGCGGTGGCGTTCGAAGACCGTGAGTACTTCGAGGAAACCTGCCGCAAGGTAATAGACAGCCGCGAGGTGTTGGTCGGTCAGCTGCAGGCCAAGGGtttt > 5:57813/1‑149 (MQ=255) gAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTCGGCGCGGCGGTGGCGTTCGAAGACCGTGAGTACTTCGAGGAAACCTGCCGCAAGGTAATAGACAGCCGCGAGGTGTTGGTCGGTCAGCTGCAGGCCAAGGGtttt > 5:84933/1‑148 (MQ=255) | GGTGGTGAAACGGCGATCAGCCTGGTGGACCGCTACGATAACCTGCTGGTTACCCAGACGCTGTCGAAGTCGCGCTCGCTGGCGGGGCTACGAGTGGGCCTGGCGGTGGGGCATCCGGACCTGATCGAGGCGTTGGAGCGGATCAAGAACAGCTTCAACTCCTATCCACTGGACCGTGCAGCTATCGTCGGCGCGGCGGTGGCGTTCGAAGACCGTGAGTACTTCGAGGAAACCTGCCGCAAGGTAATAGACAGCCGCGAGGTGTTGGTCGGTCAGCTGCAGGCCAAGGGTTTT > NC_002947/1107410‑1107703 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |