Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,530,076 | (A)7→8 | coding (155/483 nt) | ydcD → | putative immunity protein for RhsE |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,530,069 | 1 | . | A | 95.0% | 204.9 / 3.9 | 60 | coding (148/483 nt) | ydcD | putative immunity protein for RhsE |
Reads supporting (aligned to +/- strand): ref base . (1/2); new base A (27/30); total (28/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01 |
TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGC > NC_000913/1529970‑1530162 | tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑g < 2:195696/100‑1 (MQ=255) tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑g < 1:895278/100‑1 (MQ=255) tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑g > 1:1570005/1‑100 (MQ=255) gTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑a > 1:109254/1‑100 (MQ=255) gTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑a > 2:436167/1‑100 (MQ=255) tttGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑aaaaaa < 1:670324/100‑1 (MQ=255) tgttgtTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑aaaaaaaa > 2:551630/1‑99 (MQ=255) gttgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑aaaaaaaat < 2:1214670/100‑3 (MQ=255) tgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑aaaaaaaatat > 1:642053/1‑96 (MQ=255) ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATAttattta > 1:1447697/1‑100 (MQ=255) ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATAttattta < 1:599202/100‑1 (MQ=255) cTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATAttattta > 2:571954/1‑99 (MQ=255) ttGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTtata > 2:1266457/1‑100 (MQ=255) gTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAaat > 2:1178070/1‑100 (MQ=255) gATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAAAAAAAA‑AAAAAAAAATATTATTTATAAaataat < 2:109254/99‑1 (MQ=255) gATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACAAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAAtaatt < 2:657005/100‑1 (MQ=255) tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAAtaattaa > 2:275575/1‑100 (MQ=255) tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAAtaattaa > 1:1393542/1‑100 (MQ=255) aaCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAAtaattaat > 2:1454916/1‑100 (MQ=255) aCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATg > 1:337504/1‑100 (MQ=255) aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGAcc < 1:1286379/100‑1 (MQ=255) ttCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGAcc > 2:1176911/1‑99 (MQ=255) aaaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCaaaa < 1:1170897/100‑1 (MQ=255) aTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCaaaaaa < 1:948447/100‑1 (MQ=255) cTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATta > 1:277106/1‑100 (MQ=255) cTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATta < 1:1034004/100‑1 (MQ=255) aCATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATAAAAAAAAAA‑AAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagta < 1:649740/100‑1 (MQ=255) aCAACGTATATTATGGAAATAAAAAACAAAAGAGGGTTGAGATATAAATTTCA‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagt < 2:448403/99‑1 (MQ=255) aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagtag < 2:1272649/99‑1 (MQ=255) aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga > 1:549300/1‑100 (MQ=255) tCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag < 2:642051/100‑1 (MQ=255) cGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAgg > 1:1316875/1‑100 (MQ=255) gTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAggg < 1:976689/100‑1 (MQ=255) tATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGa < 1:436168/100‑1 (MQ=255) aTTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAAt > 1:1539552/1‑100 (MQ=255) ttttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAAt < 1:381531/100‑1 (MQ=255) tttATGGACATAAAACACCAATGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTCTAAAATGATTAATGACCAACAAACTAGTatagtgcat > 2:32462/1‑94 (MQ=255) tATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGt < 1:1266462/100‑1 (MQ=255) tATGGACATAAAAAACTAAAGAGCGTTGAAATATAACTATCTCA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATgg < 2:210216/99‑1 (MQ=255) aaagaaaaaaaaCACCAAAGAGAGTTAAGATAAAACAATCT‑AAAAAAAAATATTATTTATAAAAAAAAAAAAAAAAAAAAAATTAGTAGAGGGAATGGt < 2:1403027/91‑1 (MQ=255) gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt < 1:1178075/100‑1 (MQ=255) gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt > 2:713586/1‑100 (MQ=255) gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt > 1:1384231/1‑100 (MQ=255) gACATAAAACACAAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt < 2:606844/100‑1 (MQ=255) aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAttt < 1:20668/100‑1 (MQ=255) cATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtttt < 1:194820/100‑1 (MQ=255) cATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtttt > 1:827712/1‑100 (MQ=255) aTAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTta > 1:449618/1‑100 (MQ=255) aaaaCACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata < 1:571955/100‑1 (MQ=255) aaaCACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGTGGTAATGGTCATTTTATaa > 2:1354480/1‑100 (MQ=255) aaCACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg < 1:1345082/100‑1 (MQ=255) aaCACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg < 1:699974/100‑1 (MQ=255) acacCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATTACCAAAAAATTAGTTGAGGGAATGGTCTTTTTATAAg > 2:1228566/1‑99 (MQ=255) acacCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGt > 1:933306/1‑100 (MQ=255) cacCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta > 1:800046/1‑100 (MQ=255) ccAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtata < 1:1186149/100‑1 (MQ=255) gCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTta > 1:427192/1‑100 (MQ=255) tGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaataa > 2:344450/1‑100 (MQ=255) tGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat < 2:1569995/98‑1 (MQ=255) tGAGATATATCTTTCT‑G‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTaataat < 1:518329/100‑1 (MQ=255) gaTATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa < 2:186189/100‑1 (MQ=255) gaTATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa > 2:1098650/1‑100 (MQ=255) gaTATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa > 1:998740/1‑100 (MQ=255) tatatCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa < 2:1467584/100‑1 (MQ=255) tatatCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa > 2:967604/1‑100 (MQ=255) tCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGAc < 2:449617/100‑1 (MQ=255) tCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGAc < 1:617474/100‑1 (MQ=255) tttCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGc < 2:1384223/100‑1 (MQ=255) ttCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATGAGTAGAGGGAATGGTCATTTTATAAGGATAATGGTTAATAATTACAGGACGc > 2:1260512/1‑99 (MQ=255) | TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGC > NC_000913/1529970‑1530162 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |