Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,601,176 | 0 | T | C | 42.9% | 41.9 / 16.8 | 37 | E304G (GAA→GGA) | ftsX | putative ABC transporter permease |
Reads supporting (aligned to +/- strand): ref base T (19/1); new base C (0/15); total (19/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.93e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTA > NC_000913/3601031‑3601322 | ttCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGc > 2:474831/1‑150 (MQ=255) cGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGCAAccgcc < 2:126828/151‑1 (MQ=255) aaaGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGCAACCGCCGAt > 2:864258/1‑151 (MQ=255) aaaTGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGCAACCGCCGATGACAATCgc > 2:12530/1‑150 (MQ=255) aaTGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGCAACCGCCGATGACAATCgc > 2:881181/1‑149 (MQ=255) aaTGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGCAACCGCCGATGACAATCgc > 1:635300/1‑149 (MQ=255) ccGGGGCAACCCCCGCGCCCCCCCACAAAAACATCGAGCATACCAGCAAAAAAAGCAGGCATTCATAAAATGAAAAGCTTTTAAAACCAACCTTCGCCCCAAAAC‑CC‑TGTC‑CCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAAt < 2:580940/151‑1 (MQ=255) gTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAAttt > 2:722556/1‑151 (MQ=255) aaGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGa > 2:539153/1‑150 (MQ=255) gCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACa > 2:602781/1‑150 (MQ=255) cccccccccccccccacccaACCCTCGAGCATCCCAGCAGAAAGAGCGGGCATTCATCGAATAAAAAGCTATAAAAAAAAACCTCCGTCCCAAAAC‑CC‑GGTC‑CCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTg < 2:426377/135‑1 (MQ=255) ccccgcccccccccACCAAATCATCGAGCACCCCAGCAACAATAGCGGGTTTTCACCAAATGAAAAGCCTTTAATATAAAACTTCTTCCCAAAAC‑CC‑TGCC‑CCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACa < 2:928071/138‑1 (MQ=255) gCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAAtt > 2:432621/1‑150 (MQ=255) ctcccACCGAGCAAACCACCGAGCAGCCCAGAACAAATAGCAGACTTTCATCAAATAATAACCCATGTAAATAAAGCTTC‑TCCCCAAAA‑CCTTGTC‑CCTCTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAAt < 2:444497/145‑1 (MQ=255) cccccaccaACCAAATCAGGGAGCACCCCAGCAGCAATAGCAGTCTTCCAGAGAATAAAAAACCATTTAAAAAAAACTCCGTCCCAAACC‑CC‑TTTC‑CTATCCACGCACCCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAAtt < 1:285806/146‑1 (MQ=255) ccccccAGACAATCATCGAACACCCAAGCAGAAATAGCAGCCTTTCAGGAAATAATAACCCTTTAAAAAAAAATTTCGTACAGACAT‑CC‑TGTG‑CAATTACCACATCCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGa < 1:146199/148‑1 (MQ=255) ccccccACCAAATCACCGAACACCCCACCAGAAATAGCGGCCTTTCAACGAATAAAACCCCATTGATAAAAACCTTCTTTCCAAAAC‑CC‑TGTC‑CCCCTCCCAAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTaa < 2:417471/145‑1 (MQ=255) aaTCATCGAGAAAACCAGCAGAAAGAGCGGGCATTCATCAAATAAAAAGCCATTGAAATAAACCTCCGTCCCAAAAA‑CC‑TTTG‑CCCCCACCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATgc < 1:600911/150‑1 (MQ=255) aaTCACCGACCAACCCACCAGCAAAAGCAGATTTTCAAAAAATAATACCCCATTAATAAAAACCTCCGTCCCAAAAC‑CC‑TGTG‑CCCCTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATg < 1:428662/149‑1 (MQ=255) atcatcGAAAACACCAGAAGCAAACGCAGGATTTCATCAAATAATAA‑CCATTTAAAAAAAATCTCCGCACAAAAA‑CC‑TTTGCCCCCTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATgc < 1:766385/149‑1 (MQ=255) caccGAACATACCACCAGCAAAAGCAGGCATTCATCAAATAATAACCAATTAATATAAACCTTCGTCCCAAAAC‑CC‑TGTC‑CCACTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATgcg < 1:146216/145‑1 (MQ=255) accccaacaaaaaTAGCAGCCATCCATCAAATAATAACCCTTTAAAAAAAACCTTCTTCCCAAAAC‑CC‑TGGC‑CCTCTCCCCAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAg < 2:345373/133‑1 (MQ=255) aTACCAGAAGAAATAGCGGGCATTCATCGAATAAAAAGCTATTGATATCAACCTGC‑TTCCAAAAACCC‑TGTC‑CCTCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGaa < 2:913656/146‑1 (MQ=255) tACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATccc > 1:872303/1‑151 (MQ=255) aaTAGCGGCTATTCATCAAATAAAAACCAATTAAAAACAATCTTCATCCCAAAAC‑CC‑TGCC‑CCCCTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTgcg < 1:143481/150‑1 (MQ=255) aGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGccac > 2:348649/1‑151 (MQ=255) ggCATTCATAAAAAGATAAGCTATTAATATCAATCTCCGTACAAAAAC‑CC‑TGTC‑CCCCTCCCAAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCAt < 1:742419/150‑1 (MQ=255) gAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACCTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACggg > 1:264498/1‑151 (MQ=255) ccATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATg > 2:654281/1‑150 (MQ=255) ccATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATg > 1:862112/1‑150 (MQ=255) taataccaaCCTCCTTCCCAAAA‑CCTTGTC‑CCCCTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATc < 2:586959/143‑1 (MQ=255) cGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCAcc > 1:858192/1‑150 (MQ=255) ttCCGAAAA‑CC‑TGTG‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAAt > 2:854179/1‑151 (MQ=255) tCCGAAAA‑CC‑TGTG‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATc > 2:595235/1‑151 (MQ=255) gAAAA‑CC‑TGTG‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGt > 2:781989/1‑151 (MQ=255) aaaa‑CC‑TGTG‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGtt > 1:247238/1‑151 (MQ=255) cACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTa < 1:357694/151‑1 (MQ=255) | TTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTG‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTA > NC_000913/3601031‑3601322 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |