Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,350,844 | 0 | A | C | 44.4% | 84.6 / 11.9 | 55 | 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 (26/4); new base C (0/24); total (26/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.48e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCC > NC_000913/3350695‑3350994 | tttAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTaa < 1:942523/150‑1 (MQ=255) ttAGAGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAAcc > 2:884533/1‑151 (MQ=255) aGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAAccc < 1:1255596/150‑1 (MQ=255) gTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAAcccc > 2:683251/1‑150 (MQ=255) gTGGCTTTTGCCACCCACGCTTTCAGCACTTCCACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCa < 2:1075259/151‑1 (MQ=255) tGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAccc < 1:257477/150‑1 (MQ=255) cccccGCTTTCAGCCCTTCTCCTTCGTGCCGCCCCTCTTCTTCCATCTCTTCAATCCATTCCCCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGTTGTTGCAGATGGGGTACCCCCACTGCCCccgccg < 1:401752/150‑1 (MQ=255) acGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCAcc > 1:112069/1‑151 (MQ=255) gCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCtt > 1:679941/1‑151 (MQ=255) tttCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCttt > 1:445945/1‑150 (MQ=255) cAGCTCTTCTACGTCGTGACCCCCCTTTTCTTTCATCTCTTCAATCCATCCCCCGACGTTATCTCCCCAGGCGGGAAGCTC‑AGGAGACTGAATTTG‑CTGCCCCATCTGTTGGAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTa < 1:236931/148‑1 (MQ=255) accccTTCTCCGCCGGGCCCCCCCTCTTCTTTCATCTCTTCACCCCATCCCCCGACTTTATCTCCCCAGGCCGAAAGCTC‑AGGAGACTGAATTTC‑CTGCCCCGGCTTTTGAAGAGGCCGAACCCCCACTGACCCCGCCGCACCTTTAAtt < 1:165474/146‑1 (MQ=255) cTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATg > 2:58101/1‑151 (MQ=255) ccttcgtGACCCCCCTCTTTTTTCTTCTCTCCAACCCATCCCCCGACTTTATCTCCCCGGGCCGGAAGGCC‑AGGAGACTAATTTCC‑CTCCCCCATTTTTTGAAGATGGCATACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGtcctt < 1:459736/148‑1 (MQ=255) gtcgtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtccttcc > 2:1007233/1‑151 (MQ=255) tGACCCCACTCTTTTTTCATCTCTCAACCCCTTTCCCCGACTTTATCTCCCCAGGCCGGAAGGCC‑GGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTca < 2:1249459/150‑1 (MQ=255) cccccctctCCTTTCACCCCTCCAACCCATCCCCCAACTTTTTCTCCCCGGGCGGGAAGGTC‑GGGAG‑CCGAATTTGCCTCCCCCAGTTGTTGCAGATGGCGTACCCCCCCTTCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaacaa < 1:592971/146‑1 (MQ=255) ccACTCTTCTTTCATCTCTTAAACCCTTCCACCGACTTTATCTCCCCAGGCGGGAAGGTC‑AGGAGACTGAATTTG‑CTGCCCCAGCTGTTGAAGATGGCGTACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaacaa < 2:553637/149‑1 (MQ=255) cACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGc > 1:904091/1‑151 (MQ=255) acaattcttcCACCTCTTAACCCCTTCCCCCGCCGTTCTCTCCCCAGCCGGGGAGGGC‑AGGAGACTAAATTTC‑CAGACCTTGCTTTTGAAGATGGCATATCCCCGCGAACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGc < 1:377740/146‑1 (MQ=255) aCTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGcc > 2:471275/1‑151 (MQ=255) ttcttTCATCTCTAAAATCCATCCCCCGACTTTATCTCCCCAGGCGGGAAGGTC‑AGGAGACTGATTTGG‑CTGCCCCAGCTGTTGAAGAGGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCt < 2:809435/148‑1 (MQ=255) tttttCATCTCTTAAACCCATCCCCCGACTTTTTCTCCCCAGGCCGGAGGTTC‑GGGGGACTGAATTTG‑CTGCCCCAGCTTTTGCAGATGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttt < 2:1130511/148‑1 (MQ=255) tttCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt > 1:85541/1‑150 (MQ=255) tttCATCTCTTAAACCCATCCCCCGCCGTTATCTCCCCAGCCGGGAAGGTC‑AGGAGACTGATTTTG‑CTGCCCCAGCTGTTGAAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttt < 2:723542/148‑1 (MQ=255) ttccATCTCTTAAACCCATCCCCCGGCTTTATC‑TCCCAGGGCGGAAGGTCAAGGAGACTAATTTTG‑CTGCCCCAGCTTTTGAAGGGGGCTTACCCCCACCGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaaca < 1:1134686/136‑1 (MQ=255) tacttctcTTCATTCCCTCCCCCGAGGTTTCCTCCCCAGGCGGAAGGGTG‑AGGTAACTGATTTGC‑CTCCCCCATTTTTTGAAGATGCCGTCCCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt < 1:1186588/145‑1 (MQ=255) aTCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCt > 2:659787/1‑151 (MQ=255) tctcTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTg > 1:125294/1‑151 (MQ=255) aTCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc > 2:260030/1‑150 (MQ=255) cATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc > 1:509648/1‑147 (MQ=255) tcccccGACTTTCTCTCCCCGGGCGGGAAGGTG‑AGGAAACTGATTTTG‑CTGCCCTTGTTTTTGGAATTGGCGAACCCCCACGCCACCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAtt < 1:1109129/146‑1 (MQ=255) ccccgccgTTATCTCCCCGGGCGGAAAGGTC‑AGGAGACTGATTTTG‑CTGCCCCAGCTGTTGAAGATGGCGTACCCCCTCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAtt < 2:152295/146‑1 (MQ=255) ccccgattttctctcCCCAGGCGGAAAGTTG‑GGGAGACTGATTTGG‑CTCCCCCTCCTGTTGAAGTTGGCCAAACCCCATTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGa < 2:445956/140‑1 (MQ=255) cGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTc > 1:9763/1‑151 (MQ=255) gACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTcc > 2:457878/1‑151 (MQ=255) tttctcCCCAGCCCGGAGGTTC‑GGGACATAAATTTTG‑CTGCCCAAGTTTTTGAAGAGGCCATACCCCCTCCGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCa < 1:642041/148‑1 (MQ=255) tATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCAc > 2:464284/1‑151 (MQ=255) cTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCt > 1:336619/1‑151 (MQ=255) ttCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTg > 2:1246488/1‑151 (MQ=255) ccccggcccGGAAGGAA‑GGGAGACATAGTTTG‑CTCACCAATGTTTGGCAGAGGGCGAACCCCCCCTGCCGCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCt < 1:339571/142‑1 (MQ=255) cccGGGCCGGAGGGTC‑AGGAGACTAATTTTG‑CTGACCCAGTTTTTGCAGGGGGCATACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTg < 1:599365/149‑1 (MQ=255) cGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGccc > 1:819561/1‑151 (MQ=255) ggAAGGTC‑AGGACACTAAATTTG‑CTCCCCCAGTTGTTGCAGATGGCATACCCCCACGCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGccc < 2:426611/150‑1 (MQ=255) ggTC‑AGGAGACTAATTTTT‑CTGCCCCTGCTGTTGCAGATGGCTTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGc < 1:806567/149‑1 (MQ=255) aggaaaGGAGACTGATTTTG‑CTCCCCCTGCTGTTGCAGGGGGCTAACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATa < 1:1041651/138‑1 (MQ=255) c‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatcat > 1:1212587/1‑151 (MQ=255) ggAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtct > 2:1033189/1‑151 (MQ=255) gagaCTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtctt > 2:281104/1‑151 (MQ=255) cTGATTTTG‑CTGCCCCAGTTGTTGCAGATGGCGTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttct < 2:784011/149‑1 (MQ=255) aTTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAac > 1:1253500/1‑150 (MQ=255) ttttgaagtggccGTACCCCAACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTc < 1:693015/137‑1 (MQ=255) ggggCATACCCCCACGCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTg < 1:660729/143‑1 (MQ=255) tGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCt > 1:923262/1‑151 (MQ=255) tGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCt > 2:178234/1‑151 (MQ=255) ggCTTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCa < 2:944000/147‑1 (MQ=255) acaaCCCCCCGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCg < 1:254197/148‑1 (MQ=255) aCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTcc > 2:134701/1‑151 (MQ=255) | TTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTC‑AGGAGACTGAATTTG‑CTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCC > NC_000913/3350695‑3350994 |
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: ‑ |