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.7 | 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 > 1:702955/1‑143 (MQ=255) tATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACCGGCGTCATTCTAATGACGCTATCCGGAAATGCTACGTTAGCCGTTGATTATCCTGCAACGCTCGCCAAg > 2:453564/1‑143 (MQ=255) cAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGAAtt > 2:128630/1‑143 (MQ=255) gTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATACTGCCAGGCGTGCAAAGAATTTTTTGtga > 1:363462/1‑141 (MQ=255) tGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCGCAAAGAATTTTTTGTAATc < 2:472675/143‑1 (MQ=255) tttCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGAATTTTTTGTAATCcgtc < 1:653491/143‑1 (MQ=255) ttCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCGCACAGAATTTTCTGTAATCcgtcg > 1:646977/1‑143 (MQ=255) aTGAAATTTCCTGATTACAACGGACGTAGCCAACACGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGACTTTTTTGTAATCCGTCGTTGc > 2:447998/1‑143 (MQ=255) ggAAAGTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCa < 1:601463/142‑1 (MQ=255) gagaTTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCaa < 2:97782/140‑1 (MQ=255) gAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCaa > 1:352217/1‑143 (MQ=255) tGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCGCAAAGAATCTATAGTACACCGTCGTTGCAAGCACtcttt > 1:33537/1‑143 (MQ=255) tACAACGGACGTAGCCAACAGGCGAAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAg < 1:413353/143‑1 (MQ=255) cgACATAGCCAACATGCGTATTCTTTAACAGGCGTTATTATATTGACGTGATCAGTAAATGCTACGTTATCCGTTGATTATCCTTCCTCGCTCGCAAAGTATTTTTTGTACTCCGTCGTTTCAAGCACTTTTTCAAGCTGTGc < 1:571610/142‑1 (MQ=255) aGCCAACGAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTTATCATCCTGCTACTCTCGCAAAGAATTTTTTGTAATCCCTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAg < 2:193972/143‑1 (MQ=255) aGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATAACGCTATCAGTAAATGCTACGTTATCCCTTCAGTATCCTGCAACGCTCGCAAAGAATTTCTTGTGACCCGTCGTTGCACGCACGTTTTCCAGCCCTGCcacgcc > 1:106600/1‑140 (MQ=255) aaGCGCAGGCATGAACAGGCGTCATAATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGt > 1:594165/1‑143 (MQ=255) gcAGGCATGACCAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGc < 1:633022/143‑1 (MQ=255) cAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGAATTTTTTGTAATCCGTCGTCGCAAGCATTTTTGCAAGCTGTGCCACAAGCTGACGTAAGcc > 2:204562/1‑143 (MQ=255) gCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTg > 1:440397/1‑143 (MQ=255) gCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTg < 2:200445/143‑1 (MQ=255) aTGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCg < 1:761144/143‑1 (MQ=255) tGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGa > 2:409435/1‑143 (MQ=255) aCGCTATCAGTAAATGCTACGTTATCCGCTTACGAGCCTGCGACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGAcc < 2:514199/143‑1 (MQ=255) cAGTAAATGCCACGTTACCCGTTGATTATCCTGCCACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGAAGACg < 2:732765/143‑1 (MQ=255) attAAATGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTGTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTTCTCGTCTGTAAAGCGACCGAAGACgg < 1:286923/141‑1 (MQ=255) tAACTGCTACGTTAGCCGTTGATTATCCTGCAACGCTCGCAAGGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGAAGACGGTa < 1:278948/143‑1 (MQ=255) tAAATGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGAAGACGGTa < 2:272202/143‑1 (MQ=255) tGCTACGTTATCCGTTGATTATCCTGCAACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGGCCGAAGACGGTACTAt > 2:513013/1‑143 (MQ=255) aTCCGTTGATTATCCTGCGACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGAAGACGGTACTATCGATGTCga > 2:16870/1‑143 (MQ=255) tCCGTTGATTATCCTGCAACGCTCGCAAAGAATTTTTGGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGAAGACGGTACTATCGATGTCgag < 1:323309/143‑1 (MQ=255) gATTATCCTGCGACGCTCGCAAAGAATTTTTTGTAATGCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGAAGACGGTACTATCGGTGGCGAGAACAcc > 2:330827/1‑143 (MQ=255) ttATCCTGCGACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGAAGACGGTACTATCGATGTCGAGAACACCaa < 2:353236/143‑1 (MQ=255) ttATCCTGCAACGCTCGCAAAGACTTTTTTATAATCCGGCGTTGCAAGCACTTTTTCAAGCTGTGCAACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGAATACGGTACTATCGATGTCGAGAACACCaa < 1:50881/143‑1 (MQ=255) tATCCTGCAACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGAAGACGGTACTATCGATGTCGAGAACACCaat < 1:292782/143‑1 (MQ=255) aTCCTGCAACGCTCGCAAAGAATTTTTTGTAATCCGTCGTTGCAAGCACTTTTTCAAGCTGTGCCACAAGCTGACGTAAGCCTTGCTCGTCCTCGTCTGTAAAGCGACCGCAGACGGCACTATCGATGTCGCGAACACCacta > 1: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: ‑ |