Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 73,581 | 0 | A | C | 40.4% | 97.0 / 16.4 | 57 | V314G (GTA→GGA) | thiP | thiamine/thiamine pyrophosphate ABC transporter permease |
Reads supporting (aligned to +/- strand): ref base A (29/5); new base C (0/23); total (29/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.85e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTC > NC_000913/73433‑73730 | cTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAcc > 2:1286268/1‑151 (MQ=255) gCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAc < 2:1224944/148‑1 (MQ=255) cAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATc < 1:858179/151‑1 (MQ=255) cgcgACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACCCCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGc < 1:521997/151‑1 (MQ=255) cgACTGCTCCATAGCAGCATCATGGCCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGcc > 1:519648/1‑150 (MQ=255) cTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACg > 2:475813/1‑151 (MQ=255) ccaccagggTCAGCCCTCCCCAAAAACCCCCGGCCGCCGCCGAAACCCGCAAGGGGGTCCACAGGCCGGGCCACAGCCCCGTTTGGGCCAGCACTCCCGGCAACTGGCGATTTACCCCGGCGACGATC‑ACCGCCAGTAACGGTGGCAGCaa < 1:710273/144‑1 (MQ=255) catGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAAcagc > 1:106082/1‑151 (MQ=255) ggTCAGCACTCCGCACAATCCCCCTCCCGCCAGCCAAATACGCAACGGGTTCCACAGCGCCTGCCACAGCACCGGTTGTCCCAGCACTTCCGGCAACGGGCGATTTCCCCCATCGACGATC‑ACCGCCAGTAACGGTGgca < 2:911388/140‑1 (MQ=255) ggTCAGCACTCCCCACAATCCCCCTGCCCCCAGCGCAATACGCAACGGGCCCCACAGCCCCTCCCACAGCCCCGGTTGGGCCAGCCCTTCCGGCACCGGGCGTTTTCCCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAAca < 2:1330871/146‑1 (MQ=255) cTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGcc > 2:422118/1‑149 (MQ=255) aCGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGc > 2:12189/1‑150 (MQ=255) cGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGca > 1:936252/1‑150 (MQ=255) cGCAAAATCCCCCTGCCCCCAGCGCAATACGCAACGAGGTCCACGGCGCCTGCCCCAGCCCCGTTTGTGCCACCACTGCCGGCAACTGGCGATTTCCCCCATCGACAATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCa < 2:1202728/147‑1 (MQ=255) aatcccccggcccccAGCGCAATACGCACCGGGGCCCACACCGCCTGCCACAGCCCCGTTTGTGCCAGCCCTCCCGGAAACGGGCGTTTTCCCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTa < 2:347165/142‑1 (MQ=255) aatccccccgccgccAGCACAATACGCAACGAGGCCCACAGCGCCTGCCACAGCCCCGGTTGTGCCACCACTCCCGGCAACTGGCGTTTTCCCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACaa < 2:291685/139‑1 (MQ=255) aatcccccTGCCGCCAGCGCAATACGAAACGGGGCCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTCCCGGCAACGGGCGTTTTCCCCCACCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTa < 2:1176922/145‑1 (MQ=255) accccccgcccgccAGCCCAATCCCCAGCGAGCTCCCCAGCGCCCGCCACAGCCCCGGTTGGCCCAGCCCTTCCGGCAACGGGCGATTTCCCCCATGGACGATA‑ACCGCCAGTAACGGTGGCAGCAACagcagcagc < 1:979483/128‑1 (MQ=255) aTACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAAtt > 1:1287142/1‑149 (MQ=255) acacCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAac > 2:716869/1‑150 (MQ=255) ccccTGCCCCCAGCGCAATACGCAACGGGGCCCACAGCCCCTGCCACAGCACCGGTTGGGCCAGCCCTTCCGGCAACGGGCGATTTCCCCCACCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAac < 2:123222/147‑1 (MQ=255) acCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACAcc > 1:2857/1‑151 (MQ=255) acCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACAcc > 2:406304/1‑151 (MQ=255) acCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACAcc > 2:1180404/1‑151 (MQ=255) ccTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCg > 1:936082/1‑151 (MQ=255) ccTGCCCCCAGCGCAATACGAAACGAGCCCCACAGCGCCTCCCACAGCCCCGTTTGCGCAAGCACTCCCGGCAACTGGCGATTTCCCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTaa < 2:206109/146‑1 (MQ=255) gccgccAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCgtgt > 1:1102507/1‑151 (MQ=255) cgccAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCg > 2:990134/1‑151 (MQ=255) cgccAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCg > 1:780751/1‑151 (MQ=255) aaTCCGCAAGGGGGCCCCCGGCGCCTGCCACACCCCCGTTTGGGCCGGCCCTTCGGGAAACTGGCGATTTCCCCCACCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCg < 2:357948/150‑1 (MQ=255) aTACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGc > 1:225785/1‑151 (MQ=255) tacccacaagcccccACCGCTCCTGCGACGGACCCGGTTGTGCCGGATGTTCGGGGACGTGAGGATTGCCCCCACCGAC‑ATCAACCGCCAGTATCGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGAAAATGc < 1:475851/134‑1 (MQ=255) gCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATg > 1:847112/1‑151 (MQ=255) caCGAGGCCCACAGCCCCTGACACAGCACCGGTTGTGCCAGCTCTTTAGGGAACTGGCGATTTACCCCATCGAGGATC‑ACCGCCAGTAACGCTATGAGCAACAGCAGCAGCTGCAGCACAATTAAGACCGTGACCCAAACGCGGCGATg > 1:774278/2‑149 (MQ=255) gAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGAc > 2:990022/1‑151 (MQ=255) gAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGAc > 2:1323334/1‑151 (MQ=255) gcccaacaGCGCCTGCCACAGCACCGGTTGTCCCAGCACTCCCGCCATCGGGCGATTTCCCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGa < 2:797639/146‑1 (MQ=255) acaGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTc > 2:774261/1‑151 (MQ=255) acaGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTc > 1:132749/1‑151 (MQ=255) aGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCg > 1:653199/1‑151 (MQ=255) ccggCCACGGCTCCGGTGGTGCAACCCCTGCCGGCAACGGGGGATTTCCCCCATCGACAATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGt < 2:1229680/147‑1 (MQ=255) ccTCCCACAGCCCCGGTTGTGCCAGCCCTCCCGGCAACGGGCGTTTTCCCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTcc < 2:1286479/146‑1 (MQ=255) tGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTcgc > 1:822716/1‑151 (MQ=255) ccACAGCCCCGTTTGTGCCAGCACTCCCGGCAACGGGCAATTTCCCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTcgcg < 2:152565/150‑1 (MQ=255) ccACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCgcgc > 1:1267920/1‑151 (MQ=255) cacaGCCCGGGTTGTGCAAGCACTCCCGGCAACGGGCGATTTCCCCCACCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCt < 1:560019/125‑1 (MQ=255) ccGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTg > 2:318354/1‑151 (MQ=255) ccGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCAACGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTg > 2:122896/1‑151 (MQ=255) cacgatctccgggCACCTGGCTGTTTCCCCCCTAGCCGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTgcagc < 1:104575/135‑1 (MQ=255) caccACTCCCGGCAACTGGCATTTTCCCCCACCGACGATC‑CCCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCgtg < 2:533863/148‑1 (MQ=255) ccctcccGGCAACGGGCGATTTCCCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCgtg < 2:1165357/143‑1 (MQ=255) ttCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCg > 1:645596/1‑151 (MQ=255) cgggCACCGGGCGATTCCCCCCTGCGACCACC‑CCCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAAGTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGcc < 1:1058631/146‑1 (MQ=255) ccGGAAACTGGCTTTTTCCCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCg < 2:206996/149‑1 (MQ=255) ggCGATTTCCCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGc < 2:578939/151‑1 (MQ=255) ggCGATTTACCCCACCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATg < 2:1319596/149‑1 (MQ=255) cGATTTCCCCCACCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGcc < 1:1232603/150‑1 (MQ=255) aCCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTc > 1:1319224/1‑150 (MQ=255) | CTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATC‑ACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTC > NC_000913/73433‑73730 |
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: ‑ |