Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,227,996 | 0 | C | G | 53.8% | ‑6.7 / 11.6 | 26 | R433P (CGG→CCG) | PP_4603 | transposase |
Reads supporting (aligned to +/- strand): ref base C (0/12); new base G (0/14); total (0/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GACAGCGCAGCCCGGAAGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGG > NC_002947/5227896‑5228137 | caccgcgcTGCCCGGACGGGCGGGGTGCTCTCACACCAAAAAAATGGCAACCCATTGTTTGATAAGCACTCCAGCGGCCACCAGTTGATACCCCTGTTTCGGGGGTTTGCTAAGCACAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTc < 2:2827834/146‑1 (MQ=1) caccgcgcAGCCCGTAAGGGCTGGGTGATCTCTAATCAAAAAAAATGCACCCCATGTTTTTCTCCGCACACACGGGGCCCCCAGTGGATACCCCCGTTGCGGGGGTTTGCCACCAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTc < 2:6015366/146‑1 (MQ=0) caccgcgcAGCCCGGATGGGGTGGTTTATCTCACAAAAAAAAAATTGCAACCCTTTTTTTTCTAACCACTCACGCGGCAACCATTTGATACCCCTTTTCCGGGGATTTGCTAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTc < 2:195965/146‑1 (MQ=0) gccgCCGGGAGGGGCGGGGTTATCTCTCATCACAAAAATTGCACCCCTTTGTTTTCTACGCCCCCACGGGCCCACCCTTTTATACCCCCTTTGCGGGGGTTTGCGAAGAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAActc < 2:174470/147‑1 (MQ=0) gccgCCCGGATGGGGTGGTTGATATCTCAACAAACAAAATGCAACCCATTGTTTTCTACCCCCTCACGCGGCCACCAGTTGATACCCCTTTTGCGGGGCTTTGCGAAGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAActc < 2:174378/147‑1 (MQ=1) ccggcgggggTGTGTTATCTCCCATAAAAAAAATTGAAACCATTTTTTTTCTAAGCACTCCCGGGGACACCTTTTTATACCCCTTTTGCGGGGTTTTGCCACCAAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCaa < 2:1647701/144‑1 (MQ=0) gatgggCTGGTTTATCTCTAATAAAAAAAAATGAAACCCTTTTTTTTCTAACCACTCACGCGGCCACCTGTTGATACCCCTGTTTCGGGGCTTTGCTAAGAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGAc < 2:2714995/147‑1 (MQ=1) actcaaaaaacaaaATGCAACCCTTTTTTTTCTAAGCACTCACGCGGCCACCTGTTGATACCCCCGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGc < 2:6150908/143‑1 (MQ=1) cacaaTGCACCCCTTTGTTTTCTCCGCCCCCCCGCGGCCACCCGTTGCCCCCCCCGTTGCGGGGCTTTGCGAAGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTg < 2:2667775/147‑1 (MQ=1) caccccTTTGTTTCCAAACCACCCACGCGGCCCCCCGTTGATCCCCCCGTTCCGGGGCTTTGCGACGAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa < 2:1086053/147‑1 (MQ=1) cacccaaTTTTTTTCTAACCACCCCCGCGCCCCCAATTCTATAACCCGTTTGCCGGGATTTCCTAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa < 1:2191346/144‑1 (MQ=0) ccttttttttCTAAGCACTCACGCGGCCACCTGTTGATACCCCAGTTGCCGGGCTTTGCGAAGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt < 1:7937406/144‑1 (MQ=1) ccttttttttCTAAGCACTCACGCGGCCACCATTTGATACCCCTGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt < 2:4487889/144‑1 (MQ=1) cctttGTTTTCTACGCACTCACGCGGCCACCTTTTTATACCCCTTTTTCCGGGCTTTGCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt < 1:4834784/147‑1 (MQ=1) cctttGTTTTCTAAGCACTCAAGCGGCCCCCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt < 2:8345389/147‑1 (MQ=1) tttttttCCAAGCACTAAAGCGGCCACCAGTTCAAACCCCAGTTTCGGGGCTTTGCTACGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTa < 1:5514451/147‑1 (MQ=1) atcaatCAGGCGGCCACCTTTTTATACCCCTGTTTCGGGGTTTTGCTAAGAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGt < 1:6090861/145‑1 (MQ=1) acccacgCGGCCACCCGTTGATACCCCTCTTTCCGGGATTTGCTACGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGcc < 1:503951/144‑1 (MQ=1) cacgcgcCCACCTGTTGAGCCCCCTGTGTCGGGGCTTTGCCACCAAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAAc < 1:7672437/147‑1 (MQ=1) cgcccaccaGTTCATACCCCAGTTTCCGGGATTTCCTAACCAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGAt < 2:2526519/147‑1 (MQ=1) ccctgTTGACCCCCCTGTTCCGGGGTTTTGCTAACAAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCa < 1:5227752/146‑1 (MQ=1) ccccgTTTATCCCCCTTTTTCGGGGCTTTGCGACGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCa < 1:1456407/146‑1 (MQ=1) tttatcccccAGTTTCCGGGCTTTGCTAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGc < 1:1524664/144‑1 (MQ=0) tttaTACCCCTTTTTCCGGGCTTTGCTAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCTGTGGc < 2:3562734/147‑1 (MQ=1) cctgTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATgg < 1:5553673/147‑1 (MQ=1) cccgTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATgg < 2:9069027/147‑1 (MQ=1) | GACAGCGCAGCCCGGAAGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGG > NC_002947/5227896‑5228137 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |