| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,601,176 | 0 | T | C | 47.7% | 40.8 / 24.6 | 44 | E304G (GAA→GGA) | ftsX | putative ABC transporter permease |
| Reads supporting (aligned to +/- strand): ref base T (21/2); new base C (0/21); total (21/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.26e-10 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCC‑ATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATG > NC_000913/3601026‑3601325 | tattatTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCC‑ATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACtt > 2:21820/1‑151 (MQ=255) attCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCC‑ATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCg < 2:153363/150‑1 (MQ=255) gTAAAGTGGCATAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCC‑ATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAAccgcc < 1:839518/150‑1 (MQ=255) aaTGTTGCCCCGGGCAAAGCCCCGCCGCCCCCCAACAATACACAGAGAATACCAGCAAAAATAGCAGTCATTAATCAAATGAAAAGCT‑TTTGAAAACAACCTTC‑TTCCCAAAA‑CCAGGCCCCACTCCCGAAACCGCCGATGACAATCgcag < 1:128323/151‑1 (MQ=255) ccGTGGAAAGCCCCGCTCCCCCCCAGCAAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCAAATGATAAGCC‑ATTGATAACAACCTTCGTCCCAAAAC‑CC‑TGTCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAAt < 2:872874/151‑1 (MQ=255) gggCAAGCCACGCTCCCCCCCAGCCAATCATCGAGCATCCCAGCAGAAATGGCAGGCTTTCATCGAATAATAAGCC‑ATTGATATAAAACTTCTTTCCAAAAC‑CC‑TGTGCCACTCCCCAAACCGCCGATGACAATCGCAGCACCAGAAtt < 2:311584/148‑1 (MQ=255) ccACCCTCCCCCCCAGCAAATCACCGACCATCCCAGCAAAAATAGCAGGCATTCATCAAATGATAACCC‑ATTGATACCAACCTCCGTCCAAAAAC‑CC‑TGTGCCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACaaa < 1:825906/151‑1 (MQ=255) ccgCCCCCCACCCAATCATCAACCACCCAAGCAAAAATAGCAGGCATTCATCAAATGATAACCC‑ATTGATATCAACCTTCGTCCCAAAAC‑CC‑GGTCCCACTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTa < 1:872721/148‑1 (MQ=255) aCCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCC‑ATTGATATCAAACTTCGTTCCGTAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGac > 2:945444/1‑150 (MQ=255) aGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCC‑ATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGacaaca > 2:892074/1‑150 (MQ=255) aaTCACCGACAATACCAGCAGAAATAGCAGGCTTTCACCAAATGAAAACCA‑ATTAATACAAACCTTCTTCCCAAAAC‑CC‑TGTCCCTCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATgcg < 2:476379/151‑1 (MQ=255) catcGAGCATCCCAGCAAAAAAAGCAGTCATCCAAAGAATGA‑AAGCCAATTAAAATCAATCTCCGTCCCAAAAC‑CC‑TGTGCCCCTCCCAAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATgcgc < 1:210214/149‑1 (MQ=255) atcGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCC‑ATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGcc > 1:747474/1‑149 (MQ=255) gCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCC‑ATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAAt > 2:322904/1‑151 (MQ=255) gCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCC‑ATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAAt > 1:82886/1‑151 (MQ=255) aCCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCC‑ATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATccc > 2:222522/1‑150 (MQ=255) aCCAGCAGAAATAGCAGGCATTCATCGAATAATAACCC‑TTTAATATAAACCTTCTTCCCAAAAC‑CC‑TGTCCCCCTCCCAAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATc < 1:669437/148‑1 (MQ=255) agaagaaaTAGCAGGCATTCAACGAATGAAAACCA‑ATTGATATCAAACTTCGTACCAAAAC‑CC‑TGTCCCACTCCCAAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGaaa < 1:366279/136‑1 (MQ=255) gcagcaATAGCAGGCATTCATCGAATGATAAGCC‑ATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCcagc > 1:625988/1‑149 (MQ=255) cagcaATAGCAGGCATTCATCGAATGATAAGCC‑ATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGt > 2:296074/1‑151 (MQ=255) agaaaTAGCAGGCATTCATCAAATAAAAAGCC‑ATTAAAAACAAACTTCTTTCAAAAAC‑CC‑TGTCCCACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCcag < 2:353712/143‑1 (MQ=255) acaaaTAGAGGGAATTCACCAAATAAAAACCC‑ATTAATATAAACCTTCTTCCCAAAAC‑CC‑TGTCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAg < 2:893021/133‑1 (MQ=255) aaTAGCAGGCTTTCATAAAATAAAAAGCC‑ATTAATACCAACCTTCGTCCCAAAAC‑CC‑TGTCCCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGt < 2:647280/147‑1 (MQ=255) aaGCAGGCTTTCATAAAATAAAAGGCT‑ATTAATA‑CAAACCTTCTTCCAAAAACCC‑TGGCCCCCTCCCGAAGCCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGcca < 2:723739/150‑1 (MQ=255) agcagcCATTCATGAAATAAAAACCC‑TTTGATATCAACCTCCTTCCCAAAAC‑CC‑TGTCCCACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGcca < 1:479258/150‑1 (MQ=255) aGGCATTCATCGAATGATAAGCC‑ATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCAt > 1:116311/1‑151 (MQ=255) aTTCATCGAATGATAAGCC‑ATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAATTTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATAc > 2:332246/1‑149 (MQ=255) aaTGATAACCC‑ATTGATATAAACCTTCGTCCAAAAAC‑CC‑TGTCCCCCTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAAc < 2:658208/147‑1 (MQ=255) aaTGAAAACCC‑ATTGATAAAAACCTTCTTCCCAAAAC‑CC‑TGTCCCCCTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACgg < 2:790636/149‑1 (MQ=255) aTGATAAGCC‑ATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGCTGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACgg > 2:831675/1‑148 (MQ=255) tAAGCC‑ATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAgg > 2:674740/1‑151 (MQ=255) catttaatacaaaaCTCCGTCCCAAAAC‑CC‑TGCCCCCCCCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGa < 1:267747/140‑1 (MQ=255) cAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGt > 2:537134/1‑151 (MQ=255) aaaCTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTc > 2:27683/1‑151 (MQ=255) aaaCTCCGTCCAAAAAC‑CC‑TGCCCCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTc < 1:114014/151‑1 (MQ=255) ttCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCAc > 1:544453/1‑151 (MQ=255) gTTCCGAACC‑CC‑TGTGCCCCTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCa > 2:274473/1‑150 (MQ=255) cGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATc > 1:780315/1‑149 (MQ=255) gAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGt > 2:231469/1‑151 (MQ=255) gAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGt > 2:123703/1‑151 (MQ=255) aaac‑CC‑TGTCCCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGtt < 2:210352/151‑1 (MQ=255) aaaa‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGtt > 2:142648/1‑151 (MQ=255) cc‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTg > 2:167072/1‑151 (MQ=255) cc‑TGTCCCACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTc < 2:522911/149‑1 (MQ=255) ccACTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGt < 2:900483/146‑1 (MQ=255) ttCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATg > 1:297155/1‑151 (MQ=255) | TATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCC‑ATTGATATCAAACTTCGTTCCGAAAA‑CC‑TGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATG > NC_000913/3601026‑3601325 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |