| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,234,716 | 0 | C | T | 20.0% | 59.7 / 10.1 | 30 | intergenic (‑261/+23) | higB/rlmG | mRNA interferase toxin of the HigB‑HigA toxin‑antitoxin system/23S rRNA m(2)G1835 methyltransferase, SAM‑dependent |
| Reads supporting (aligned to +/- strand): ref base C (11/13); new base T (3/3); total (14/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 = 3.44e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTAATCCAGGTGACGGTTGGC > NC_000913/3234580‑3234858 | gATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGt > 5:147421/1‑148 (MQ=255) ttttGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTAcc < 1:167781/148‑1 (MQ=255) ttttGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTAcc > 2:39693/1‑148 (MQ=255) tGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACcgac < 1:39693/148‑1 (MQ=255) ttttATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTAAACGCCTTATCCGGc > 1:137569/1‑123 (MQ=255) ttttATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTAAACGCCTTATCCGGc < 2:137569/123‑1 (MQ=255) tGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATGCGGCGTAAACGCCTTATCCg > 6:103994/1‑116 (MQ=38) tGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATGCGGCGTAAACGCCTTATCCg < 5:103994/116‑1 (MQ=38) aTTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCa < 8:29412/136‑1 (MQ=255) aTTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCa > 7:29412/1‑136 (MQ=255) cGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCa < 1:136781/78‑1 (MQ=255) cGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCa > 2:136781/1‑78 (MQ=255) ccGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTAcc > 8:84410/1‑78 (MQ=255) ccGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTAcc < 7:84410/78‑1 (MQ=255) cTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCCTTCAGTACCACAAATTTATTATTCGTCGCAATAgtg < 2:67/119‑1 (MQ=255) cTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCCTTCAGTACCACAAATTTATTATTCGTCGCAATAgtg > 1:67/1‑119 (MQ=255) tGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTtattat > 5:7131/1‑104 (MQ=255) tGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTtattat < 6:7131/104‑1 (MQ=255) aGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAATTTTATTATtcgt < 2:53452/101‑1 (MQ=255) aGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAATTTTATTATtcgt > 1:53452/1‑101 (MQ=255) gCCGGATGCGGCGTAAACGCCTTATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGc > 8:22937/1‑113 (MQ=38) gCCGGATGCGGCGTAAACGCCTTATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGc < 7:22937/113‑1 (MQ=38) ggATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTAATCCa < 3:126179/149‑1 (MQ=255) aTAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGt > 8:132807/1‑74 (MQ=255) aTAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGt < 7:132807/74‑1 (MQ=255) gTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTa > 1:126191/1‑135 (MQ=255) gTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTa < 2:126191/135‑1 (MQ=255) tttACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTAATCCAGGTGACgg > 4:111832/1‑147 (MQ=255) tttACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTAATCCAGGTGACgg < 3:111832/147‑1 (MQ=255) tACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTAATCCAGGTGACGGTTg < 6:147421/148‑1 (MQ=255) cgccgcATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTAATCCAGGTGACGGTTGGc > 1:90293/1‑148 (MQ=255) gccgcATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATcttctt > 4:11382/1‑114 (MQ=255) gccgcATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATcttctt < 3:11382/114‑1 (MQ=255) ccgcATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGaa > 2:137796/1‑125 (MQ=255) ccgcATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGaa < 1:137796/125‑1 (MQ=255) | GATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTAATCCAGGTGACGGTTGGC > NC_000913/3234580‑3234858 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |