Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 355,050 | A→G | 13.8% | V42A (GTC→GCC) | cbcW ← | choline/betaine/carnitine ABC transporter permease |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 355,050 | 0 | A | G | 13.8% | 70.6 / 3.2 | 29 | V42A (GTC→GCC) | cbcW | choline/betaine/carnitine ABC transporter permease |
Reads supporting (aligned to +/- strand): ref base A (11/14); new base G (2/2); total (13/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGC > NC_002947/354910‑355184 | gCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGa > 1:529560/1‑150 (MQ=255) aGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGCAAGTCACGGCGCTGGTGACGCCATTGATCa < 2:402264/150‑1 (MQ=255) tAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATg < 2:533206/150‑1 (MQ=255) ccAGGTTGAGGGTCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGTCCACTTGCGCTGGCCGAAATGCGCAAGGGAGGCGAATAGGGCGCTAAGAAAGAAAGGGTTGAAAAAAGTCAGGGCGCGGGTGACGGCATGCATCATGAACt > 2:286636/1‑150 (MQ=255) ccAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATGAACt < 1:569853/150‑1 (MQ=255) ggATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGGGACGCCATGGATCATGAACTCCAGGCCtt > 1:642687/1‑150 (MQ=255) gcagGAACGACAGGGCGACAAATGCGGTCAGTGCCCCCTTGCGCTGGATGAAATTCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGa > 1:16688/1‑150 (MQ=255) aaCGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCg > 2:804656/1‑150 (MQ=255) aCGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGt < 2:362827/150‑1 (MQ=255) ggCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGaa > 1:518556/1‑147 (MQ=255) ggCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGaa < 2:518556/147‑1 (MQ=255) cGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGt < 2:553708/150‑1 (MQ=255) aGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTgg < 1:804656/150‑1 (MQ=255) aGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTgg < 2:365305/150‑1 (MQ=255) cACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGc > 2:42289/1‑150 (MQ=255) gcTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTcga < 1:286636/150‑1 (MQ=255) gATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGacgaa > 2:5620/1‑150 (MQ=255) aGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGGCGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTcgac > 1:483913/1‑134 (MQ=255) aGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGGCGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTcgac < 2:483913/134‑1 (MQ=255) gggCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACt < 2:305893/150‑1 (MQ=255) gcggcgAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTgg < 2:16688/150‑1 (MQ=255) ggcgAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGcc > 1:198000/1‑150 (MQ=255) gAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAg < 1:42289/150‑1 (MQ=255) tAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGGCGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCt > 1:720135/1‑150 (MQ=255) gaacgaacGGGTTGAAGAAAGTCAGGGCGCTGGTGGCGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGAt < 2:720135/150‑1 (MQ=255) aacgaacgGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGa > 1:670672/1‑148 (MQ=255) aacgaacgGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGa < 2:670672/148‑1 (MQ=255) gTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGc < 1:5620/150‑1 (MQ=255) gTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGc > 1:36759/1‑150 (MQ=255) | GCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAAGTCAGGGCGCTGGTGACGCCATGGATCATGAACTCCAGGCCTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGC > NC_002947/354910‑355184 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |