Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,624,517 | 0 | G | T | 41.4% | 75.5 / 26.6 | 58 | E492* (GAA→TAA) | ppk | polyphosphate kinase, component of RNA degradosome |
Reads supporting (aligned to +/- strand): ref base G (7/27); new base T (24/0); total (31/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.48e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGTC > NC_000913/2624375‑2624666 | aGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCa < 1:55950/151‑1 (MQ=255) aGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCa < 1:662706/151‑1 (MQ=255) aGAAAACGGTAAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCa < 1:547517/151‑1 (MQ=255) gTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTcc < 1:339463/151‑1 (MQ=255) gTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTcc < 2:370316/151‑1 (MQ=255) gAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCg < 1:86193/150‑1 (MQ=255) gCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACAttt < 1:342450/151‑1 (MQ=255) gTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGa < 1:305043/151‑1 (MQ=255) tATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTTAATTTTTTTAAAACCCAT‑A‑CCGccccggggcgttttttt > 2:377007/1‑137 (MQ=255) tATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTTAATTTTTTTAAAAACCAT‑A‑ACGTCCCGTGGCAttt > 2:925342/1‑147 (MQ=255) aTGCACACATCGGGGCGGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACTAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGatta < 2:898286/151‑1 (MQ=255) aTGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCTCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGatt < 1:823548/150‑1 (MQ=255) gCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTAACTTTTTTTAAAAACCAT‑A‑CCGTCCGGTGGCATTTGattatt > 2:972903/1‑151 (MQ=255) cacaccTGGTGACCGGGACCTTTAACGAAAAATCCGCGCGTCTTTATACTGACTATTCGTTGCTGCCCGCCGATGCGCGCATCCCCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGattatt < 1:713750/150‑1 (MQ=255) gggACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTTAATTTATTTGAAAACCAA‑A‑CCGTCCGGGGGCCTTTGATTTTTTTATTGTAt > 1:835697/1‑151 (MQ=255) gggACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTAACTTTATTTACAACCCAT‑A‑CCGTCCGGGGACATTTGATTATTTTATTGt > 1:172520/1‑149 (MQ=255) ggACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATc < 1:940156/151‑1 (MQ=255) ggACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGT‑TTTTAATTTTTTTTAAACCCCT‑A‑CCCTCCGGGGACATTTGTTTATTTTATgg > 2:919331/1‑147 (MQ=255) gACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTTACTTTTTTTAAAAACCAC‑A‑CCGTCCGGGGGCATTTGATTATTTTATGGTATcg > 2:107799/1‑151 (MQ=255) ccGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATcgcc < 1:433594/151‑1 (MQ=255) ggAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCa < 1:546023/151‑1 (MQ=255) aaCTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCaaa < 2:654617/151‑1 (MQ=255) aaCTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTTACTTTTTTTAAAACCCAT‑A‑CCGTCCCGGGGCATTTGATTATTTTATTGTATCGCCGCaa > 2:493947/1‑150 (MQ=255) cTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGATTTTAATTTTTTTTAAAACCC‑T‑A‑CCGGCCGGGGGGAATTTTTTTTTTAATTGTATTGCCCCcaaa > 2:386182/1‑150 (MQ=255) tttAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTTAATTTTTTTAAAAACCCC‑A‑CCCTCCGGGGGCATTTTATTTTTTTATGGGATCGCCGCcaaa > 2:894849/1‑149 (MQ=255) tAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTc < 1:813317/150‑1 (MQ=255) aaCGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTccc < 2:304142/151‑1 (MQ=255) cGAAAAAACCGCGCGTCTTTAGTCTGACTATTCGTTGCTGACCGCCGATGCGCGCATCTCCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCccgc < 1:792911/151‑1 (MQ=255) aaaaacCCCCGCGTCTTTATACTGACTATCCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCccgcc < 1:220282/150‑1 (MQ=255) cgcgcgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTAACTTTTTTTAAAACCCCT‑A‑CCCTCCGGGGGCACTTTTTTTTTTTATTGTATCTCCCCAAAACcccccccccccctt > 1:322217/1‑146 (MQ=255) gcgcgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTTTTGAAAACCCCC‑A‑CCCGCCGGGGGCATTTTTTTTTTTTATGGTATCGCCGCAAAACTCCCGCCGCCTCtt > 1:415017/1‑149 (MQ=255) cgcgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTTACTTTTTTTAAAAACCATAA‑CCGCCCGG‑GGCATTTGATTATTTTATGGTATCCCCGCCAAAATCcccccccccttttt > 2:374856/1‑144 (MQ=255) gcgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTAt < 1:95861/151‑1 (MQ=255) gTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTT‑ATTCAATTTATTGAAAAACCAT‑ACCCCCCCGGCGA‑AATTGCTTTTTTTATGGTTTAACCCCCCAACTCCCGCCCCCTACTGttt > 2:805673/1‑147 (MQ=255) tataCTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGTGG‑TTTTTATTTTTTGAAAAACCAC‑A‑CCCTCCCGGGACATTTTTTTTTTTTATTGGATTTCCCCAAAAAACCCCCCCCCTAtttttttgaa > 1:563421/1‑140 (MQ=255) cTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGTT‑TTTACCTTTTTTTAAAACCCCA‑C‑CCCCCCGGGGCAATTTG‑TTTTTTTATGGGTTCGCCCCAAAACTCCCCCCCCCTAttttttttaatgggggg > 1:16318/1‑136 (MQ=255) tATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTTTTTAAAACCCAA‑A‑CCGCCCGGGGGCATTTGATTATTTTATGGTATCGCCGCCAAAATCCCCCCCCCCATTGTTTGAAATGGTGGGcc > 1:969369/1‑148 (MQ=255) ttGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGG‑TTTTACTTTTTTTAAAACCCCT‑A‑CCCCCCGGGGGCATTTTATTTTTTTATTGTATCGCCGCAAAACTCCCCCCGCCCCTTTTTTGAAATGGGGGGCCGCGAGGTc > 1:533587/1‑150 (MQ=255) gCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTTAATTTTTTTAAAAACCAA‑A‑CCCTCCGGGGGCATTTGTTTATTTTATGGTATCGCCCCCAAACTCCCCCCCCCCATTGTATGAAAAGGGGGACCGCGGGAACGcc > 2:624779/1‑151 (MQ=255) tGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTAACTTTATTTAAAAACCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTTTATGAAATGGTGGACCGCGAGATCGcc > 2:897135/1‑149 (MQ=255) gACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTTACTTTTTTTAAAACCCAT‑A‑CCGCCCGGTGGCATTTGATTTTTTTATTGTATTGCCGCCAAACTCCCCCCCCCTATTGTATGAAAATGTGGACCGCGGGGTCGCCCAc > 1:895130/1‑151 (MQ=255) gACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCaa < 1:7312/150‑1 (MQ=255) aCCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTTACTTTATTTAAAACCCAA‑A‑CCGCCCGGTGGCATTTTATTTTTTATTGGAATCGCCGCCAAAATCCCGCCCCCCATTTTATAAAAATGTGGGACCCGAGACCGccca > 2:104738/1‑147 (MQ=255) cgccgATGCGCGCATCACCAACGAAGTACGGC‑GGGGA‑TTTTACTTTTTTTAAAAACCAT‑A‑CCCTCCGGGGGCATTTTATTTTTTTATGGGATCGCCGCAAAACTCCCCCCGCCCCTTGTTTGAAATGGGGGACCGCGGGATCGCCCAccc > 2:655918/1‑148 (MQ=255) tGCGCGCATCACCAACGAAGTACGGC‑GGGGT‑GTTAACTTTTTTTAAAAACCAT‑A‑CACTGCGGGGGGAATTTATTATTTAATTGGATCTGCCCAAAAAACCCGGCCCCCCTTTTATGAAAaggggggacccggggttg > 2:886663/1‑125 (MQ=255) cgcgcATCACCAACGAAGTACGGC‑GGGTA‑TTTTAATTTTTTTAAAAACCAT‑A‑ACGTCCGGGGGCATTTGTTTTTTTTATAGGATCGCCGCAAAAATCCCGCCCCCTCTTGTATGAAATGGGGGACCGCGAGAACGCCAACGCGCAGCaa > 2:401489/1‑149 (MQ=255) gcATCACCAACGAAGTACGGC‑GGGTG‑TTTTATTTTTTTTTAAACACAA‑A‑CCCTCCCGGGGCATTTGGTTATTTAATTGTATTGCCGCAAAACTCCCCCCCCCCATTTTTTTAAAATGGGGAGCGCGGGGGCGCCAACCCCCCGCAAggg > 1:624026/1‑149 (MQ=255) cATCACCAACGAAGTACGGC‑GGGTA‑TTTTACTTTTTTTAAAAACCAT‑A‑ACGTCCCGGGGCATTTTATTTTTTTATTGTATCGCCCCAAAAAACCCCCCCCCCCTTTTATGAAAAGGTGGACCCCGAGGTCGCCCACCCCCCGCAAggg > 2:603726/1‑148 (MQ=255) ccAACGAAGTACGGC‑GGGGA‑TTTAACTTTATTGTAAAACCAT‑A‑CCGCCCGGGGGCATTTGATTATTTTATTGTATTGCCGCAAAAATCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCCGCAAGGGCTGccc > 1:74359/1‑149 (MQ=255) cAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGt < 1:932659/151‑1 (MQ=255) aaCGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTg < 2:698068/151‑1 (MQ=255) gAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTa < 2:649177/151‑1 (MQ=255) aGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATc < 2:733779/151‑1 (MQ=255) gggTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACAGTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGaa < 2:585908/151‑1 (MQ=255) aaCTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAAt > 2:312962/1‑150 (MQ=255) aCTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATa > 2:71426/1‑150 (MQ=255) tttATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAAcc > 2:180087/1‑151 (MQ=255) tGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGt < 2:449433/150‑1 (MQ=255) gAAACCCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGTc < 1:226399/150‑1 (MQ=255) | AGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTA‑TTTAACTTTATTGAAAACCCAT‑A‑CCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGTC > NC_000913/2624375‑2624666 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |