| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | REL606 | 3,145,726 | G→A | intergenic (+57/‑138) | dnaG → / → rpoD | DNA primase/RNA polymerase sigma factor |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 3,145,726 | 0 | G | A | 98.6% | 272.3 / ‑3.7 | 70 | intergenic (+57/‑138) | dnaG/rpoD | DNA primase/RNA polymerase sigma factor |
| Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (35/34); total (35/35) | |||||||||||
| 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.99e-01 | |||||||||||
AGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTA > REL606/3145653‑3145797 | aGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGc < 1:1391721/76‑1 (MQ=255) gCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCa < 1:1157538/76‑1 (MQ=255) tGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAagcagc < 2:1923226/76‑1 (MQ=255) tGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAagcagc > 1:1491501/1‑76 (MQ=255) tGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAagcagc > 2:507840/1‑76 (MQ=255) tGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAagcagc < 2:1497246/76‑1 (MQ=255) ggCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCg < 1:768230/76‑1 (MQ=255) gAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCa > 2:1819643/1‑76 (MQ=255) gAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCa > 1:1431403/1‑76 (MQ=255) aaaaaGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCaa > 1:404367/1‑76 (MQ=255) aaaGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAat < 2:1822995/76‑1 (MQ=255) aaaGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAat < 2:1670279/76‑1 (MQ=255) aatGATTTACCGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAatat < 2:1312351/74‑1 (MQ=255) aGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAatat < 2:408065/76‑1 (MQ=255) aGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAatat < 2:291321/76‑1 (MQ=255) aGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAatat > 1:1399278/1‑76 (MQ=255) tGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAatatat > 2:1438819/1‑76 (MQ=255) gATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAAtatata < 1:1946400/76‑1 (MQ=255) gATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAAtatata < 1:1403852/76‑1 (MQ=255) aTTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATaa > 1:762769/1‑76 (MQ=255) aTTTAACGGCGTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATaa < 2:1284531/76‑1 (MQ=255) tAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTa < 2:859781/76‑1 (MQ=255) cGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGc > 2:662278/1‑76 (MQ=255) cTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCt > 2:253649/1‑76 (MQ=255) tAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCg > 2:1529073/1‑76 (MQ=255) tAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCg > 1:58758/1‑76 (MQ=255) tAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCg < 1:307081/76‑1 (MQ=255) tAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCg < 1:1926282/76‑1 (MQ=255) aaGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGt > 1:231872/1‑76 (MQ=255) aaGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGt < 1:1078657/76‑1 (MQ=255) aGTGCGGAAGCGCGACCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTa < 2:1625436/76‑1 (MQ=255) gCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAAtt < 2:489881/76‑1 (MQ=255) gCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAAtt > 2:1421616/1‑76 (MQ=255) gCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAAtt > 2:1615975/1‑76 (MQ=255) ccGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTa > 2:1795789/1‑76 (MQ=255) cGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTAt > 2:1545208/1‑76 (MQ=255) cGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTAt > 2:613184/1‑76 (MQ=255) cGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTAt < 2:863915/76‑1 (MQ=255) gAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATc < 1:1958362/76‑1 (MQ=255) agagCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGt > 2:1875650/1‑76 (MQ=255) gagCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGtt < 2:1130525/76‑1 (MQ=255) gagCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGtt > 2:1941060/1‑76 (MQ=255) gagCGATCGGGAAG‑CCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTg > 2:926621/1‑76 (MQ=255) agcgagcgGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTg < 2:1333345/76‑1 (MQ=255) agCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTg < 1:507936/76‑1 (MQ=255) agCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATAATCGTTg > 2:1662276/1‑76 (MQ=255) gCGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTgg > 2:703616/1‑76 (MQ=255) cGATCGGGAAGCCCCCGACAGCCGCTCTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGc > 2:923001/1‑76 (MQ=255) cGATCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGc > 2:163998/1‑76 (MQ=255) aTCGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCgg < 1:682777/76‑1 (MQ=255) cGGGAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTa < 2:1059935/76‑1 (MQ=255) gggAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTaa < 2:759865/76‑1 (MQ=255) ggaaaGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTaa < 1:38171/73‑1 (MQ=255) gAAGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAaac < 1:1818826/76‑1 (MQ=255) aaGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAaaca < 2:312703/76‑1 (MQ=255) aaGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAaaca > 2:603446/1‑76 (MQ=255) aGCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAaacaa < 1:438685/76‑1 (MQ=255) gCCCCCGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAaacaac > 2:1017928/1‑76 (MQ=255) ccccGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCg < 2:1631286/76‑1 (MQ=255) cccGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGt > 2:81542/1‑76 (MQ=255) ccGACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGtt > 1:802087/1‑76 (MQ=255) gACAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTgg > 2:1985956/1‑76 (MQ=255) aCAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGa < 2:582544/76‑1 (MQ=255) cAGCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGAt > 2:1589110/1‑76 (MQ=255) gCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGAttt > 1:1188410/1‑76 (MQ=255) gCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGAttt < 2:372991/76‑1 (MQ=255) gCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGAttt < 2:628199/76‑1 (MQ=255) gCCGCACTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGAttt < 2:14929/76‑1 (MQ=255) aCTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCg > 2:556896/1‑76 (MQ=255) aCTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCg < 2:435824/76‑1 (MQ=255) aCTGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCg > 1:1739609/1‑76 (MQ=255) tGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGtt > 2:366122/1‑76 (MQ=255) tGAGAAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGtt > 2:1308255/1‑76 (MQ=255) gagaAGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTa > 1:737483/1‑76 (MQ=255) | AGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTA > REL606/3145653‑3145797 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |