Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | CP000730 | 1,002,200 | T→G | M173R (ATG→AGG) | mecA2 → | competence negative regulator MecA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP000730 | 1,002,200 | 0 | T | G | 95.2% | 144.4 / ‑4.1 | 42 | M173R (ATG→AGG) | mecA2 | competence negative regulator MecA |
Reads supporting (aligned to +/- strand): ref base T (1/1); new base G (11/29); total (12/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.95e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AGGTTCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATATGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTATGCT > CP000730/1002063‑1002336 | aGGTTCTAAACGTCAAAAGTCTTCAGCACGTAACAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATTTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGt < 2:337178/141‑1 (MQ=255) aGGTTCTAAACGTCAAAAGTCTTCAGCACGTAACAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATCTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGt < 2:337181/141‑1 (MQ=255) tAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATgg < 1:4705/141‑1 (MQ=255) aaaaGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACttatt > 2:438921/1‑141 (MQ=255) cTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCACTCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATg > 1:208173/1‑141 (MQ=255) cAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTg > 2:239317/1‑141 (MQ=255) aGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATATGGTTGATGGTACTTATTATTATGCTGt > 1:380093/1‑141 (MQ=255) aGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGt < 1:281247/141‑1 (MQ=255) aGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGt < 1:281256/141‑1 (MQ=255) gTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATAtt < 1:405776/141‑1 (MQ=255) aaaTACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATCTTATGAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGa < 1:264024/141‑1 (MQ=255) aaTACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGAt > 1:72174/1‑141 (MQ=255) aTACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATa > 2:94271/1‑141 (MQ=255) aTACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATa > 2:205831/1‑141 (MQ=255) cAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTc > 1:46600/1‑141 (MQ=255) aaTCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATATGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGa < 2:380093/141‑1 (MQ=255) cATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCa > 1:187826/1‑141 (MQ=255) ttAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCattaatgacttc < 1:4946/141‑12 (MQ=255) gATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGAt < 1:387939/141‑1 (MQ=255) aagaTGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACa < 1:80679/141‑1 (MQ=255) aagaTGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACa < 2:272067/141‑1 (MQ=255) agaTGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTACTGATAGTTACAg < 2:16092/141‑1 (MQ=255) agaTGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAg < 2:204354/141‑1 (MQ=255) tattaATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAAtt < 2:226936/141‑1 (MQ=255) taattaTGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCt < 1:205831/141‑1 (MQ=255) attaTGCATATCATAGCAATACATTTACTCTAGTTGTTGAATTTTTTTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTg < 1:239317/141‑1 (MQ=255) gCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCaa < 1:215888/141‑1 (MQ=255) aTCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAacag < 2:46600/141‑1 (MQ=255) ccAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAacagac < 2:168838/141‑1 (MQ=255) ccAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAacagac < 2:291831/141‑1 (MQ=255) ccAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAacagac < 2:161747/141‑1 (MQ=255) ccAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAacagac < 1:145018/141‑1 (MQ=255) cAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATATTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAacagaca > 1:221223/1‑141 (MQ=255) aaTAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAacagacag < 2:252601/141‑1 (MQ=255) aCTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGacagaaca > 1:149407/1‑141 (MQ=255) cTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGacagaacag < 1:349356/141‑1 (MQ=255) tACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGacagaacaga < 1:419/141‑1 (MQ=255) gagTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGtt < 2:155238/141‑1 (MQ=255) agTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGttt < 1:136305/141‑1 (MQ=255) tGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTc > 1:187089/1‑76 (MQ=255) tGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTc < 2:187089/76‑1 (MQ=255) gAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGtttattta < 1:246171/141‑1 (MQ=255) tatagGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTATGCt > 1:394266/1‑141 (MQ=255) | AGGTTCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATATGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTATGCT > CP000730/1002063‑1002336 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |