Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,328,951 | A→G | S10P (TCC→CCC) | mreB ← | cell wall structural complex MreBCD, actin‑like component MreB |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,328,951 | 0 | A | G | 95.5% | 260.2 / 1.5 | 67 | S10P (TCC→CCC) | mreB | cell wall structural complex MreBCD, actin‑like component MreB |
Reads supporting (aligned to +/- strand): ref base A (1/2); new base G (40/24); total (41/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.55e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01 |
TTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCG > REL606/3328878‑3329021 | ttCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGgaa < 2:866807/76‑1 (MQ=255) ttCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGgaa > 1:1156165/1‑76 (MQ=255) ttCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGgaa > 2:320427/1‑76 (MQ=255) aCGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGc > 2:820798/1‑76 (MQ=255) aCGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGc > 2:1612931/1‑76 (MQ=255) aCGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGc > 2:1168427/1‑76 (MQ=255) cGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGcc < 1:489397/76‑1 (MQ=255) gATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCa > 2:373138/1‑76 (MQ=255) aTGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCAc > 1:198797/1‑76 (MQ=255) tgtcTTGTTCTTTTAGTTAAATGAGGGTGTTCGCAGTTCCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACg < 2:1387399/73‑1 (MQ=255) ccTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGaa > 1:564137/1‑76 (MQ=255) ccTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGaa < 1:411470/76‑1 (MQ=255) ccTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGaa > 1:969454/1‑76 (MQ=255) ccTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGaa < 2:1203979/76‑1 (MQ=255) cTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGaaa > 1:1331221/1‑76 (MQ=255) gTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAAttt > 1:858630/1‑76 (MQ=255) tCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAAtttt > 1:1337760/1‑76 (MQ=255) ccTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAACttttt > 1:293831/1‑76 (MQ=255) ccTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAAttttt > 1:293819/1‑76 (MQ=255) cTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAAtttttt > 1:902170/1‑76 (MQ=255) tttACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCa > 2:149792/1‑76 (MQ=255) tACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAAc > 1:1609562/1‑76 (MQ=255) tACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAAc > 1:1612778/1‑76 (MQ=255) tACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAAc > 1:933182/1‑76 (MQ=255) tACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAAc > 1:752102/1‑76 (MQ=255) tACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAAc > 1:665169/1‑76 (MQ=255) tACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAAc > 1:1041827/1‑76 (MQ=255) cATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACac > 2:433383/1‑75 (MQ=255) tAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATAc < 1:1064735/76‑1 (MQ=255) aaTGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTa < 1:684000/76‑1 (MQ=255) aTGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTaa < 1:178053/76‑1 (MQ=255) gAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAgg > 2:505743/1‑76 (MQ=255) gAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAgg < 1:1591985/76‑1 (MQ=255) aGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAggg > 2:1101476/1‑76 (MQ=255) gggTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGa < 2:1421194/76‑1 (MQ=255) gTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATa < 1:550462/76‑1 (MQ=255) tATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATaa > 2:1524256/1‑76 (MQ=255) aTTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAAt > 2:170315/1‑76 (MQ=255) aTTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAAt > 2:163756/1‑76 (MQ=255) cGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCt > 2:224528/1‑76 (MQ=255) cGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCt > 1:571906/1‑76 (MQ=255) gCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTg < 1:1130432/76‑1 (MQ=255) gCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTg > 1:1374797/1‑76 (MQ=255) taGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGa > 2:1025146/2‑76 (MQ=255) aGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGaa < 1:135914/76‑1 (MQ=255) aGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGaa < 2:1168382/76‑1 (MQ=255) gtgcccAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGaaa < 2:1439796/73‑1 (MQ=255) gTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGaaa > 1:1374087/1‑76 (MQ=255) gTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGaaa < 2:1067699/76‑1 (MQ=255) tACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAg > 2:894834/1‑76 (MQ=255) aCCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGc > 1:1474558/1‑76 (MQ=255) aCCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGc > 1:699312/1‑76 (MQ=255) aCCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGc > 2:177639/1‑76 (MQ=255) cccAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCt < 1:44379/76‑1 (MQ=255) aGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTggg < 1:1287688/76‑1 (MQ=255) aGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTggg < 1:1328818/76‑1 (MQ=255) aGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTggg < 1:204216/76‑1 (MQ=255) tCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCg < 2:584661/76‑1 (MQ=255) tCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCg > 1:1558427/1‑76 (MQ=255) cAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCgg < 2:1180116/76‑1 (MQ=255) aaTGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGga < 2:19207/76‑1 (MQ=255) tGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaa < 1:189576/76‑1 (MQ=255) tgACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaa < 2:46695/75‑1 (MQ=255) ggACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaa > 1:497844/1‑76 (MQ=255) aCAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaaga < 1:1430807/76‑1 (MQ=255) cAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaagaa > 2:1626276/1‑76 (MQ=255) cAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaagaa > 2:1614445/1‑76 (MQ=255) cAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaagaa > 1:148916/1‑76 (MQ=255) cAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaagaa > 2:1260576/1‑76 (MQ=255) aaGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaagaaa > 2:158316/1‑76 (MQ=255) gTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAAt < 1:1622236/76‑1 (MQ=255) cATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATcc > 2:1360436/1‑76 (MQ=255) aTTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCg < 2:1174774/76‑1 (MQ=255) | TTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCG > REL606/3328878‑3329021 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |