Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,859,104 | 0 | T | G | 46.5% | 40.4 / 17.7 | 43 | M671R (ATG→AGG) | mutS | methyl‑directed mismatch repair protein |
Reads supporting (aligned to +/- strand): ref base T (16/7); new base G (0/20); total (16/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.43e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGC > NC_000913/2858954‑2859254 | gTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTAt < 1:425441/151‑1 (MQ=255) ccTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTAtgg > 2:304082/1‑150 (MQ=255) cTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTAtggtg > 1:302095/1‑151 (MQ=255) tataTGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTAtggtgg > 1:305874/1‑151 (MQ=255) atGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTAtggtgg > 2:132901/1‑148 (MQ=255) aCCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGaa > 2:858774/1‑150 (MQ=255) gATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAat > 2:626754/1‑151 (MQ=255) cTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTAc > 1:903734/1‑151 (MQ=255) gATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACAt > 1:235692/1‑151 (MQ=255) ggCCTACATCGGCAGTTATGTACCGGCAAAAAAAGTCGAGATTGGACCTATCGATCGCATTTTTCCCCGCGTGGGCGCGGCAGATGCCCTGGCGCCCGGGCGCTCACCTTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACAt < 2:7953/148‑1 (MQ=255) ccTACATCGGCAGTTATTTCCCGGCAAAAAAGGTCGAGATTGGCCCTATCGACCGCTTCTTTCCCCGCGGGGGCGCGGCAGATGACCTGGCGCCGGGGCCCTCAACTTTTAGGGTGAAGATGACTAAAACCGCCAATATTTTACATAACGc < 2:945741/151‑1 (MQ=255) acacCGGCATTTATGTCCCGGCAAAAAAAGTGGAGATTGCACCTACCGACCGCATTTTTCCCCGCGGGGGCGCGGAAGATGCCCGGGCTCCCGGGCGCCCACCCTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCa < 1:395562/150‑1 (MQ=255) tATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGt > 1:461256/1‑150 (MQ=255) tATGCACGGGCAAAAAAAGTCGAGTTTGGCCCTATCGACCGCTTTTTTCCCCGGGTAGGCGCGGCAAATGCCCTGGCGCCCGGGCGCTCACCTTTTATGGTGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGt < 2:663441/150‑1 (MQ=255) tGTCCCGGCAAAAAAAGTCGAGTTTGGCCCTATCGTCCGCTTCTTCCCCCGCGGGGGCGGGGCAGATGCCCGGGCGTCGGGGCCCCAACCCTTTAGGGGGGAGATGAATAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCt < 1:958488/150‑1 (MQ=255) gTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTgg > 1:110890/1‑151 (MQ=255) gTACCGGCAAAAAAGGTGGAGATTGGCCCTATCGATCCCTTCTTTCCCCGGGTGGGCGCGGCAGATGCCCGGGCGTCCGGGCGCCCAACCTTTATGGTGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACa < 1:367074/145‑1 (MQ=255) taccgcaaaaaaaaGTCGAGATTGGCCCTATCGCCCGCTTTTTTCCCCGCGTAGGCGCGCCAGAGGACCGGGCGTCCGGGCGCCCACCCTTTAGGGTGGAAATGATTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCt < 2:860816/139‑1 (MQ=255) tACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTgg > 2:4333/1‑150 (MQ=255) cGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGtt > 1:889454/1‑151 (MQ=255) ggaacaAAAAGTGGAGATTGCCCCTACCGATCGCTTTTTTCCCCGGGGGGGCGGGGAAGATGCCCGGGCGCCCGGGCCCTAACCCTTTAGGGTGGAGATGACTAAAGCCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGtt < 1:699038/147‑1 (MQ=255) aaaaaaaGTCGAGATTGACCCTAGCGATCGCATCTTTCCCCGCGTAGGCGCGGCAGATGCCCTGGCGCCCGGGCGCTCACCTTTTAGGGGGGGGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGgt < 2:899636/142‑1 (MQ=255) aaaaGTCGAAGTTGGCCCTCCCGACCCCATTTTTCCCCGGGTGGGCGCGGCAGAGGCCCGGGCGCCGGGGCCCTAACCCTTTATGGTGGAGATGCTGGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGAt < 1:130784/151‑1 (MQ=255) cGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATg > 1:480985/1‑146 (MQ=255) ggACCTATCGACCGCTTCTTCCCCCGCGTGGGCGCGGCAGATGCCCTGGCGCCGGGGCGCTCAACCTTTAGGGTGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGt < 2:627913/150‑1 (MQ=255) cttcccttgaattCTTTCCCCGCGAGGGGGCGGAGGTGGACCGGGCGGCCGGCCCCACAACTTTTTTGGGGGAGAGGAAGAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGt < 1:911195/134‑1 (MQ=255) cgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCAcc > 2:100471/1‑150 (MQ=255) cATCTTTCCCCGCGGCGGCGCGGCAGATGCCCTGGCGCCCGGGCGCTCACCCTTTGTGGGGGAGATTACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt < 2:562112/149‑1 (MQ=255) cATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTAc > 2:425766/1‑151 (MQ=255) aTCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACg > 1:179654/1‑151 (MQ=255) tttACCCGCGTGGGCGCGGCAGAGGACCGGGCGTCGGGGCGCTCAACTTTTAGGGGGGAGATGACGGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTAc < 2:57313/147‑1 (MQ=255) cccGCGTGGGCGCGGCAGATGCCGGGGCGCCCGGCCGCCCACCTTTTATGGTGGAAATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCt < 2:305882/151‑1 (MQ=255) ccGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTg > 2:838516/1‑151 (MQ=255) ggCGCGGCAGATGCCCTGGCGTCGGGGCCCTCAACTTTTAGGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCg < 2:354570/147‑1 (MQ=255) ggCGCGGCAGATGCCCGGGCGCCGGGGCCCTCACCTTTTAGGGGGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTgg < 2:947231/151‑1 (MQ=255) gcgcGGCGCATGCCCGGGGGTCCGGTCGCTCCACTTTTATGGGGGAGAGGATTGAATCCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTg < 2:903749/149‑1 (MQ=255) gCAGGTGACCGGGCGGCCGGACGCCCAACTTTTTGGGTGGAAGTCACAAAAGCCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtgg < 1:938827/150‑1 (MQ=255) aTGCCCGGGCGTCGGGGCCCCAAACTTTTAGGGTGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtgg < 1:645950/146‑1 (MQ=255) tGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtgggcgtg > 1:618209/1‑150 (MQ=255) tGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTgc > 2:938464/1‑151 (MQ=255) ctggcctcccGGGCGCTCAACTTTCAGGGTGGGGATACTTTAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTgc < 2:892947/140‑1 (MQ=255) ccTGGCGTCCGGGCGCTCAACCTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTgc < 2:350030/148‑1 (MQ=255) ggCGCCCGGGCCCTCAACCTTTGGGGTGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGaa < 2:302103/151‑1 (MQ=255) gcgggccccTCACCTTTTAGGGGGGAGAGTAATGAATCCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAAt < 2:273549/142‑1 (MQ=255) ggCGCTCAACCTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGc < 2:914856/151‑1 (MQ=255) cgcTCAACTTTTAGGGTGGAGATGACAAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGc < 1:603672/149‑1 (MQ=255) tttAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGAtt < 2:894415/150‑1 (MQ=255) ttgttggtggAAATGTCTAAATCCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGa < 1:815833/138‑1 (MQ=255) ttATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTa > 2:316272/1‑150 (MQ=255) tggtggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGc > 1:774832/1‑151 (MQ=255) | GTACCTATATGCGCCAGACCGCACTGATTGCGCTGATGGCCTACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGC > NC_000913/2858954‑2859254 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |