Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,378,802 | 0 | A | C | 49.0% | 72.8 / 14.5 | 51 | V48G (GTA→GGA) | ampC | penicillin‑binding protein; beta‑lactamase, intrinsically weak |
Reads supporting (aligned to +/- strand): ref base A (25/1); new base C (0/25); total (25/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.09e-13 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGG > NC_000913/4378657‑4378948 | tCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑accgcc > 2:602256/1‑151 (MQ=255) tCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑ATTACCC‑T‑GATAAATT‑accgcc > 2:759644/1‑151 (MQ=255) cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑accgcc > 1:499622/1‑150 (MQ=255) gCAATAGCGTCGCCACCAACCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑accgccacc < 2:972362/150‑1 (MQ=255) cAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑accgccaccgc > 2:683396/1‑151 (MQ=255) gTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATAccc > 1:814009/1‑151 (MQ=255) gTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATAccc > 2:734220/1‑151 (MQ=255) gTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATAccc > 1:391647/1‑151 (MQ=255) tCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCCg > 2:141746/1‑151 (MQ=255) ccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCCGGGATCt > 2:204777/1‑151 (MQ=255) cgcaataaattGTTTTGCTGCCCGACCCAAACAAAAACAATTTTTGTTGTGTGAGGGTTTGTTTTTGGGGGTTGCCCGCATA‑CCCCCGG‑TAAAAGAAAAAAG‑TTTCCCCTT‑AAAAAATC‑CCCCCCCCCGCCATACCCGGGATCTTttgtt < 1:783327/140‑1 (MQ=255) ccccATTAATTTTTTCGCTGCCCAACCCAAACAAAAAAAGCGTTGGTTGTGTGCGGGGGTTTTTTTTGGCGATCCCCGCACGCCCCCAGA‑AAAAGAAATAGGT‑TTTCCCG‑T‑AATAAATA‑CCCGCCCCCGCCATACCCGGGATCTTttgttg < 2:692850/149‑1 (MQ=255) cGCCAATAATTTTTTTCCTGCCCGAACATACCTCAAAAAACTTTTGCTGTGTGACGGGTTTTTTTTGGGCGATCTCCGCATCGCCCAAGA‑TAAAGTAA‑AAGG‑TTTTCCCTT‑AATAAATT‑CCCGCCACCGCCATACCCGGGATCtttt < 1:387525/147‑1 (MQ=255) gTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGctct > 1:707514/1‑150 (MQ=255) ttttttGCGGACCAAACCAACCAAAAACAACGTTTGCTGTGGGGGGGGTTGTTTTTTGGCGAGGCCCCCCCCCCCCCAGA‑AAAAGAAAAATTT‑TTTCCCA‑T‑AATAATTC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATaa < 2:32283/147‑1 (MQ=255) tGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGc > 1:810593/1‑151 (MQ=255) accaaccttaacaaaaaCAACTTTTGGTGTGTGGGGGGTTGTTTTTGGGCGATGCCCGCATACCCCAAGA‑AAAAGTAATAGTT‑TTTCCCC‑T‑AAAAAATC‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCgggg < 2:72501/131‑1 (MQ=255) ccgacccaaccaaaaaCAGCTTTTGTGGTGTGAGGGGTTTTTTTTGGGCGATGCCCGCATACCCCAAGA‑AAAAGTAATAATT‑TTTCCCT‑T‑AAAAATTT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTaa < 2:896170/134‑1 (MQ=255) ccaaaCCAAACACAACCATCGTTGGTTGTGGGAGGGGCTGCTTTTTGGCGATGCCCGCATAGCCCCAGG‑AAAATTAATAGTT‑TTTCCCC‑T‑AATAAATT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtg < 1:201432/147‑1 (MQ=255) ccGAACCTAACTAAAACACCGTTTGCTGGGTGACGGGGTTTTTTTGGGCGATTTCCGCACAGCCCGAGG‑AAAAGTAATAAGG‑TTTACCC‑G‑AATAAATT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGt < 1:119968/145‑1 (MQ=255) acccaaacaaaaaCAACTTTTGGTGGGTGAGGGGCTTTTTTGTGGGGATGCCCGCATCCCCCGAGA‑AAAAGTAATAGGT‑TTTCCCC‑T‑AAAAAATT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCg < 2:125021/141‑1 (MQ=255) ccTAACACAAACAACTTTTGCTGTGTGGGGGGTTGTTTTTGGGGGATGCCCGCATA‑CCCCAGGAAAAAGTAATAAGT‑TTTCCCC‑G‑AAAAAATT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGAt < 1:916946/151‑1 (MQ=255) tCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcaca > 2:99911/1‑150 (MQ=255) aaaaCAATGTTTCTTGTGTGGGGGGCTGTTTTTTGGTGTTGCCCTTTAATCCCCAGG‑TAAAGTAA‑AAGG‑TTTCCCC‑T‑TAAAAATTCCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat < 1:765855/150‑1 (MQ=255) aaaCAACTTTTGTTGTGTGAGGGTCTTTTTTTTGGCGATCTCCCCATACCCCCAGG‑TAAAGTAATAGTT‑TTTACCT‑T‑AATAATTC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat < 2:65104/150‑1 (MQ=255) caccttttGTTGTGTGGGGGGCTTTTTTTGGGGGATGTCCCCATCGCCCCAGG‑AAAAGTAAAATGT‑TTTCCCC‑T‑AAAAATTT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATc < 2:399453/145‑1 (MQ=255) ttGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAtt > 2:872201/1‑151 (MQ=255) gtggtgtgACGGGTTTTTTTTGGGCGTTTCCCGCTTCCCCCCAGA‑AAAATTAATATTT‑TTTCCTC‑T‑AAAAATTT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAttt < 2:577782/147‑1 (MQ=255) gtgtgCCGGGCTTTTTTTTGGCGATGCCCGCATACCCCCAGG‑‑AAAATAAAAAGGTTTTCCCC‑T‑AATAAATT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgt < 1:6009/149‑1 (MQ=255) tgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCAGGGATCTTTTGTTGCTCTATAAGCGGGGCAATTGTGCGATGCACAATATCGTTGATttgttg > 2:774400/1‑150 (MQ=255) cGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc > 2:595069/1‑151 (MQ=255) tttttGGCGATGTCCCCACACCCCGAGT‑AAAAGTAATATGT‑TTTCCCT‑A‑AATAAATC‑CCCGCAACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAt < 1:583530/150‑1 (MQ=255) ttttgGGGGATGTCCGCATACCCCCAGT‑AAAAGAAAAAAGT‑TTTCCCC‑T‑AAAAATTT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCaa < 1:123504/148‑1 (MQ=255) tgggCGAGGCCCGAAAACCCCCAGA‑AAAAGTAATAGGT‑TTTCCCA‑T‑AAAAAATC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgtg < 1:869698/148‑1 (MQ=255) cGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGc > 2:90673/1‑150 (MQ=255) cGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGc > 1:269511/1‑150 (MQ=255) tGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa > 1:514181/1‑148 (MQ=255) tCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAgg > 2:796283/1‑151 (MQ=255) cGCATGCCCCCGGG‑AAAAGTAATAAGT‑TTTCCCC‑A‑AAAAATTT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCg < 1:914394/151‑1 (MQ=255) tAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGta > 2:225010/1‑150 (MQ=255) cccccGGT‑TAAAGAAAAATGT‑TTTCCCC‑T‑AATAATTT‑CCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt < 1:213305/150‑1 (MQ=255) ccccAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt > 2:756685/1‑150 (MQ=255) ccccAGG‑AAAATAAATAATT‑TTTCCCC‑T‑AATAATTC‑CCCCCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatta < 2:433286/151‑1 (MQ=255) gg‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaat > 1:722090/1‑148 (MQ=255) aaataaaaaAG‑TTTTCCC‑TGGATAATTC‑ACCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg < 2:828166/148‑1 (MQ=255) ataaGT‑TTTACCT‑T‑GATAAATA‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag < 1:759729/148‑1 (MQ=255) aaGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCAGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTACATAAAAAGACGCAGAGcg > 2:213338/1‑150 (MQ=255) ttttCCC‑T‑AAAAAATA‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGc < 2:907818/145‑1 (MQ=255) tttACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtt > 2:239765/1‑151 (MQ=255) ttACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGttt > 2:133969/1‑151 (MQ=255) tttcccttaaTAATTT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtcg < 2:813928/140‑1 (MQ=255) tACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtttt > 2:539597/1‑151 (MQ=255) t‑GATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAAc > 1:158114/1‑150 (MQ=255) aTAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt > 1:831768/1‑150 (MQ=255) aaaTT‑CCCCCCACCGCCATTCCCGTGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt < 1:183323/148‑1 (MQ=255) aattcaccgccaccgccATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATa < 2:499571/144‑1 (MQ=255) tt‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg > 2:112720/1‑149 (MQ=255) | TCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGG‑TTTACCC‑T‑GATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGG > NC_000913/4378657‑4378948 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |