Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,429,990 | 0 | A | T | 35.0% | 64.1 / 17.0 | 40 | M43L (ATG→TTG) | cycA | D‑alanine/D‑serine/glycine transporter |
Reads supporting (aligned to +/- strand): ref base A (8/18); new base T (14/0); total (22/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.20e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGAGTCCCGATACAGATTCGTCGTATCATAGACTGACTAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGCCATTGG‑TACGGGGTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑TTTTCGTTTATATGATCA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTG‑CTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAG > NC_000913/4429776‑4430215 | aTGAGTCCCGATACAGATTCGTCGTATCATAGACTGACTAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGGGGCCATTGG‑ATGGGGGTGGTTTTT‑GGGGTC‑GGCCAAACGaatta > 2:1470270/1‑233 (MQ=255) cGATACAGATTCGTCGTATCATAGACTGACTAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGCCATTGG‑TACGGGGTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCggg < 1:995924/239‑1 (MQ=255) gATTCGTCGTATCATAGACTGACTAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGCCCCTGG‑CACGGGGGTTTTTTT‑GGGGTTTGGCACAACAC‑TTTGCCCTTCCCGGGcctcgt > 1:1114331/1‑235 (MQ=255) cgtATCATAGACTGACTAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGCCATTGG‑TACGGGGTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑tttt < 1:1001280/239‑1 (MQ=255) actgactAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGTCATTGG‑AACCGGGGTGTTTTT‑GGGGGCTGGCGAAAAAA‑GTAG‑GCATTTCGGGGGGTCGAGAC‑TTTTCGTTTATATg > 2:1001280/1‑239 (MQ=255) actAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGCCCTTTG‑TACGGGGGTGTTTTT‑GGGGTCTGGCAAAACGA‑TTAG‑CCCTTCCGGGCCGTCGGTCA‑TTTTTGTTTTTTTGTTCa > 1:696384/1‑239 (MQ=255) ctAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGCCATTGG‑TACGGGGTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑TTTTCGTTtata < 1:1470270/232‑1 (MQ=255) aGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGCCATTGG‑TACGGGGGTGTTTAT‑TGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGGCGATCA‑TTTTTGTTTATATGATCA‑TTGGt > 1:50284/1‑239 (MQ=255) gTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGCCCTTGG‑GACGGGGGTGTTTTT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGTCGGGCCGGCGGTCC‑TTTTCGGTTATATGATAA‑TTTGGTTTAt > 1:211165/1‑239 (MQ=255) gTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGCCATTGG‑TACGGGGTTGTTTAT‑GGGGGCTGGCAAAAAGA‑TTAG‑CCTTGCCGGGCGCGGGATCG‑TTTTTGTTTTTATGGTCATTTGGttttt > 1:1252155/1‑238 (MQ=255) agagCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGCCCTTGG‑GAAGGGGGTGTTTTT‑TGGGGGTGGCAAAACGA‑ATAG‑CCTTTCCGGGCCGTCGAGCA‑TTTTCGTTTTTATGATTC‑TTGGGTTTAt > 1:509647/1‑237 (MQ=255) cTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGCCATTGG‑TACGGGGTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑TTTTCGTTTATATGATCA‑TTGGTTTTATGctct < 1:1485676/237‑1 (MQ=255) gTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGCCATTGG‑TACGGGGTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑TTTTCGTTTATATGATCA‑TTg < 1:361027/208‑1 (MQ=255) tACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGGCCTTGG‑GACGGGGGTGTTTTT‑GGGGTCTGGCAAAACGA‑ATAG‑CCTTGCCGGGCCGCCGGGCC‑TTTTCGTTTTTATTATAA‑TTGGTTTTATTCTCTTTTTCGTGGTTCGGGCAAta > 2:554659/1‑238 (MQ=255) tACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGCCATTGG‑TACGGGGTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑TTTTCGTTTATATGATCA‑TTGGTTTTATGCTCTTTTTCGTGATGCggg < 2:518412/234‑1 (MQ=255) tACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGCCATTGG‑TACGGGGTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑TTTTCGTTTATATGATCA‑TTGGTTTTATGCTCTTTTTCGTGATGCg < 2:973123/232‑1 (MQ=255) ggTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGGGGCATTGGCTGCGGGG‑TGTTTTA‑GGGGCGTGGCCAAAAGA‑TCAG‑CCCTGCCGCGGCGCCGTTCA‑TTTTCCGTTTTTTGATCA‑TTCGTTTTTTGTTCTTTTTTGTTGTGCGGGCAATGGGGGGATTGCTG‑CTTTCGa > 2:896646/1‑239 (MQ=255) atgatCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCCGCTTATTGCCATTGGC‑GGTGGCCTTCC‑GCTGGGGGTGTTTTTGGGGGTTTGGC‑AAACGATTTCG‑CCTTCCCGGG‑CGTCCGTCATTTTTCGTTTTATTGTTCT‑TTGGTTTTTCGCT‑TTTTTTGTGTTTGGGGCCAGTGGGGGATTGGTGCCTTTCGTACCTGGATTA‑ACATCtttttatgga > 2:1920/1‑233 (MQ=255) atCAGGCTCCGGCGGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGCCATTGG‑TACGGGGTTGTTTAT‑GGGGTCTGGCaa < 1:1323985/109‑1 (MQ=255) cTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGCCATTGG‑TACGGGGTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑TTTTCGTTTATATGATCA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTG‑CTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCt < 1:554659/236‑1 (MQ=255) ggCTGAACAGTCGCTACGGCGCAATCTCGCAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGCCATTGG‑TACGGGGTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑TTTTCGTTTATATGAt < 1:322228/144‑1 (MQ=255) ggCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGCCCTTGG‑TACGGGGGTGTTTTT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑TTTTCGTTTATATGATCA‑TTGGTTTTATGCTCTTTTTTGTGATGCGGGGAATGGGGGGATTGCTG‑CTTTTGAATCTGGAATACAAATCTTTTTGTGGCTTCggcttc > 2:812618/1‑230 (MQ=255) ggCCCAATCTCACAAACCGAAATATTCAGCTTATTGCCATTGGC‑GGTGCCATTGG‑TAGGGGGTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGCTCA‑TTTTCGTTTATATGATCA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTG‑CTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACGCGGGCCGTGGGCa < 2:1322159/239‑1 (MQ=255) aCATATTCAGCTTATTGCCATTGGC‑GGGGCCATTTG‑TACGGGGGTGTTTTT‑GGGGTCTGGCCAAACGA‑ATAG‑ACCTGCCGGGCCGGCGATCA‑TTTTTGTTTTTATGATTA‑TTGGTTTTTTTCTCTTTTTTGTGATGCGGGCAATGGGGGGATTTGTG‑CTTTTGTATCTGGTATACAAAACTTTTATTGACTTCGCTTCCGATTTACTCCGGGCGTGGGCAGGGTATTTTACAGGccgt > 2:1559828/1‑235 (MQ=255) atTCAGCTTATTGCCATTGGC‑GGTGCCATTGG‑TACGGGGTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑tt < 1:714911/90‑1 (MQ=255) cTTATTGCCATTGGC‑GGTGCCATTGG‑TACGGGGGTGTTTTT‑GGGGGCTGGCCAAAAGA‑ATAG‑GCTTGCCGGGGCGTCGGTCA‑TTTTTGTTTATATGATCA‑TTTGTTTTTTTCTCTTTTTTTTGATGCGGGGAATGGGGGGAATGTTG‑CTTTCTGATCTTGAATACAAAACTTTTAGTGAACTCGCTTCCGGTTTACTCCGGGCGGGGGGAGGGTATTTTAACGGCCGGGCTTACttggt > 2:1602997/1‑235 (MQ=255) ttATTGCCATTGGC‑GGTGCCATTGG‑TACGGGGTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑TTTTCGTTTATATGATCA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTG‑CTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTgg < 1:1602997/236‑1 (MQ=255) ttGCCATTGGC‑GGTGCCATTGG‑TACGGGGGTGTTTTT‑TGGGGCTGGCAAAAAGG‑TTAG‑CCCTGCCGGGGCGGCGATCA‑TTTTCTTTtata > 2:714911/1‑90 (MQ=255) tGCCATTGGC‑GGTGCCATTGG‑TACGGGGGTGTTTAT‑g > 1:190167/1‑37 (MQ=255) ccATTGGC‑GGTGCCATTGG‑TACGGGGTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑t < 1:1550374/76‑1 (MQ=255) ccATTGGC‑GGTGCCATTGG‑TACGGGGTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑TTTTCGTTTATATGATCA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTG‑CTTTCGAATCTg < 1:1270907/152‑1 (MQ=255) tGCCGTTGG‑TACGGGGGTGTTTAT‑GGGGGCGGGCAAAACGA‑ATAG‑CCTTGCCGGGCCGGCGGTCT‑TTTTTGTTTATATGGTCA‑TTGGGTTTATGCTCTTTTTTGTGATGCCGGCAATGGGGGAATTGCTG‑CTTTCGAATCTTGGATAAAAAACTTTTTGTTACTTCGCTTCCGGTTTTTTCGGGCCGTGGGGAGGAGATTTTTCCGGCGGGGCTTTCTTGGTCTGCTGGGgtgttaacc > 2:382662/1‑234 (MQ=255) tGCCATTTG‑TACGGGGGTGTTTTT‑GGGGGCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑TTTTTGTTTATATGATCA‑TTGGTTTTTTGGTCTTTTTCGGGATGCGGGCAATGGGGGGATTGCTG‑CTTTTGAATCTTGGATACAAAACTTTTTGTGACTTCGCTTCCGGTTTTCTCGGGCCGTGGGGAGGATATTTTACCCGCTGGGCTTTATGGGTCTCGTGGGGTGTAACCg > 2:113016/1‑239 (MQ=255) aCGGGGTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑TTTTCGTTTATATGATCA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTG‑CTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTa < 2:1195629/225‑1 (MQ=255) aCGGGGTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑TTTTCGTTTATATGATCA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTG‑CTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTa > 1:1195629/1‑225 (MQ=255) aCGGGGTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑TTTTCGTTTATATGATCA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTG‑CTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGCAACCGGTATGGCAg > 2:1328423/1‑238 (MQ=255) ggTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑TTTTCGTTTATATGATCA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTG‑CTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTg < 2:1088669/196‑1 (MQ=255) ggTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑TTTTCGTTTATATGATCA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTG‑CTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTg > 1:1088669/1‑196 (MQ=255) ggTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑TTTTCGTTTATATGATCA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTG‑CTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCgg > 1:481015/1‑226 (MQ=255) ggTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑TTTTCGTTTATATGATCA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTG‑CTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCgg < 2:481015/226‑1 (MQ=255) | ATGAGTCCCGATACAGATTCGTCGTATCATAGACTGACTAAAGGCCGTAGAGCCTGAACAACACAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGC‑GGTGCCATTGG‑TACGGGGTTGTTTAT‑GGGGTCTGGCAAAACGA‑TTAG‑CCTTGCCGGGCCGTCGATCA‑TTTTCGTTTATATGATCA‑TTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTG‑CTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAG > NC_000913/4429776‑4430215 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |