Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 1,708,758 | G→C | 54.7% | E614Q (GAA→CAA) | rsxC → | SoxR [2Fe‑2S] reducing system protein RsxC |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,708,758 | 0 | G | C | 54.7% | ‑6.6 / 9.5 | 11 | E614Q (GAA→CAA) | rsxC | SoxR [2Fe‑2S] reducing system protein RsxC |
Reads supporting (aligned to +/- strand): ref base G (3/2); new base C (2/4); total (5/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.67e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.05e-01 |
CCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCC > NC_000913/1708609‑1708904 | cccGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAGCAGAACAACAGGGCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCagaag < 1:1294160/150‑1 (MQ=255) ataacgAAGCAGAACCACAGGATGATCCTCGCAAACCCGCCGTCGAAGCAGCTATTGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTa < 1:1561769/145‑1 (MQ=14) cAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGAGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTCGATCCGAGCAACGACGACGTCGAAGCCGCTATTGCCAGTGCAAAAGCACGCAAGCGGGAACa > 2:1494670/1‑150 (MQ=14) cccGTGAAAAATCTCGCAAGCAGGAACAGCAACAGGATAATCCGGAAGCAGAACAACAGGTCGATCCCCGCAAACCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaac < 1:1469556/150‑1 (MQ=11) ccAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTg < 1:206030/96‑1 (MQ=33) ccAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTg > 2:206030/1‑96 (MQ=33) aGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCa > 1:950521/1‑150 (MQ=255) aaCAGGCTAATGCGCAACCAGAACACCCGGGCGCTCCGCGCAACGCCGCCGTCGACGCCGCTATAGCCCGTGCCAACGCACGCAAGCGGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCg < 1:1231344/150‑1 (MQ=11) cGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGa < 1:1548593/100‑1 (MQ=18) cGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGa > 2:1548593/1‑100 (MQ=18) cAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa > 1:1548910/1‑100 (MQ=37) cAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa < 2:1548910/100‑1 (MQ=37) gaagaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAAGCAAATGCCTAGCCAGAAAAACAGGTCGATCCGCGCAAAGCCGCAGTCGAAGCAGCTATTGACCGAGca > 2:592628/1‑149 (MQ=14) gaagaaCAGGTCGATCCGCGCAAAGCCGCAGTCGAAGCCGATATTGCCCGAGACAAAGCACGCAAGATGGAAAAGCAAAAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCCGCAGTCGAAGCCGGTATTGACGGAGcc > 2:1269155/1‑150 (MQ=14) | CCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCC > NC_000913/1708609‑1708904 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |