| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 1,511,862 | G→A | 6.8% | G26D (GGC→GAC) | rsmI → | 16S rRNA 2'‑O‑ribose C1402 methyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 1,511,862 | 0 | G | A | 6.8% | 168.0 / 2.9 | 60 | G26D (GGC→GAC) | rsmI | 16S rRNA 2'‑O‑ribose C1402 methyltransferase |
| Reads supporting (aligned to +/- strand): ref base G (26/29); new base A (2/2); total (28/32) | |||||||||||
| 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.50e-01 | |||||||||||
CGCCGGCGAAGTTCTTTGCGCCCTGTTGCAGGCGCTCGGGCTACAATGGCGCCTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTGGCATCGATACGCCGTTGGCCGCTTGTCATGAGCACAACGAGCGTG > NC_002947/1511716‑1512011 | cGCCGGCGAAGTTCTTTGCGCCCTGTTGCAGGCGCTCGGGCTACAATGGCGCCTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCaa < 1:477147/150‑1 (MQ=255) cGAAGTTCTTTGCGCCCTGTTGCAGGCGCTCGGGCTACGATGGCGCCTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTcga < 1:1493149/150‑1 (MQ=255) gAAGTTCTTTGCGCCCTGTTGCAGGCGCTCGGGCTACAATGGCGCCTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTcgac > 2:549445/1‑150 (MQ=255) gAAGTTCTTTGCGCCCTGTTGCAGGCGCTCAGGCTACAATGGCTCCTTTGATCCGAATGACAGGTGTGCGCGGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTcgac > 2:1032886/1‑150 (MQ=255) tttGCGCCCTGTTGCAGGCGCTCGGGCTACAATGGCGCCTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGa < 2:852740/150‑1 (MQ=255) tGCGCCCTGTTGCAGGCGCTCGGGCTACAATGGCGCCTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAgc < 1:1368096/150‑1 (MQ=255) ccTGTTGCAGGCGCTCGGGCTACAATGGCGCCTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTcga > 2:550058/1‑135 (MQ=255) ccTGTTGCAGGCGCTCGGGCTACAATGGCGCCTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTcga < 1:550058/135‑1 (MQ=255) ggCGCTCGGGCTACACTGGCGCCTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGaa < 2:1141408/150‑1 (MQ=255) gcgcTCGGGCTACAATGGCGCCTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCCACCTCGACGACATGAGCGCCCGGGCGCTGAAg < 1:580527/150‑1 (MQ=255) gcgcTCGGGCTACAATGGCGCCTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGCGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAg < 2:271009/150‑1 (MQ=255) tCGGGCTACAATGGCGCCTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCATGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGc > 1:420244/1‑150 (MQ=255) tCGGGCTACAATGGCGCCTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCATGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGc > 1:420263/1‑150 (MQ=255) gCTACAATGGCGCCTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTcgacg > 1:1616430/1‑119 (MQ=255) gCTACAATGGCGCCTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTcgacg < 2:1616430/119‑1 (MQ=255) gCTACAATGGCGCCTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGc < 2:675158/150‑1 (MQ=255) cAATGGCGCCTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGAc > 1:261298/1‑150 (MQ=255) aTGGCGCCTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGACAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGt > 1:1030081/1‑150 (MQ=255) gcCTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGgcgc < 2:261298/150‑1 (MQ=255) cTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTg < 2:1581237/150‑1 (MQ=255) tttGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGa < 2:1236947/150‑1 (MQ=255) tGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTATTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATc < 1:364223/150‑1 (MQ=255) tgctattGACATGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTTAAGGTGCTGGCGGACGTGGCGCTGATcgcc < 2:348614/144‑1 (MQ=255) tCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACa < 2:130474/103‑1 (MQ=255) tCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACa > 1:130474/1‑103 (MQ=255) gAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCgccgcc > 1:634191/1‑150 (MQ=255) tGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGgc > 2:771823/1‑90 (MQ=255) tGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGgc < 1:771823/90‑1 (MQ=255) gATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGCCGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGc < 1:397215/150‑1 (MQ=255) aTGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGcc > 1:1470150/1‑150 (MQ=255) aGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCa > 2:1407024/1‑150 (MQ=255) ttCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGa < 1:128620/93‑1 (MQ=255) ttCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGa > 2:128620/1‑93 (MQ=255) tCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACttt < 1:729092/150‑1 (MQ=255) cGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTg < 1:1407024/150‑1 (MQ=255) gAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGACAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTgg < 2:1030081/150‑1 (MQ=255) tgcacaGTGGGGACGCTGTATGTGGTGACAACCCCCATCGCCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTGGCATc < 2:290579/148‑1 (MQ=255) tCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGACAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTc > 1:1498983/1‑122 (MQ=255) tCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGACAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTc < 2:1498983/122‑1 (MQ=255) ccACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCAGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTGGCATCg > 2:1559382/1‑150 (MQ=255) gccGGTGCGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTGGCATCGa > 2:1593835/3‑150 (MQ=255) aCGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTGGCATCGAt > 1:1589518/1‑150 (MQ=255) gggACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTGGCATCGATACGCCg > 2:1432384/1‑150 (MQ=255) ggACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACttt < 1:938736/134‑1 (MQ=255) ggACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACttt > 2:938736/1‑134 (MQ=255) aTGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTGGCATCGATACGCCGTTGGCCGCtt < 1:549445/150‑1 (MQ=255) tgtgGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTTCTGCAACACTTTGGCATCGATACGCCGTTGGCCGCTTg > 2:715811/1‑150 (MQ=255) gtggCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAgg < 2:1067663/52‑1 (MQ=255) gtggCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAgg > 1:1067663/1‑52 (MQ=255) gtggCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTGGCATCGATACGCCGTTGGCCGCTTGTCAt > 1:519261/1‑150 (MQ=255) gtggCAACCCCCATCGGCAACCTAGAAGACATCAGACCCCGGGCACTGAACTTGCGTGCGGACTTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGACTGCTGAACCAGTTTGGCATCGATACGCCGTTGGCGGTTTGTCAt < 2:765647/150‑1 (MQ=255) gCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTGGCATCGATACGCCGTTGGCCGCTTGTCATGAg < 2:420244/150‑1 (MQ=255) gCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTGGCATCGATACGCCGTTGGCCGCTTGTCATGAg < 2:420263/150‑1 (MQ=255) aCCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCctgc > 2:1267313/1‑103 (MQ=255) aCCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCctgc < 1:1267313/103‑1 (MQ=255) cccATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTGGCATCGATACGCCGTTGGCCGCTTGTCATGAGCACAAc > 2:776753/1‑150 (MQ=255) ccATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTGGCATCGATACGCCGTTGGCCGCTTGTCATGAGCACAACg > 2:24262/1‑150 (MQ=255) caccGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTg > 2:323420/4‑110 (MQ=255) cATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTg < 1:323420/110‑1 (MQ=255) aTCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTGGCATCGATACGCCGTTGGCCGCTTGTCATGAGCACAACGAg > 2:1331460/1‑150 (MQ=255) gCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTg < 2:1224568/105‑1 (MQ=255) gCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTg > 1:1224568/1‑105 (MQ=255) gCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTGGCATCGATACGCCGTTGGCCGCTTGTCATGAGCACAACGAGCGTg > 1:447441/1‑150 (MQ=255) | CGCCGGCGAAGTTCTTTGCGCCCTGTTGCAGGCGCTCGGGCTACAATGGCGCCTTTGATCCGAATGACAGGTGTGCGCAGTGACTGATGTGGCAGGGGCTTCGAAATCCACGGTGGGGACGCTGTATGTGGTGGCAACCCCCATCGGCAACCTCGACGACATGAGCGCCCGGGCGCTGAAGGTGCTGGCGGACGTGGCGCTGATCGCCGCCGAAGACACCCGTCATTCGATCCGCCTGCTGCAACACTTTGGCATCGATACGCCGTTGGCCGCTTGTCATGAGCACAACGAGCGTG > NC_002947/1511716‑1512011 |
| 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: ‑ |