| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,413,735 | 0 | A | T | 50.0% | 23.7 / 23.5 | 33 | Y2* (TAT→TAA) | ydaC | DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein |
| Reads supporting (aligned to +/- strand): ref base A (2/14); new base T (16/0); total (19/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.09e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01 | |||||||||||
CAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTAT‑ACATGAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTC‑GAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTAC > NC_000913/1413588‑1413879 | cAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTtat‑a < 1:55581/150‑1 (MQ=255) tGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTT‑TAATGAAAACcccc > 1:115247/1‑150 (MQ=255) gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTAT‑ACATGAACACCTc < 2:587695/148‑1 (MQ=255) aaaCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTT‑AAATTAAAACCTCCCTTATTAtttt > 1:522700/1‑150 (MQ=255) gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTT‑AATTGTAAACCTCCCTtttttttt > 2:943458/1‑142 (MQ=255) gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTAT‑ACATGAACACCTCCATTATTATTTCCAGT‑gg < 1:519808/149‑1 (MQ=255) cttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TCATTTTTT‑TCATGAAAACCTCCCTCTTTAtttttagcggtgtg > 2:705981/1‑138 (MQ=255) tttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTAT‑ACATGAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTAt < 1:445641/150‑1 (MQ=255) ttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTAT‑ACATGAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTAtt < 2:826193/150‑1 (MQ=255) cTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGATAATTTT‑T‑AAAACAAAACCCCCACTATTATTTTCCGT‑GGGTGGtttttt > 1:903827/1‑146 (MQ=255) agagCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTAT‑ACATGAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCttt < 1:382588/151‑1 (MQ=255) agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTT‑ACATTAACAACCCCCTTATTATTTTCCGT‑GGGTCGTTTTTTTCCTCttt > 1:461561/1‑149 (MQ=255) aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTT‑ACATGAAAACCTCCCTTATTTTTTTCAAT‑GGGTCGTTTTTTTCCTCttt > 1:326355/1‑146 (MQ=255) gCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTAT‑ACATGAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTC‑GAGTGCTTCttt < 1:908994/151‑1 (MQ=255) tCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTT‑AAATGAAAACCCCCCTTTTTATTTTCAGTGGGTTGGTTTTTTCC‑TCTTTT‑GAGTGGtttttttt > 1:134273/1‑149 (MQ=255) aTTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑GAATTTTTT‑ATATTAACACCCCCCTTTTTTTTTTCTAT‑GGTTTGTTTTTTTTATCTTTT‑CTGTGCttttttttttttt > 2:1010437/1‑145 (MQ=255) aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTAT‑ACATGAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTC‑GAGTGCTTCTTTTTCACTTCCACCATa < 2:957302/150‑1 (MQ=255) ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTAT‑ACATGAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTC‑GAGTGCTTCTTTTTCACTTCCACCATAAcc < 1:191500/150‑1 (MQ=255) ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTAT‑ACATGAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTC‑GAGTGCTTCTTTTTCACTTCCACCATAAcc < 2:699047/150‑1 (MQ=255) ttttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAAATTTTT‑AAATAAAAACCCCCATTATTTTTTTCTGT‑GGGTTGTTTTTTTCCTCTTTC‑GTGGGCTTCTTTTTTACTTTCCCCCTAAcccggt > 1:576519/1‑145 (MQ=255) ttttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑CAATTTTTT‑AAAATAAAACCTCCCCTTTTTTTTTTATTGGGTTGGTTTTTTTC‑TTTTTT‑TAGTGCTTCTTTTTTATTTcttccctaaaacgtt > 2:601472/1‑135 (MQ=255) ttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTAT‑ACATGAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTC‑GAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCgg < 2:535318/150‑1 (MQ=255) atatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑GAATTTTTT‑TAATTAAAACCCCCCTTATTTTTTTCAGT‑GGGTTGTTTTTTTCATCTTTT‑TGGGGCTTTTTTTTTTCTTCCCCCATAAACGGGTCGGGGttt > 1:107282/1‑150 (MQ=255) atCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑GAATTTTTT‑TCATTAACACCTCCCTTTTTTTTTTTAGT‑GGGTTGTTTTTTTCCTTTTTT‑CTGTGCTTTTTTTTTTCTTTCCCCCTAAACCGTTCgggggt > 2:931252/1‑145 (MQ=255) cTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTT‑ACATGAAAACCCCCCTTTTTATTTTCCGT‑GGGTCGTTTTTTTCATCTTTC‑GTGTGCTTTTTTTTTTCTTTCACCATAAACGGTTCGGGGtt > 1:704194/1‑146 (MQ=255) ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTT‑AAATTAACACCCCCCTTATTTTTTTCCGT‑GGGTTGTTTTTTTCATCTTTT‑TTGTGGTTCTTTTTTTCTTCCCCCCCAAACCGGTCGGGGTTTGCATCCGtt > 1:918025/1‑151 (MQ=255) gTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAAATTTTT‑AAATGAAAACCTCCCTTTTTTTTTTCCGT‑GGGTTGTTTTTTTCATTCTTT‑GATTGGTTCTTTTTTTCTTCCCTCACAAAACGTTCGGGGTTTGCCTCCGtt > 2:598349/1‑148 (MQ=255) agtagGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTAT‑ACATGAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTC‑GAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACtt < 2:1040903/151‑1 (MQ=255) tagGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTCAT‑ACATGAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTC‑GAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTc < 2:767206/150‑1 (MQ=255) aGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTTTAAAAAGAAACCCCCCTTTTTTTTTTCCGGG‑GGTT‑GTTTTTTTCATTTTTCTTGGTGTTTTTTTTTTTCTTCCCCCTTAACGGGtcggggttttg > 2:34060/1‑122 (MQ=255) cACCTTCCTTTTCAATAGTGGCGG‑TAATTTTAT‑ACATGAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTC‑GAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAat < 1:383377/151‑1 (MQ=255) aCCTTCCTTTTCAATAGTGGCGG‑AAAATTTAT‑TCAAAAACACCTCCCTTTTTATTTCCTGT‑GGTTCGTTTTTTTTTTCTTTT‑GCGTGCTTTTTTTTTTTTTCCCCCATCACCCGTTCTGGTTTTCCACCCGTTACTCtttcctccgt > 2:114013/1‑139 (MQ=255) tCAATAGTGGCGG‑TAATTTTAT‑ACATGAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTC‑GAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAAt < 1:585034/150‑1 (MQ=255) aaTTTTAT‑ACATGAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTC‑GAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTAc > 1:1039977/1‑151 (MQ=255) | CAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAATTTTAT‑ACATGAACACCTCCATTATTATTTCCAGT‑GGTTCGTTTATTCCATCTTTC‑GAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTAC > NC_000913/1413588‑1413879 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |