Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,327,541 | G→A | S110F (TCC→TTC) | mreC ← | cell wall structural complex MreBCD transmembrane component MreC |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,327,541 | 0 | G | A | 97.8% | 172.2 / NA | 46 | S110F (TCC→TTC) | mreC | cell wall structural complex MreBCD transmembrane component MreC |
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (22/23); total (22/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
ATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCA > REL606/3327471‑3327613 | aTCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGAAAccc < 1:985966/76‑1 (MQ=255) cGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGAAACCcag > 1:983074/1‑76 (MQ=255) gCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGAAACCcagc < 1:1049643/76‑1 (MQ=255) cTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGAAACCcagca > 2:289828/1‑76 (MQ=255) cTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGAAACCcagca < 1:289828/76‑1 (MQ=255) tataAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGAAACCcagcag < 2:436455/76‑1 (MQ=255) tataAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGAAACCcagcag < 1:932709/76‑1 (MQ=255) aaGGATCGTTAACCGTGGAGATAACTTGAGTTACCATTTANNGCTCATCCTGNNGTAGCGGGAAACCCAGCAGCTc > 1:1264782/1‑76 (MQ=255) aGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTcg > 2:414898/1‑76 (MQ=255) aTCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCgcgc > 2:567899/1‑76 (MQ=255) tCGTTAACNGTGGAGATAANCTGAGTCACCATTTTCTGCTCANCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCa < 2:1173447/76‑1 (MQ=255) tAACCGTGGTGATAACCTGAGTCNCCATTNTCTGCGCATCCTGANGCAGCGAGAAACCCAGCAGCTCGCGCAGAcg > 2:178093/1‑76 (MQ=255) tAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCCGAcc > 1:864564/1‑75 (MQ=255) tAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGAcg > 1:102476/1‑76 (MQ=255) tAACCGTGGAGATAACCTGAGTCACCATATTCTGATCTTCCTGACGCNGCGGGAAACCCAGCAGCTCGCGCAGAcg > 1:119390/1‑76 (MQ=255) aCCGTGGAGATAACCTGAGTCACCATTTTTTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGAcgcg < 2:839951/76‑1 (MQ=255) gTGGAGATANCCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGt > 1:36786/1‑76 (MQ=255) tgagaTAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTtc > 2:925991/2‑76 (MQ=255) gagaTAACCTGAGTCGCCTTTTGCTTCTCATCCTGACGCGGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTtct < 2:1225714/76‑1 (MQ=255) gagaTAACCTGAGTCGCCATTTTCTGCTCGTCCCGGCGCAGCGGGAATCCCAGGAGCTCGCGCAGACGCGCGTtct < 2:1068982/76‑1 (MQ=255) agaTAACTTGAGTCACCATTTTNNGCTCATCCTGNNGCAGCGGGGAACCCGGCAGCTCGCGCAGACGCGCGTtctc < 1:1257529/76‑1 (MQ=255) aTAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCt < 2:465345/76‑1 (MQ=255) aTAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCt > 2:151092/1‑76 (MQ=255) tAACCTGACTCTCCATTTTCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTg < 2:540929/76‑1 (MQ=255) tGATTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGt < 2:146113/76‑1 (MQ=255) tGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCCGACGCGCGGTCTCCTGTTTGt > 2:175733/1‑76 (MQ=255) tGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGt > 1:1129354/1‑76 (MQ=255) gTCACCATTTTCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTAtt < 2:980311/76‑1 (MQ=255) aCCATTTTCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGGTCTCCTGtttgttttgtt > 1:914104/1‑75 (MQ=255) cgttttCGGATAATCCGGAAGCAGCGGGAAACCAATGATCTCGCGCAGACGAGCGTGCTCATGTTTGTATTGTACa < 2:727199/74‑1 (MQ=255) cATTTTCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCa < 2:727225/76‑1 (MQ=255) cATTTTCTGCTCATCCTGACGCAGCAGGAAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTATCCa < 2:1374460/76‑1 (MQ=255) aTTTTCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCaa < 1:557351/76‑1 (MQ=255) aTTTTCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCaa > 2:576198/1‑76 (MQ=255) ttttCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAg < 2:504701/76‑1 (MQ=255) tttCTGCTCATCCTGACGCAGCGGGAAACCCAGNAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGc > 2:1182487/1‑76 (MQ=255) tttCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTTCTTCTGTTTGTATTGTCCAAga > 2:990498/1‑75 (MQ=255) tttCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGc > 2:307526/1‑76 (MQ=255) tttCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGc < 1:907120/76‑1 (MQ=255) tttCTGCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGGTCTCCTGTTTGTATTGTCCCAGc > 1:311829/1‑76 (MQ=255) gCTCATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATcag > 2:855425/1‑76 (MQ=255) cATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATcagcag < 2:141950/76‑1 (MQ=255) cATCCTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATcagcag < 1:296232/76‑1 (MQ=255) aTCCTGACGCAGCGGGAAACCNAGCAGNTCGTGCAGACGCGAGTTCTCNTGCGTGTATTGCCCGACCATCTGCAGt > 1:1357094/1‑76 (MQ=255) ccTGACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTc < 1:1197687/76‑1 (MQ=255) gACGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACt < 1:601660/76‑1 (MQ=255) cGCAGCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGt < 2:685587/76‑1 (MQ=255) gccgCGGGAAACCCGGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGtt < 2:46541/73‑1 (MQ=255) gCGGGAAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGttttt > 1:482737/1‑76 (MQ=255) cGGGAAACCCAGCANCTCGCGCAGACGCACGTTCNCCTGTCTGNATTGTCCAAGCATCAGCAGCTCACTGTTTTTc < 2:1158209/76‑1 (MQ=255) gggaAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCa < 1:552239/76‑1 (MQ=255) | ATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACTGTTTTTCA > REL606/3327471‑3327613 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |