| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,234,708 | 0 | T | A | 20.7% | 51.0 / 11.7 | 29 | intergenic (‑253/+31) | 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 T (11/12); new base A (3/3); total (14/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.33e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTAATCCAGGTGACGG > NC_000913/3234576‑3234853 | aGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTAc > 6:161423/1‑136 (MQ=255)aGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTAc < 5:161423/136‑1 (MQ=255) 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 < 7:22937/113‑1 (MQ=38) gCCGGATGCGGCGTAAACGCCTTATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGc > 8:22937/1‑113 (MQ=38) ggATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTAATCCa < 3:126179/149‑1 (MQ=255) aTAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGt < 7:132807/74‑1 (MQ=255) aTAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGt > 8:132807/1‑74 (MQ=255) gTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTa < 2:126191/135‑1 (MQ=255) gTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTa > 1:126191/1‑135 (MQ=255) tttACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTAATCCAGGTGACgg > 4:111832/1‑147 (MQ=255) tttACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTAATCCAGGTGACgg < 3:111832/147‑1 (MQ=255) | AGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAATCTTACCGACGACGTCCCAGCTTCACCGCTTTCAGTACCACAAATTTATTATTCGTCGCAATAGTGGTGCAGTTGCCGAAAATCTTCTTCAGTTTATGGAAGTAATCCAGGTGACGG > NC_000913/3234576‑3234853 |
| 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: ‑ |