Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,915,640 | 0 | G | A | 56.2% | 2.2 / 32.8 | 36 | V163V (GTC→GTT) | msrC | free methionine‑(R)‑sulfoxide reductase |
Reads supporting (aligned to +/- strand): ref base G (8/6); new base A (8/10); total (17/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.22e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.04e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGAAGACGGTACTATCGATGTCGAGAACACCAATA > NC_000913/1915507‑1915775 | tGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGTCCTGAACAGGCGTCATTATAATGACGCCATCAGTAAATGCTACGCTATCCGTTGATTATCCCGCGACGCTCGCa > 2:702955/1‑143 (MQ=255) tATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACCGGCGTCATTCTAATGACGCTATCCGGAAATGCTACGTTAGCCGTTGATTATCCTGCAACGCTCGCCAAg > 1:453564/1‑143 (MQ=255) cAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGAAtt > 1:128630/1‑143 (MQ=255) gTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATACTGCCAGGCGTGCAAAGAATTTTTTGtga > 2:363462/1‑141 (MQ=255) tGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCGCAAAGAATTTTTTGTAATc < 1:472675/143‑1 (MQ=255) tttCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGAATTTTTTGTAATCcgtc < 2:653491/143‑1 (MQ=255) ttCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCGCACAGAATTTTCTGTAATCcgtcg > 2:646977/1‑143 (MQ=255) aTGAAATTTCCTGATTACAACGGACGTAGCCAACACGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGACTTTTTTGTAATCCGTCGTTGc > 1:447998/1‑143 (MQ=255) ggAAAGTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCa < 2:601463/142‑1 (MQ=255) gagaTTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCaa < 1:97782/140‑1 (MQ=255) gAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCaa > 2:352217/1‑143 (MQ=255) tGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCGCAAAGAATCTATAGTACACCGTCGTTGCAAGCACtcttt > 2:33537/1‑143 (MQ=255) tACAACGGACGTAGCCAACAGGCGAAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAg < 2:413353/143‑1 (MQ=255) cgACATAGCCAACATGCGTATTCTTTAACAGGCGTTATTATATTGACGTGATCAGTAAATGCTACGTTATCCGTTGATTATCCTTCCTCGCTCGCAAAGTATTTTTTGTACTCCGTCGTTTCAAGCACTTTTTCAAGCTGTGc < 2:571610/142‑1 (MQ=255) aGCCAACGAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTTATCATCCTGCTACTCTCGCAAAGAATTTTTTGTAATCCCTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAg < 1:193972/143‑1 (MQ=255) aGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATAACGCTATCAGTAAATGCTACGTTATCCCTTCAGTATCCTGCAACGCTCGCAAAGAATTTCTTGTGACCCGTCGTTGCACGCACGTTTTCCAGCCCTGCcacgcc > 2:106600/1‑140 (MQ=255) aaGCGCAGGCATGAACAGGCGTCATAATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGt > 2:594165/1‑143 (MQ=255) gcAGGCATGACCAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGc < 2:633022/143‑1 (MQ=255) cAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGAATTTTTTGTAATCCGTCGTCGCAAGCATTTTTGCAAGCTGTGCCACAAGCTGACGTAAGcc > 1:204562/1‑143 (MQ=255) gCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTg > 2:440397/1‑143 (MQ=255) gCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTg < 1:200445/143‑1 (MQ=255) aTGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCg < 2:761144/143‑1 (MQ=255) tGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGa > 1:409435/1‑143 (MQ=255) aCGCTATCAGTAAATGCTACGTTATCCGCTTACGAGCCTGCGACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGAcc < 1:514199/143‑1 (MQ=255) cAGTAAATGCCACGTTACCCGTTGATTATCCTGCCACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGAAGACg < 1:732765/143‑1 (MQ=255) attAAATGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTGTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTTCTCGTCTGTAAAGCGACCGAAGACgg < 2:286923/141‑1 (MQ=255) tAACTGCTACGTTAGCCGTTGATTATCCTGCAACGCTCGCAAGGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGAAGACGGTa < 2:278948/143‑1 (MQ=255) tAAATGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGAAGACGGTa < 1:272202/143‑1 (MQ=255) tGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGGCCGAAGACGGTACTAt > 1:513013/1‑143 (MQ=255) aTCCGTTGATTATCCTGCGACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGAAGACGGTACTATCGATGTCga > 1:16870/1‑143 (MQ=255) tCCGTTGATTATCCTGCAACGCTCGCAAAGAATTTTTGGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGAAGACGGTACTATCGATGTCgag < 2:323309/143‑1 (MQ=255) gATTATCCTGCGACGCTCGCAAAGAATTTTTTGTAATGCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGAAGACGGTACTATCGGTGGCGAGAACAcc > 1:330827/1‑143 (MQ=255) ttATCCTGCGACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGAAGACGGTACTATCGATGTCGAGAACACCaa < 1:353236/143‑1 (MQ=255) ttATCCTGCAACGCTCGCAAAGACTTTTTTATAATCCGGCGTTGCAAGCACTTTTTCAAGCTGTGCAACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGAATACGGTACTATCGATGTCGAGAACACCaa < 2:50881/143‑1 (MQ=255) tATCCTGCAACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGAAGACGGTACTATCGATGTCGAGAACACCaat < 2:292782/143‑1 (MQ=255) aTCCTGCAACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGCAGACGGCACTATCGATGTCGCGAACACCacta > 2:353759/1‑143 (MQ=255) | TGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGAAGACGGTACTATCGATGTCGAGAACACCAATA > NC_000913/1915507‑1915775 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |