| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,622,816 | 0 | G | A | 24.2% | 63.1 / 13.3 | 33 | intergenic (+75/+1) | mgtS/mgrR | small protein MgtS/small regulatory RNA MgrR |
| Reads supporting (aligned to +/- strand): ref base G (24/1); new base A (0/8); total (24/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.48e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTA > NC_000913/1622684‑1622951 | aataatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCGaaaaaaa > 2:416544/1‑140 (MQ=255) tttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAActggc < 2:201287/140‑5 (MQ=255) ttcggcccGCTTTTTTCAGCCAAAAGGGGGATAACTAATGAACGGAAATACCCCCCCACCTACCCGGCCCT‑TTGTTACGGTTGTGAAAAGGGGG‑CT‑TATTTTTTAGGAGGGCAAAAAAAAACCGCCAGTAAACCGGCGGt < 2:123310/134‑1 (MQ=255) gCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATg > 1:516058/1‑140 (MQ=255) cgcgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCt > 1:1086133/1‑140 (MQ=255) cgcgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCt > 1:24737/1‑140 (MQ=255) gcgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCtt > 2:205470/1‑140 (MQ=255) gcgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAA‑CGGCGGTGAATGCTTg > 1:602294/1‑140 (MQ=255) ccAGCCAAAAGGGGGATAAATAAGGAAGGGAGA‑AACCCCCCCCCCAACCCGCCCCTTTTTTACGTTTGTGAAAAGGGGG‑CT‑TTTTTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGa < 2:1245012/139‑1 (MQ=255) acccAAAAAGGGGATGAAAAATGAAGGGAAATACCCCCCCCCCAACCCGCCCCT‑TTTTTACGTTTGTGAAAAGGGGG‑CT‑TTTTTTTTAGGCGGGAAAAAAAAACCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATa < 2:41588/138‑1 (MQ=255) aaaaTGGGATAAAAAATAAACGAAAATACCCCCCCCCCAACCCGCCCCT‑TTTTTACATTTGTAAAAAGGGGG‑CT‑TTTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAttt < 2:24737/139‑1 (MQ=255) aaTGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTtgt > 1:1245012/1‑140 (MQ=255) aTAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGc > 1:1156988/1‑140 (MQ=255) aaaccccTCCCCTACCCGCCCCT‑TTTTTACGGTTTTGAAAAGGGGG‑CT‑TTTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCa < 1:71959/136‑1 (MQ=255) tCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCt > 1:897142/1‑140 (MQ=255) tCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCt > 1:363426/1‑140 (MQ=255) tCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCt > 1:1237980/1‑140 (MQ=255) ccccaacccGCCCCT‑TTGTTACAGTTGTGAACAGGGGG‑CT‑TGTTTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTg < 2:897142/133‑1 (MQ=255) cACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTg > 1:1362980/1‑140 (MQ=255) cACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTg > 2:205414/1‑140 (MQ=255) ctccccGCCCCT‑TTGTTACTGTTGT‑TAAAAGGGGTCT‑TTTTTTTTAGGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCAc < 1:205470/136‑1 (MQ=255) aCCGGCCCC‑TTGTTACAGTTGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGa > 1:343250/1‑140 (MQ=255) aCCGGCCCC‑TTGTTACAGTTGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGa > 2:1254217/1‑140 (MQ=255) cccc‑TTGTTACAGTTGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCG‑AAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGa > 1:1373103/1‑140 (MQ=255) ttGTTACAGTTGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATc > 1:1312507/1‑140 (MQ=255) ttGTTAAAGTTGTGAAAAGGGGG‑CT‑TAATTTTTAGGAGGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATc < 2:516058/140‑1 (MQ=255) tttaggTTTTGTAAAAGGGG‑CCTTTTTTTTTATGGCGGAAAAAAAAACCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGt < 2:1237980/135‑1 (MQ=255) aagtttttgtACAGGGGG‑CT‑TTTTTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGa < 2:406410/133‑1 (MQ=255) ttGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGAcaca > 1:406410/1‑140 (MQ=255) gtgtACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGt > 1:665489/1‑140 (MQ=255) gtgtACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGt > 1:1086465/1‑140 (MQ=255) gggg‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAg > 2:1176377/1‑140 (MQ=255) g‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGttt > 1:1328917/1‑140 (MQ=255) gttttttttttATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGttt < 2:1312507/134‑1 (MQ=255) aTTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa > 1:1251769/1‑140 (MQ=255) ttttATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaat < 2:1086465/140‑1 (MQ=255) cGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTa > 1:477560/1‑140 (MQ=255) | AATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGG‑CC‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTA > NC_000913/1622684‑1622951 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |