Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,624,517 | 0 | G | T | 49.1% | 63.4 / 36.1 | 59 | E492* (GAA→TAA) | ppk | polyphosphate kinase, component of RNA degradosome |
Reads supporting (aligned to +/- strand): ref base G (8/22); new base T (29/0); total (37/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.27e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAG‑ATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGTC > NC_000913/2624371‑2624666 | gTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTAATTAAA‑TTTTTTTaaa > 2:342925/1‑150 (MQ=255) aaaGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTAA‑TTTATTTaaaaa > 1:57625/1‑149 (MQ=255) aaGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTTAA‑TTTTTTTaaaaaac > 2:116376/1‑148 (MQ=255) aGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAA‑TTTATTTAAAAacc > 1:968007/1‑150 (MQ=255) gTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTc‑c < 1:762417/151‑1 (MQ=255) aaGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTAC‑TTTTTTTTAAAACCATACCGTC‑Cggg > 1:971745/1‑150 (MQ=255) gtggtgCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGTTTTTAA‑TTTTTTTAAAACCCCAACCCCC‑CGGTg > 2:865549/1‑150 (MQ=255) aggtgCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGAGTTTAA‑TTTATTTAAAAacccaaccaccc > 2:86397/2‑136 (MQ=255) gtgCGATATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGTCTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGAGCTTTTAG‑GTTTTTTAAAAACCATCCCGTc‑ctggggc > 1:806627/1‑143 (MQ=255) tgCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAAC‑TTTTTTTAAAAACCCTACCGCC‑CGGGGGCAtt > 2:750307/1‑151 (MQ=255) aTGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTTAC‑TTTATTGAAAAACCATACCCTC‑CGGGGGCATTTGatt > 1:63645/1‑150 (MQ=255) acaTCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTAC‑TTTTTTTAAAAACCACACCGTC‑CGGGGGCATTTTTTTAttt > 2:132887/1‑148 (MQ=255) cGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTa < 1:1032119/151‑1 (MQ=255) cGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTa < 2:826671/151‑1 (MQ=255) ggACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTTAC‑TTTTTTTAAAAACCACACCCTC‑CCGGGGCAttttttttttttattg > 2:711390/1‑134 (MQ=255) gACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTTAC‑TTTTTTTAAAAACCAAACCCTC‑CGGGGAGATTTTATTTTTTTATTGTATcg > 1:813525/1‑151 (MQ=255) ggAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCa < 1:937619/151‑1 (MQ=255) aCTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAAA‑TTTTTTTAAAACCCCAACCCTC‑CGGGGACAATTGTTTATTTTATTGTATtgccgcca > 1:869547/1‑147 (MQ=255) cTTTAACGACAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGCCATTTATC‑TTTTTTTTAAACCCCAACCCTC‑CGGGGGGCTTTGATTTTTTTCTTGTATCTGCGCgaaa > 2:256218/1‑150 (MQ=255) ttAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAAC‑TTTTTTTAAAACCCCTACCCCC‑CGGGGGCAATTTATTTTTTTATTGTATTGCCCCaaaa > 2:150276/1‑148 (MQ=255) tAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAC‑TTTTTTTAAAAACCACAACCTC‑CGGTGGCATTTTATTTTTTTTTGGTATCGCCCCcaaaaacc > 1:585752/1‑147 (MQ=255) aaCGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTccc < 1:673027/151‑1 (MQ=255) aaCGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTcc < 1:468178/150‑1 (MQ=255) cGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGCTTTAAC‑TTTATTTTCAACCCACACCCTC‑CCGGGAAATTTTTTTTTTTTATCGTATTGCGGCAAAACTccccc > 2:176942/1‑149 (MQ=255) gAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCccgc < 2:590361/150‑1 (MQ=255) gAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGGTTTTACTTTTTTTTAAAACCCA‑ACCGTC‑CCGGGGCATTTTTTTTTTTTATTGTATTGCCGCcaaaaacccccc > 2:551421/1‑143 (MQ=255) aaaaaCCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTAA‑TTTTTTTAAAAACCATACCGTC‑CCGGGGCATTTGTTTATTTTATGGTATCGCCGCCAAAATCCcgccgc > 1:307418/1‑151 (MQ=255) aaaaaCCGCGCGTCTTTAGACTGGCTATTCGTTGTTGACCGCCGTTGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAGACCCGTACCGGC‑CGGTGACTTTTGATTATTTAATGGTATCGCCGCAAAACTCccgc < 2:650386/148‑1 (MQ=255) gcgcgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTTAC‑TTTTTTTAAAAACCCTACCCCC‑CGGCGACATTTGATTATTTTATTGTATTGCCCCCAAACTCCCCCCCCCCAtttt > 2:922544/1‑149 (MQ=255) cgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTAt < 2:895740/150‑1 (MQ=255) gTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTGTTTTAC‑TTTTTTTAAAAACCATACCATC‑CGGTGAAATTTTTTTTTTTAATGGAATTGCCGCAAAACACCCCCCGCCCCTTGt > 1:876621/1‑147 (MQ=255) ttATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGaa < 1:115574/148‑1 (MQ=255) ttATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAAC‑TTTTTTTAAAACCCATACCCGC‑CCGGGGCATTTTGTTTTTTTATTGGATTGCCCCAAAAAACCCCCCGCCCATTTTATGaaa > 1:550494/1‑149 (MQ=255) taCTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTTAC‑TTTATTTAAAACCCCAACCCGC‑CCGGGACATTTTATTTTTTTATGGGAATGCCGCAAAACTCCCgccgcccgttttttgtaa > 1:407029/1‑134 (MQ=255) taCTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTTAA‑TTTATTTAAAACCCCAAACGCCACGGTG‑CATTTTGTTTTTTTATGGGTTACCCCCAAAAACCCCGCCCCCCATTGTATGAAAtggtg > 2:1065585/1‑151 (MQ=255) actaactaTTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAAtggtg < 2:691510/150‑1 (MQ=255) tGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAAC‑TTTATTGAAAAACCATACCCTC‑CGGGGGCATTTGATTTTTTTATGGGATCGCCGCAAAACTCCCCCCGCCTATTGTATGAAAtggtgg > 1:741491/1‑149 (MQ=255) aCTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTTAC‑TTTATTTTAAACCCATACCCTC‑CGGTGGCATTTGGTTTTTTTATTGTATCGCCGCAAAAAACCCCCCGCCCATTGTATGAAAtggtgg > 2:437135/1‑147 (MQ=255) aCTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACcg < 1:96591/151‑1 (MQ=255) cTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGAcc < 2:954634/149‑1 (MQ=255) ttCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTTAC‑TTTTTTGTAAAACCCTACCGTC‑CCGGGACATTTGATTATTTTATTGTATTGGCGCCAAAATACCCCCGCCCACTGTATGaaaatgggggccg > 2:382009/1‑137 (MQ=255) cGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCgag‑a < 1:630736/150‑1 (MQ=255) cGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAG‑At < 1:29361/151‑1 (MQ=255) cgccgATGCGCGCATCACCAACGAAGTACGGCGGGTGTTTTACATTTTTTTTAAACCCAT‑TCGCC‑CGGGGGAAATTTTTTATTTTATTGTATTTACCCAAAACTCCCGCCGCCCCCTGTATTGAATTGCGGAACGCGGGAATCGCCacac > 2:862013/1‑147 (MQ=255) cgATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATAGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAG‑ATCGCCAACGCgcag < 2:865547/151‑1 (MQ=255) cgATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTAC‑TTTTTTTTAAAACCCTACCCTC‑CGGGGGCATTTGATTTTTTTATTGTATTGCCCCAAAAATCCCCCCCCCCCTTTTATTAAAAAGTGGACCGCGGG‑AGCGCCCAccc > 2:583049/1‑145 (MQ=255) cgcATCACCAACGAAGTACGGCGGGTATTTTAC‑TTTATTTAAAACCCCTACCCTC‑CGGGGGCATTTTATTATTTTATGGGATCGCCGCAAAAATCCCGCCCCCCCTTTTTTGAAAATGGGGGCCGCGAG‑ATCGCCAACCCGCAGCAAagg > 1:333263/1‑150 (MQ=255) aCCAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAG‑ATCGCCAACGCGCAGCAAGGGCTGCCCa < 1:579747/151‑1 (MQ=255) aCCAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAG‑ATCGCCAACGCGCAGCAAGGGCTGCCCa < 1:797712/151‑1 (MQ=255) ccAACGAAGTACGGCGGGTATTTTAC‑TTTATTTAAAACCCAAACCCCC‑CGGGGGCAATTTTTTTTTTTATTGTATTGCCCCAAAAATACCCCCGCCTATTTTATGAAAAGGTGGGCCGCGGG‑AACGGCCACCCCCAGCCAAGGccggcc > 2:1048211/1‑144 (MQ=255) ccAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAG‑ATCGCCAACGCGCAGCAAGGGCTGCCCAg < 1:190472/151‑1 (MQ=255) cAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAG‑ATCGCCAACGCGCAGCAAGGGCTGCCCAGt < 1:296449/151‑1 (MQ=255) cggcggGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAG‑ATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCAccc < 2:813401/151‑1 (MQ=255) cggGGATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAG‑ATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTg > 1:711269/1‑150 (MQ=255) gggTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAG‑ATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGaa < 1:757480/151‑1 (MQ=255) aC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAG‑ATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATa > 2:298118/1‑150 (MQ=255) ttGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAG‑ATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGt > 2:331468/1‑151 (MQ=255) tGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAG‑ATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGTc > 2:443569/1‑151 (MQ=255) tGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAG‑ATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGTc > 1:530136/1‑151 (MQ=255) | GTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAC‑TTTATTGAAAACCCATACCGTC‑CGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAG‑ATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGTC > NC_000913/2624371‑2624666 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |