Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,350,853 | 0 | A | C | 41.1% | 86.6 / 14.7 | 56 | S725A (TCA→GCA) | 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 (26/7); new base C (0/23); total (26/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.85e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTC > NC_000913/3350707‑3351003 | tGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑ccc < 1:257477/150‑1 (MQ=255) cccccGCTTTCAGCCCTTCTCCTTCGTGCCGCCCCTCTTCTTCCATCTCTTCAATCCATTCCCCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGTTGTTGCAGATGGGGTACCCCCACTGC‑CCccgccg < 1:401752/150‑1 (MQ=255) acGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCAcc > 1:112069/1‑151 (MQ=255) gCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCtt > 1:679941/1‑151 (MQ=255) tttCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCttt > 1:445945/1‑150 (MQ=255) cAGCTCTTCTACGTCGTGACCCCCCTTTTCTTTCATCTCTTCAATCCATCCCCCGACGTTATCTCCCCAGGCGGGAAGCTC‑AGGAGACTGAATTTG‑CTGCCCCATCTGTTGGAGAGGGCGTACCCCCACTGC‑CCCCGCCGCACCTTTa < 1:236931/148‑1 (MQ=255) accccTTCTCCGCCGGGCCCCCCCTCTTCTTTCATCTCTTCACCCCATCCCCCGACTTTATCTCCCCAGGCCGAAAGCTC‑AGGAGACTGAATTTC‑CTGCCCCGGCTTTTGAAGAGGCCGAACCCCCACTGA‑CCCCGCCGCACCTTTAAtt < 1:165474/146‑1 (MQ=255) cTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATg > 2:58101/1‑151 (MQ=255) ccttcgtGACCCCCCTCTTTTTTCTTCTCTCCAACCCATCCCCCGACTTTATCTCCCCGGGCCGGAAGGCC‑AGGAGACTAATTTCC‑CTCCCCCATTTTTTGAAGATGGCATACCCCCACTGC‑CCCCGCCGCACCTTTAATTTTATGtcctt < 1:459736/148‑1 (MQ=255) gtcgtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGtccttcc > 2:1007233/1‑151 (MQ=255) tGACCCCACTCTTTTTTCATCTCTCAACCCCTTTCCCCGACTTTATCTCCCCAGGCCGGAAGGCC‑GGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTACCCCCCCTGC‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTca < 2:1249459/150‑1 (MQ=255) cccccctctCCTTTCACCCCTCCAACCCATCCCCCAACTTTTTCTCCCCGGGCGGGAAGGTC‑GGGAG‑CCGAATTTGCCTCCCCCAGTTGTTGCAGATGGCGTACCCCCCCTTC‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaacaa < 1:592971/146‑1 (MQ=255) ccACTCTTCTTTCATCTCTTAAACCCTTCCACCGACTTTATCTCCCCAGGCGGGAAGGTC‑AGGAGACTGAATTTG‑CTGCCCCAGCTGTTGAAGATGGCGTACCCCCACGGC‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaacaa < 2:553637/149‑1 (MQ=255) cACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGc > 1:904091/1‑151 (MQ=255) acaattcttcCACCTCTTAACCCCTTCCCCCGCCGTTCTCTCCCCAGCCGGGGAGGGC‑AGGAGACTAAATTTC‑CAGACCTTGCTTTTGAAGATGGCATATCCCCGCGAA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGc < 1:377740/146‑1 (MQ=255) aCTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGcc > 2:471275/1‑151 (MQ=255) ttcttTCATCTCTAAAATCCATCCCCCGACTTTATCTCCCCAGGCGGGAAGGTC‑AGGAGACTGATTTGG‑CTGCCCCAGCTGTTGAAGAGGGCGTACCCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCt < 2:809435/148‑1 (MQ=255) tttttCATCTCTTAAACCCATCCCCCGACTTTTTCTCCCCAGGCCGGAGGTTC‑GGGGGACTGAATTTG‑CTGCCCCAGCTTTTGCAGATGGCGTACCCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttt < 2:1130511/148‑1 (MQ=255) tttCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt > 1:85541/1‑150 (MQ=255) tttCATCTCTTAAACCCATCCCCCGCCGTTATCTCCCCAGCCGGGAAGGTC‑AGGAGACTGATTTTG‑CTGCCCCAGCTGTTGAAGAGGGCGTACCCCCACTGC‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttt < 2:723542/148‑1 (MQ=255) ttccATCTCTTAAACCCATCCCCCGGCTTTATC‑TCCCAGGGCGGAAGGTCAAGGAGACTAATTTTG‑CTGCCCCAGCTTTTGAAGGGGGCTTACCCCCACCGC‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaaca < 1:1134686/136‑1 (MQ=255) tacttctcTTCATTCCCTCCCCCGAGGTTTCCTCCCCAGGCGGAAGGGTG‑AGGTAACTGATTTGC‑CTCCCCCATTTTTTGAAGATGCCGTCCCCCCCCTGC‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt < 1:1186588/145‑1 (MQ=255) aTCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCt > 2:659787/1‑151 (MQ=255) tctcTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTg > 1:125294/1‑151 (MQ=255) aTCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc > 2:260030/1‑150 (MQ=255) cATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc > 1:509648/1‑147 (MQ=255) tcccccGACTTTCTCTCCCCGGGCGGGAAGGTG‑AGGAAACTGATTTTG‑CTGCCCTTGTTTTTGGAATTGGCGAACCCCCACGCC‑ACCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAtt < 1:1109129/146‑1 (MQ=255) ccccgccgTTATCTCCCCGGGCGGAAAGGTC‑AGGAGACTGATTTTG‑CTGCCCCAGCTGTTGAAGATGGCGTACCCCCTCTGC‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAtt < 2:152295/146‑1 (MQ=255) ccccgattttctctcCCCAGGCGGAAAGTTG‑GGGAGACTGATTTGG‑CTCCCCCTCCTGTTGAAGTTGGCCAAACCCCATTGC‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGa < 2:445956/140‑1 (MQ=255) cGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTc > 1:9763/1‑151 (MQ=255) gACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTcc > 2:457878/1‑151 (MQ=255) tttctcCCCAGCCCGGAGGTTC‑GGGACATAAATTTTG‑CTGCCCAAGTTTTTGAAGAGGCCATACCCCCTCCGC‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCa < 1:642041/148‑1 (MQ=255) tATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCAc > 2:464284/1‑151 (MQ=255) cTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCt > 1:336619/1‑151 (MQ=255) ttCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTg > 2:1246488/1‑151 (MQ=255) ccccggcccGGAAGGAA‑GGGAGACATAGTTTG‑CTCACCAATGTTTGGCAGAGGGCGAACCCCCCCTGC‑CGCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCt < 1:339571/142‑1 (MQ=255) cccGGGCCGGAGGGTC‑AGGAGACTAATTTTG‑CTGACCCAGTTTTTGCAGGGGGCATACCCCCACTGC‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTg < 1:599365/149‑1 (MQ=255) cGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGccc > 1:819561/1‑151 (MQ=255) ggAAGGTC‑AGGACACTAAATTTG‑CTCCCCCAGTTGTTGCAGATGGCATACCCCCACGCC‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGccc < 2:426611/150‑1 (MQ=255) ggTC‑AGGAGACTAATTTTT‑CTGCCCCTGCTGTTGCAGATGGCTTACCCCCCCTGC‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGc < 1:806567/149‑1 (MQ=255) aggaaaGGAGACTGATTTTG‑CTCCCCCTGCTGTTGCAGGGGGCTAACCCCCCCTCC‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATa < 1:1041651/138‑1 (MQ=255) c‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatcat > 1:1212587/1‑151 (MQ=255) ggAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtct > 2:1033189/1‑151 (MQ=255) gagaCTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtctt > 2:281104/1‑151 (MQ=255) cTGATTTTG‑CTGCCCCAGTTGTTGCAGATGGCGTACCCCCACTCC‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttct < 2:784011/149‑1 (MQ=255) aTTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAac > 1:1253500/1‑150 (MQ=255) ttttgaagtggccGTACCCCAACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTc < 1:693015/137‑1 (MQ=255) ggggCATACCCCCACGCC‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTg < 1:660729/143‑1 (MQ=255) tGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCt > 2:178234/1‑151 (MQ=255) tGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCt > 1:923262/1‑151 (MQ=255) ggCTTACCCCCACTCC‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCa < 2:944000/147‑1 (MQ=255) acaaCCCCCCGGC‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCg < 1:254197/148‑1 (MQ=255) ttcacccACTAACCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCGCGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGc < 2:1032594/142‑1 (MQ=255) accccaaatca‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGa < 2:1139471/129‑1 (MQ=255) cccccACTCC‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTc < 1:683227/149‑1 (MQ=255) aCCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTcc > 2:134701/1‑151 (MQ=255) aCTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACa > 1:592031/1‑150 (MQ=255) a‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTc > 2:1183223/1‑151 (MQ=255) | TGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGA‑CCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTC > NC_000913/3350707‑3351003 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |