Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,350,844 | 0 | A | C | 47.4% | 51.7 / 12.4 | 38 | L728V (TTA→GTA) | arcB | aerobic respiration control sensor histidine protein kinase, cognate to two‑component response regulators ArcA and RssB |
Reads supporting (aligned to +/- strand): ref base A (16/4); new base C (0/18); total (16/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.25e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGC‑GTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGC > NC_000913/3350695‑3350986 | tttAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGC‑GTAAc > 1:293278/1‑151 (MQ=255) tAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGC‑GTAAccc > 1:102724/1‑151 (MQ=255) ttttGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGC‑GTAACCCCACTGAc > 1:369766/1‑151 (MQ=255) tGCCTCCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGC‑GTAACCCCACTGAccc < 2:315583/150‑1 (MQ=255) cccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGC‑GTAACCCCACTGACCCcgccgc > 1:312318/1‑151 (MQ=255) cccacCCTTTCAGCCCTTCTACGTCGTGACCCCACTCTTCTTTCATCTCTTCAATCCATCCACCGACGTTATCTCCCCAGGCCGGAAGTCCGGGAGCCTGATTTTGCGGCCCCAGCTGTTGCAGATGGC‑GTACCCCCACTGCCCCcgccgc < 2:200650/151‑1 (MQ=255) cGCTTTCAGCACTTCTACGTCGTGCCGCCCCTCTTCTTTCATCTCTCCAACCCATCCCCCGACGTTATCTCCCCAGGCCGGAAGGTCAGGAGACTAATTTTGCTGCCCCAGCTGTTGCAGAGGGC‑GTACCCCCACTGACCCCGCCGCAcc < 2:672033/150‑1 (MQ=255) gCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGC‑GTAACCCCACTGACCCCGCCGCACCtt > 2:408889/1‑151 (MQ=255) ttCACCATTTCTCCGTCGGGCCGCCCCTTTTCTTTCATCTCTACAACCCATCCCCCGTCTTTTTCTCCCCAGGCCGGAAGGTCAGGAGACTAATTTTGCTCACCCAGCTGTTGCAGAGGGC‑GTAACCCCCCTGCCCCCGCCGCACCtt < 2:239653/148‑1 (MQ=255) cACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGC‑GTAACCCCACTGACCCCGCCGCACCTTTAAtttt > 2:484845/1‑150 (MQ=255) ttCTCCGTCGTGACCCCCCTCTTTTTCCATCTCTCCAACCCATCCCCCGACTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGATTTTGCGGCCCCAGCTTTTGCAGATGGC‑GAACCCCCCCTGACCCCGCCGCACCTTTAAttt < 2:291613/146‑1 (MQ=255) ttCTACGTCGTGACCCCCCTCTTCTTTCATCTCTCAAACCCATCCCCCGACGTTATCTCCCCAGGCGGAAGGGTCAGGAGACTAAATTTGCTGCCCCAGCTGTTGAAGAGGGC‑GTACCCCCACGGACCCCGCCGCACCTTTAATTTTAt < 1:407570/149‑1 (MQ=255) cgtcgtGACGCCCCTCTTCTTTCATCTCTTAAACCCATCCCCCGACGTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGATTTTGCTCCCCCAGTTGTTGCAGAGGGC‑GTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGtcct < 2:455/149‑1 (MQ=255) cgtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGC‑GTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTc > 1:429734/1‑151 (MQ=255) gCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGC‑GTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAAt > 2:264854/1‑151 (MQ=255) cccctctTCTTTCATCTCTCCAACCCATCCACCGACTTTATCTCCCCAGGCGGGAAGGGCAGGAGACTGAATTTGCTGCCCCAGTTGTTGCAGATGGC‑GTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaacaa < 1:516678/146‑1 (MQ=255) ccctctTCTTTCATCTCTCCAACCCTTCCCCCGCCGTTTTCTCCCCAGGCGGAAAGGTCAGGAGACTGATTTTGCTGCCCCAGTTTTTGCAGATGGC‑GTACCCCCACGCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGc < 2:395154/149‑1 (MQ=255) ttcttTCATCTCTACAATCATTCCCCCCCCTTTTTCTCCCAAGGCGGGAAGATCAGAAGAGTAATTTTCCCCCCCCATTTTTTGTA‑ATGGCTGCACCCCCCTTTCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtt < 2:267040/149‑1 (MQ=255) ttacttctttccAATCCATCCCCCGGCTTTATCTCCCCGGGCGGGAGGGTCGGGAGTCTAAATTTGCTGCCCCAGTTTTTACAGATGGT‑GTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtt < 1:366781/135‑1 (MQ=255) ttCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGC‑GTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgcgc > 1:519057/1‑150 (MQ=255) ccATCCCCCGACGTTATCTCCCCGGCCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGTTGTTGCAGATGGC‑GTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGa < 1:675040/151‑1 (MQ=255) cATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGC‑GTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc > 1:147556/1‑147 (MQ=255) aTTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGC‑GTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAtt > 2:678009/1‑151 (MQ=255) ccccGACGTTTTCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGC‑GTACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTc < 1:59716/147‑1 (MQ=255) ccgcctttATCTCCCCGGCCCGGAAGGTCAGGAGACTGATTTTGCTCCCCCAGCTTTTGCAGATGGC‑GTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAt < 2:132734/139‑1 (MQ=255) actttATCTCCCCGGGCCGGAGGTTCAGGAGACTGATTTTCCTGACCCAGTTTTTGCAGATGGC‑CAAACCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTcc < 2:658949/147‑1 (MQ=255) ttATCTCCCCAGCCCGGAGGGTCAGGAGACGAAATTTGCTACCCCAGTTTTTGCAGATGGC‑GTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCa < 2:102723/148‑1 (MQ=255) tccccAGGCGGGAAGCACAGGAGTCTAATTTTTCTAACCCAGTTTTTGCAGATGGC‑GTAACCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCAc < 2:438540/145‑1 (MQ=255) ggCGGGAAGGCAAGGAGACTGATTTTGCTCACCCATTTGTTGCAGAGGGC‑GTACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCC‑TTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATa < 2:68167/149‑1 (MQ=255) cAGGAGCCTGATTTTGCTGCCCCAGCTTTTGCAGAGGGC‑TTACCCCCGCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCgg < 2:659567/145‑1 (MQ=255) cAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGC‑GTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatcat > 1:312662/1‑151 (MQ=255) agaCTGAATTTGCTGCCCCAGCTTTTGCAGATGGC‑TTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttc < 1:572632/151‑1 (MQ=255) gAATTTGCTGACCCAGCTGTTGCAGATGGC‑GTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCaa > 2:559421/1‑150 (MQ=255) ttttcaaaCCCAGTTGTTAGAGAGGGC‑ATACCCCCGCTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAaca < 2:89932/143‑1 (MQ=255) cccccacttGTTGAAGATGCC‑GAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTaa < 2:440800/143‑1 (MQ=255) aCCCAGCTGTTGCAGATGGC‑GTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAc > 1:231576/1‑151 (MQ=255) cTGTTGCAGATGGC‑GTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCg > 1:471696/1‑151 (MQ=255) gacaagggtgcaaCCCCTCCGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGt < 2:249471/133‑1 (MQ=255) aTGGC‑GTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGc > 1:427554/1‑151 (MQ=255) | TTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGC‑GTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGC > NC_000913/3350695‑3350986 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |