Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,356,201 | 0 | T | G | 52.9% | 34.1 / 51.8 | 34 | T304P (ACC→CCC) | sapA | antimicrobial peptide transport ABC transporter periplasmic binding protein |
Reads supporting (aligned to +/- strand): ref base T (2/14); new base G (18/0); total (20/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.62e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.33e-06 |
GTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCG > NC_000913/1356055‑1356343 | gTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑ATTTGCGGGG‑tt > 1:622883/1‑150 (MQ=255) gCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGG‑GTAAATGccc > 1:428019/1‑147 (MQ=255) cccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATa < 1:106562/151‑1 (MQ=255) cGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCg < 1:472460/151‑1 (MQ=255) gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGG‑GTTAATGCCCGAAAGGCGGCGTTctttcc > 1:272594/1‑146 (MQ=255) tAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGTG‑TTTAATACCAGATAGGGGACGTTCATccc > 1:360684/1‑150 (MQ=255) tAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGG‑GTTAATGCCCGAAAGGGGGCGTTCATtcc > 2:584527/1‑145 (MQ=255) cgccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGtt < 2:5918/151‑1 (MQ=255) gccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGGGGGGTTTAAATGCC‑GATAGGCGGCGGTTATTCCAgg > 1:330948/1‑140 (MQ=255) ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCgg < 2:110455/149‑1 (MQ=255) ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGgt < 2:510702/151‑1 (MQ=255) aataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGG‑TTTAATGCCAGATAAGCGACGTTTATTCCCGGACGCAGCGGTTAAAGCAAGCGCggg > 1:282402/1‑146 (MQ=255) ataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTTAAAGCCAGATAAGCGGCGGTTCTCCCCGGGCGCAGCGTTTAAAGCCAACGCggg > 1:644599/1‑145 (MQ=255) taGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGG‑TTAAAT‑CCAAAAAGGGGGGGGTCATCCCAAGGAGCAGCGGTAAAAGACAACGCgggg > 2:652783/1‑144 (MQ=255) tGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGG‑TTAAATGCCAGATAGGCGACGGTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGGCGTCACGc > 1:288487/1‑149 (MQ=255) gggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCaa < 2:424634/150‑1 (MQ=255) gCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGG‑TTAAAAGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGGTAAACGCAAGCGCGGGTCGTCACGCAAAAAgg > 1:425183/1‑150 (MQ=255) tGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGGGGGG‑TTAAATACCCAATAGGCGGCACTCTTCCCACGGAGCAGCGGTAAACGGCAGGGCGGGGGGGCCCGCAAAAAgggaagacg > 2:664005/1‑139 (MQ=255) ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGctggctgg < 2:196333/150‑1 (MQ=255) ttattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGG‑TTAAATGCCAGATAGGCGACGTTCATCCCCGGACGCAGCGTTAAACGCAAGCGCGGGGCGTCACGCAAAATGGATAGctggctgggcg > 1:625890/1‑148 (MQ=255) ttattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGG‑GTTAAAGCCCGATAGGCGACGTTCATCCCAGGGCGCAGCGGTAAAAGCAAGAGCGGGGCGTCACGCCAAAAGGATAGatgggtgg > 1:57275/1‑148 (MQ=255) tattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGG‑TTAAATTCCAGATAAGCCACGGTCATCCCCGGCCCCAGCGTTAAACGCAAGCGCGGGGGGTGACGCAAAAAGGATAACtggctggccg > 1:650982/1‑148 (MQ=255) ttaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGG‑GTAAATGCCAGAAAAGCGACGGTCATCCCCGGACGCCGCGTTAAAAGCAAGCGCGGGGCGGCACGCAAAAAGGATAGCtggctggccg > 2:144759/1‑146 (MQ=255) cGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGG‑GTAAATGCCCGGAAGGGCGCGTTCCTCCCACGGCGCCGGGGTTAAAGCAAGCGCGGGGCGTCACCCAAAAAGGGCAAGCGGCTGGCCGcagg > 2:602454/1‑147 (MQ=255) gCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGaa < 2:470326/151‑1 (MQ=255) aTGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGG‑GTAAATGCCCGATAGGGGACGGTTATTCCAGGGCGCAGCGGTTAACGCAAACGCGGGGGGGTCCGCAAAAAGGAAAGGCGGGTGGGAGCCGGGCCAGCCCGAAAGTCGCa > 2:243177/1‑150 (MQ=255) aTGGCGGACAGCGGGATTATTTAGCGGCGG‑GTTTGCGGGG‑GTAAATGCCCGATAAGCGACGTTTATCCCCGGACGCAGGGGTTAACGCAAACCCGGGGGGTCACGCCAAAAGGGGAGGCGGCTGGGAGGAGGGCCGGGCAGAAAGgcgc > 2:692411/1‑149 (MQ=255) cGGACAGCGGGATTATTTAGCGGCGG‑TTTGGGGGGG‑TTTAAAGCCAGATAGGCGGCGTTCATCCCCGGGCGCAGCGTTAAACGCAAGCGCGGGGCGTCACGCCAAAAGGATAGgtgggtgggcg > 1:517529/1‑121 (MQ=255) cGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGcagcag < 2:517464/127‑1 (MQ=255) aCAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCgg < 2:317390/151‑1 (MQ=255) cGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAg < 2:644639/151‑1 (MQ=255) ggATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGa < 1:758466/151‑1 (MQ=255) gATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAg < 2:733631/151‑1 (MQ=255) ttGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCg < 1:201024/150‑1 (MQ=255) | GTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCG > NC_000913/1356055‑1356343 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |