| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 624,598 | 0 | A | T | 36.4% | 64.2 / 11.1 | 46 | intergenic (‑88/‑287) | fepB/entC | ferric enterobactin ABC transporter periplasmic binding protein/isochorismate synthase EntC |
| Reads supporting (aligned to +/- strand): ref base A (19/9); new base T (1/15); total (20/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.11e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAA‑GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATT > NC_000913/624458‑624737 | aGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAA‑GGGCGTAATGACAAAGTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGg‑aga < 2:412312/141‑1 (MQ=255) tGTTAATAGATGGGGGTTGCGGTAGTGGGGGGCGAGTCTCAAAAATCAGCTTCCTTTTATTAATATGGTTCT‑GGGCTTATTGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑AGTGGGTTTGGGTg < 2:102060/141‑1 (MQ=255) tGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAA‑GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGTg > 2:260871/1‑141 (MQ=255) aTAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAA‑GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAgggg > 1:685911/1‑141 (MQ=255) tAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAA‑GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGa > 2:230814/1‑141 (MQ=255) tAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAA‑GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGa > 2:943737/1‑141 (MQ=255) gggCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAA‑GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑AGTGGGTTAGGGTGAGGGGAAca‑gc > 1:525671/1‑141 (MQ=255) gggCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAA‑GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGAAca‑gc > 2:1129097/1‑141 (MQ=255) gggCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAA‑GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGAAca‑gc > 1:1956450/1‑141 (MQ=255) gggCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAA‑GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGAAca‑gc > 1:1898139/1‑141 (MQ=255) gtggggggCGAGTCACACAAATCAGCTTCCTTTTTTTAATAAGGTTAA‑GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑TGTGGGTTTGGGTGAGGGGAACA‑GCCAGCACTGGTGCg < 2:1944173/137‑1 (MQ=38) gagCGGGGCGAGTCTCACAAAACAGCTTCCTGTTATTAATAAGGTTAA‑GGGCGTAATGACAAATTCGACAAAGCGCACAACCCGTCCCCCCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCg < 2:1944381/141‑1 (MQ=255) gCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAA‑GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTa > 2:946980/1‑141 (MQ=255) tctcACAAATCAGCTTCCTGTTATTAATAAGGTTAA‑GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑AGA‑GGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTc < 1:1260230/141‑1 (MQ=255) ctcACAAATCAGCTTCCTGTTATTAATAAGGTTAA‑GGGCGTATTGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTc < 2:1105693/141‑1 (MQ=255) ctcACAAATCAGCTTCCTGTTATTAATAAGGTTAA‑GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCTCCTCGCCCCTTTGGGG‑TGTGGGTTTGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTc < 1:1970322/141‑1 (MQ=255) ctcAAAATTCAGCTTCCTTTTATTAATATGTTTAT‑GGGCGTAATGAAAATTTCCCAAAAGCGCACAATCCTTCCCCCCGCCCTTTTGGGG‑AGTGGGTTCGGGGGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTc < 1:350548/141‑1 (MQ=255) ccaacaAATCAGCTTCCTTTTATTAATATGGTTAC‑GGGCGTAATTAAAAATTCGAAAACGCGCACAAACCGTCCCCTCGCCCGTTGGGGG‑TGTGGGTTCGGGGGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTc < 1:945333/138‑1 (MQ=255) ataaaaaaTCAGATTCCTTTTATTAATATGTTTAT‑GGGCTTAATTAAAATTTCGACAAAGCGCACAATCCGTCCCATC‑CCCCTTTGGGGTTGTGGGTTTGGGTGTGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTc < 1:472195/136‑1 (MQ=255) cacaAATCAGCTTCCTGTTATTAATAAGGTTAA‑GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCAc > 1:288590/1‑141 (MQ=255) cacaAATCAGCTTCCTGTTATTAATAAGGTTAA‑GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCAc > 1:288905/1‑141 (MQ=255) cacaAATCAGCTTCCTGTTATTAATAAGGTTAA‑GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCAc > 1:46254/1‑141 (MQ=255) ctccgttttataaAT‑AGGTTAACGGGCTTATTCACAAATTCCAAAACGCCCACAAGCCGCCCCCGCGCCCGTGGGGGG‑TGTGGTTTCGGGTGCGGGCAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTc < 1:529365/134‑1 (MQ=255) ctcccTGTTATTAATACGTTTAC‑GGGCTTAATCAAAATTTCCACAACGCCCACAATCCGTCCCCCCCCCCTTTGGGGG‑TGTGGGTTTGGGTGTGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTc < 1:1481237/138‑1 (MQ=255) tcgtGTTATTAATAAGGTTAA‑GGGCGTAATTACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGT‑TGAGGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCAc < 2:222327/138‑1 (MQ=255) cTGATAGTAATAAGGTAAA‑GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCGCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCACACCCCAGCCCTCAccc < 2:1103589/141‑1 (MQ=37) gTTATTAATAAGGTTAA‑GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTg > 1:54649/1‑141 (MQ=255) tggTTAT‑GGGCTTAATGAAAATTTCTAAAATGCCCACAATCCGACCCCCCGCCCTTTTGGGT‑TGTGGTTTTGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGagag < 1:89221/140‑1 (MQ=255) tggTTAT‑GGGCGTAATCAAAATTTCCAAAACGCGCAAAAACCGTCCCCACGCCCCTTTGGGG‑AGGGGTTTTGGGTGTGGGCAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGagag < 1:1367410/140‑1 (MQ=255) cggTTAT‑GGGCTTAATGACAAATCCGACAACGCGCAAAATCCGCCCCCCCGCCCTTTTGGGT‑TGTGGGTTTGGGTGCGGGCAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGagag < 1:1762056/140‑1 (MQ=255) cggTTAC‑GGGCTTAATCACAAATCCCAAAACGCCAACAAGCCGTCCCCCCGCCCCTTTGGGT‑TGTGGTTTTGGGGGCGGGCAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGagag < 1:1762048/140‑1 (MQ=255) cgtttccgggCTTAATCAAAAATGCCACAACGCCCACAACCGGCCCCCCCGCCCGTTGGGGG‑TGTGGGTTTGGGGGCGGGCCACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGagag < 1:1781191/134‑1 (MQ=255) gTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg > 1:494959/1‑141 (MQ=255) taTTCGACAAAGCGCACATTCCGTCCCATCGCCCGTTTGGGT‑TGTGGGTTTGGGTTTGGGGAACA‑GCCAGCACTGGTGCGAAAATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACTGCGCAGGAc < 2:1618113/140‑1 (MQ=255) aaTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGAc < 2:1152563/141‑1 (MQ=255) ttCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTGGGGG‑AGTGGGTTTGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACAt < 1:943737/141‑1 (MQ=255) aaGCGCACAATCCGTCCCCACGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGCACCA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGACGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTg < 2:684539/141‑1 (MQ=255) cgcACAATCCGTCCCCTCGCCCCTTTGGGG‑AGTGGGTTCGGGTGAGGGCAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTgcgc < 2:292111/141‑1 (MQ=255) gTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAAt < 2:1956450/141‑1 (MQ=255) ccccTCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATcc > 1:1463077/1‑141 (MQ=255) cccTCGCCCCTTTGGGG‑TGAGGGTTAGGGTGTGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCa < 2:1483212/141‑1 (MQ=255) ccTCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCAt > 1:169738/1‑141 (MQ=255) tCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCa > 1:873668/1‑141 (MQ=255) tCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCa > 2:1945775/1‑141 (MQ=255) tttGGGG‑AGAGGGTTAGGGTGAGGGGAACAGGCCGGCACTGGCGCGAATATTTACCCTCACCCCGGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATg > 1:1031498/1‑141 (MQ=18) tttGGGG‑AGAGGGTTAGGGTGAGGGGAACAGGCCGGCACTGGCGCGAATATTTACCCTCACCCCGGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATg > 1:1031540/1‑141 (MQ=18) ttGGGG‑AGCGGGTTAGGGTGCGGGCAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCCCATTGCGCTTATGCGAATCCATCAATAATGCt < 2:322273/141‑1 (MQ=255) gagGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGACGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAtt > 2:409492/1‑141 (MQ=255) gagGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAtt > 2:409439/1‑141 (MQ=255) gagGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAtt > 2:409468/1‑141 (MQ=255) gagGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAtt > 1:1815022/1‑141 (MQ=255) gagGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAtt > 1:1814976/1‑141 (MQ=255) gagGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAtt > 1:1814687/1‑141 (MQ=255) | AGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAA‑GGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGG‑AGAGGGTTAGGGTGAGGGGAACA‑GCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATT > NC_000913/624458‑624737 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |