| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,171,892 | 0 | C | G | 42.1% | 45.8 / 19.3 | 38 | intergenic (+136/‑165) | rrfB/murB | 5S ribosomal RNA of rrnB operon/UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding |
| Reads supporting (aligned to +/- strand): ref base C (16/6); new base G (0/16); total (16/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.62e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAC‑GTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCG > NC_000913/4171777‑4172034 | gggggaaggcgggcccTTTCTTTTTTTGTTTTGTTGGGGGGGGACCGCTCCCGTGAGTAGGAAAAACCCGCGGGGGGGGGATTTGAAC‑GTTGCAAAGCAACGCCCG‑GGGGGGGGGGGGGCAGG‑CCGCCCGCCATAAACTGCCAGGCAt < 1:1106806/134‑1 (MQ=2) gCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAC‑GTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTAc > 1:528141/1‑151 (MQ=35) ctctCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAC‑GTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTGAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTAc > 2:627366/1‑150 (MQ=25) tctcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAC‑GTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACa > 2:1170538/1‑150 (MQ=35) tctcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAC‑GTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACa > 1:844177/1‑150 (MQ=35) ctcCGGGGTAGGACAAACCCGCGGGGGGGGGTTTTAAAC‑TTTGGAAAGCACCGGCCG‑GGAGGGTGGGGGGCGGG‑ACCCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACa < 2:1532407/149‑1 (MQ=11) cccgtGAGTAGAAAAATTCCGCGGGGAGGGGTTTTGACC‑GTGGGGAAGCACCGCCCG‑GGGGGGGGGGGGGCGGG‑CCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTa < 1:1349218/143‑1 (MQ=11) ccggAGTAGGAAAAACCCGCGGGGAGCGGATTGGAAG‑GTGGAGAACAAGCGGCCC‑GGGGGGGGGGGGGCGGG‑ACCCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACa < 1:1572916/144‑1 (MQ=11) ccggAGCAGTACAACCCCGCGGGGAGCGGTTTTGAAC‑GTTGCGACGCATCGGCCC‑GGAGGGGGGCGGGAAGG‑CCGTCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACaaa < 2:1491192/146‑1 (MQ=11) ccTGAGTAGGCAAACCCCGCCGGGAGC‑GATTTAAACTGTGGCGAAACAACGCCCC‑GGGGGGGGGGGGGCAGA‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAAc < 2:642538/150‑1 (MQ=14) ccTGAGTAGCAAAACCCCGGGGGGGGCGGTTTTAAGC‑GTGTCAAAGAAGCGGCCG‑GGGGGGTGGGGGGGGGG‑CCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAAc < 1:1294034/150‑1 (MQ=9) cTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAC‑GTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAActc > 2:281710/1‑151 (MQ=39) gaacaaccccgcggGGGGGGGTTTTGAAC‑TTGGCGAAGAAGCGGCCG‑GGGGGGGGGGGGGCAGG‑CCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTg < 1:1460766/142‑1 (MQ=14) aaccccgcggGGGGGGAATTTGAAG‑GTGGCAAACCAACGGCCG‑GGGGGGGGGCGGGAAGG‑CCCCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTgtc < 1:304978/144‑1 (MQ=14) aaccccgcggGGGGCGTATTTGAAG‑GTGGCGAAGCAGCGGCCG‑GGGGGGGGGGGGGCGGG‑CCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTgtcg < 1:657791/145‑1 (MQ=17) aaaCCCGCGGGGAGCGGTTTTGAGC‑TTTGAAAAACAGGGGCGG‑GGGGGGGGGGGGGCAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTgt < 1:1113290/147‑1 (MQ=14) aaTCCGCCGGGAGCGGATTTGAAC‑GTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTgtcgtc > 2:54217/1‑150 (MQ=255) ccgccgGGAGCGGATTTGAAC‑GTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCatat > 2:1649949/1‑151 (MQ=255) gccgGGAGCGGATTTGAAC‑GTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATc > 2:217024/1‑150 (MQ=255) gggagggtttttGAAC‑GTTGCGAAGCAACGGCCC‑GGAGGGGGGGGGGCAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACa < 2:531464/142‑1 (MQ=255) gggAGCGGATTTGAAC‑GTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTAc > 2:536658/1‑150 (MQ=255) gggAGCGGATTTGAAC‑GTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTAc > 2:546037/1‑150 (MQ=255) gggAGCGGATTTGAAC‑GTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTAc > 1:1117492/1‑150 (MQ=255) ggAGCGGATTTGAAC‑GTTGCGAAGCAACGGCCG‑GGGGGGTGGCGGGCAGG‑CCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACaa < 2:309168/151‑1 (MQ=255) gggCGGATTTGACC‑GTTGCGAAGCACCGGCGG‑GGGGGGGGGGGGGAAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAg < 2:1107329/149‑1 (MQ=18) ggATTTGAAC‑GTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCAt > 2:1096127/1‑151 (MQ=255) gattttGAAC‑TTTGCGAAGCACCGGCCG‑GGAGGGTGGGGGGCAGG‑CCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCa < 2:821377/147‑1 (MQ=18) gaaTTTGAAC‑GTTGC‑AACCACAGGCCCGGGGGGGGGGCGGGCGGG‑CCCCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCAt < 1:180641/149‑1 (MQ=18) tttGAAT‑GTTGGGAAGCAACGCCCG‑GGGGGGGGGCGGGCAGG‑CCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATcc < 2:1306533/150‑1 (MQ=18) tttGAAC‑GTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATccc > 1:69317/1‑151 (MQ=255) ttGACC‑GTTGCAAAGCAACGGCCG‑GGAGGGTGGGGGGCAGG‑CCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATccc < 1:1541648/150‑1 (MQ=255) ttGAAC‑GTTGGAAACCACCGGCCG‑GGGGGGTGGCGGGAAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATccc < 2:620836/150‑1 (MQ=18) gTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCcac > 2:183295/1‑149 (MQ=255) ttgaaaaGCACCGGCCC‑GG‑GGGGGGGGGGAAGGTACCCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCca < 2:374046/142‑1 (MQ=21) ttGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAAATAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCcac > 1:1158904/1‑148 (MQ=255) gAAGCAACGGCCC‑GGAGGGTGGCGGGCAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACg > 2:371250/1‑151 (MQ=255) ggCCG‑GGGGGGGGGGGGGCGGG‑CCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGt < 1:1080791/145‑1 (MQ=25) gccggggggggTGGGGGGCAGG‑GCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTaaa < 2:1255995/140‑1 (MQ=255) gggggggggggggCAGG‑CCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCt < 2:593079/147‑1 (MQ=255) ggggggggggggCAGG‑CCCCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCg < 1:439247/149‑1 (MQ=255) | GTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAC‑GTTGCGAAGCAACGGCCC‑GGAGGGTGGCGGGCAGG‑ACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAACTAGCCTCG > NC_000913/4171777‑4172034 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |