Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 648,716 | A→T | V61E (GTG→GAG) | mrdB ← | cell wall shape‑determining protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 648,716 | 0 | A | T | 97.5% | 150.1 / ‑3.2 | 40 | V61E (GTG→GAG) | mrdB | cell wall shape‑determining protein |
Reads supporting (aligned to +/- strand): ref base A (1/0); new base T (22/17); total (23/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01 |
CAAAATAATACAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCG > REL606/648643‑648787 | cAAAATAATACAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCacctcc > 1:1565667/1‑76 (MQ=255) cAAAATAAGACAGAGGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCacctcc < 2:1288915/76‑1 (MQ=255) aaaaTAATACAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCAcctcca > 2:1023002/1‑76 (MQ=255) taataCAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCatgat > 1:97763/1‑76 (MQ=255) aaaaCAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCatgatg > 2:46641/4‑76 (MQ=255) taCAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGAc > 1:1709468/1‑76 (MQ=255) aaCAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGAc < 2:1067344/75‑1 (MQ=255) cAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCAACTCCATGATGACCa < 1:988909/76‑1 (MQ=255) tgatATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGAcc < 1:1363458/76‑1 (MQ=255) atatAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGAACCa > 2:358561/1‑76 (MQ=255) taGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATcg < 2:183474/76‑1 (MQ=255) taGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATcg > 1:341581/1‑76 (MQ=255) agagATAGGGGGCTCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGAGCCATcgc < 1:1179806/76‑1 (MQ=255) gaTAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGAt > 1:827286/1‑76 (MQ=255) aTAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGAtt > 1:184587/1‑76 (MQ=255) aTAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGAtt > 1:87513/1‑76 (MQ=255) aTAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGAtt > 1:900476/1‑76 (MQ=255) ggggCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATAGCGATTTGGc > 2:1264606/1‑76 (MQ=255) ggCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCg < 1:1023002/76‑1 (MQ=255) cccAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGAt > 1:288905/1‑76 (MQ=255) ccAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGAtt > 2:1444849/1‑76 (MQ=255) ccAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGAtt > 1:1549559/1‑76 (MQ=255) aGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGAtttt < 2:226305/76‑1 (MQ=255) aGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGAtttt > 2:1316528/1‑76 (MQ=255) gCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTa > 1:188304/1‑76 (MQ=255) cTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACg < 2:520929/76‑1 (MQ=255) ttCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGc > 1:1400113/1‑76 (MQ=255) tCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCt > 2:1149667/1‑76 (MQ=255) aTAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTcc > 1:240817/1‑76 (MQ=255) aTAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTcc < 1:1776870/76‑1 (MQ=255) aTAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTcc > 2:1462550/1‑76 (MQ=255) aaaCGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCcat > 2:1453386/1‑76 (MQ=255) gTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGcc > 2:1396966/1‑76 (MQ=255) ggaATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAatat < 1:413914/76‑1 (MQ=255) ggaATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAatat < 1:75544/76‑1 (MQ=255) ggaATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAatat < 1:617883/76‑1 (MQ=255) gaATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATc < 1:618253/76‑1 (MQ=255) tGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGAc < 2:1311868/76‑1 (MQ=255) tGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGAc > 2:251739/1‑76 (MQ=255) ccATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCt < 1:774554/76‑1 (MQ=255) ccATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCt > 2:352981/1‑76 (MQ=255) aTCACGTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTgg < 2:1806274/76‑1 (MQ=255) tCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGgc < 2:123178/76‑1 (MQ=255) cacctcCATGATGNCCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGgcg > 1:1241199/1‑76 (MQ=255) cacctcCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGgcg > 2:397041/1‑76 (MQ=255) cacctcCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGgcg > 2:865337/1‑76 (MQ=255) | CAAAATAATACAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCG > REL606/648643‑648787 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |