Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,223,047 | 0 | C | G | 55.2% | ‑5.5 / 18.3 | 29 | R433P (CGG→CCG) | PP_4603 | transposase |
Reads supporting (aligned to +/- strand): ref base C (0/13); new base G (0/16); total (0/29) | |||||||||||
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.97e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCA‑TTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGGCGC > NC_002947/5222945‑5223191 | gctactgcgcAGCCCGGATGGGCTGGGTGATCTCACACCAAACAAAATGAAACCCT‑TTGTTTTCTAAGCACTCAAGCGGCCACCATTTGATACCCCTGTTTCGGGGATTTGCTAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGtt < 2:3695839/144‑1 (MQ=1) ccggtcgggCTGGGTGATCTCTCATCAAACAAAATTCAACCCT‑TTTTTTTCTAACCAAACACGCGGCCACCATTTTATACCCCTGTTTCGGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCaa < 2:843140/144‑1 (MQ=1) gggtGGTTGATCTCTCATAAAACAAAATGCACCCCA‑TTGTTTTCTAAGCACTCACGCGGCCACCAGTTGATACCCCTGTTTCGGGGCTTTGCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTg < 2:4657356/147‑1 (MQ=1) gttaTCTCTCATCAAACAAAATCCAACCCA‑TTGTTTTCTAAGCACTCAAGCGGCCACCTTTTGATACCCCAGTTTCCGGGCTTTGCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAAc < 2:5067557/147‑1 (MQ=1) gttaTCTCTCAACAAAAAAAAGGCAACCAT‑TTTTTTTCTACGCACTCACGCGGCCCCCAGTTGATACCCCAGTTTCCGGGATTTGCGAAGAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAAc < 1:4932760/147‑1 (MQ=1) tcccatcaAAAAAATTGCAACCCT‑TTGTTTTCTACGCACTCCCGGGGCCCCCCGTTGATCCCCCTTTTTCGGGGGTTTGCGAAGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGc < 1:2983572/147‑1 (MQ=1) aaaaaccaaaATGCAACCCT‑TTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCGGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGc < 1:5006394/144‑1 (MQ=0) cacacaaaATCCACCCCC‑TTTTTTTCAAAGAACTCACGCGGCCCCCAGTGGATACCCCATTTGCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGctt < 1:4951501/147‑1 (MQ=1) accaaaGGCAACCCC‑TTGTTTTCTAAGCACTCCCGGGGCCACCTGTTGAAACCCCTGTTTCCGGGCTTTGCGAAGCACCGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGcttctt < 2:1153003/147‑1 (MQ=1) aaaaaccAACCCATTTTTTTTCTAACCACACACGCGGCCACCTTTTCAACCCCCTGTTGCGGGGTTTTGCCACCCAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTg < 1:1723362/144‑1 (MQ=0) caccccT‑TTGTTTTCTAAGCACTCAAGCGGCCACCTGTTGATACCCCTGTTGCCGGGCTTGGCGAAGCAATGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa < 2:2646873/147‑1 (MQ=1) caccccT‑TTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa < 2:1318986/147‑1 (MQ=1) caccccA‑TTGTTTTCTACGCACCCACGCGGCCACCCGTTGATCCCCCCGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa < 2:2895060/147‑1 (MQ=1) ccttttttttCTACCCACCCCCGCGCCCACCCGTTGATCCCCCTGTTTCGGGTCTTGGCGACCAAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt < 1:850857/144‑1 (MQ=1) ccttttttttCTACCAACACACGCGGCCACCTTTTTAAACCCCTTTTGCGGGTCTTTGCGAAGCAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt < 1:4124101/144‑1 (MQ=1) cctttGTTTTCTACGCACCCCCGCGGCCCCCCGTTGCTACCCCAGTGGCGGGGCTTTGCGACGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt < 2:1671363/147‑1 (MQ=1) cctttGTTTTCTAAGCACTCCAGGGGCCACCCGTTGAGACCCCTGGTTCCGGGGGTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt < 2:5329324/147‑1 (MQ=1) cctttGTTTTCAACGCACTCACGCGGCCACCTGTTGATACCCCTGTTTCGGGGCTTTGCTACGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt < 2:4026177/147‑1 (MQ=1) cccttGTTTTCTCCGCACCCACGCGGCCACCCGTGGAGACCCCCGTTTCGGGGCTTTGCGAAGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt < 2:3246502/147‑1 (MQ=1) cttttttttCTATCCACACAGGCGCCCACCCGTCCTTCCCCCTTTTGCGGGGCTTGGCGACCAAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTg < 1:4162289/145‑1 (MQ=0) cttttttttCTAAGCACTCACGGGGCCACCAGTTGATACCCCCGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTg < 1:3482870/145‑1 (MQ=0) agaaCTCAAGCGGCCACCTTTTTATACCCCAGTTTCGGGGGTTTGCTAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGt < 2:895083/147‑1 (MQ=1) acacacgcgcCCACCTGTTGAAACCCCTGTTTCGGGGCTTTGCGATGAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGcc < 1:4761147/144‑1 (MQ=1) cacgCGGCCCCCTTTTGACACCCATGTTGCGGGGCTTTCCGAACAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAAc < 1:5099522/147‑1 (MQ=1) tttaTACCCCCGTTTCCGGGCTTTGCGAAGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGc < 2:4562375/147‑1 (MQ=0) atcccccAGTTGCCGGGCTTGGCGACGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCtt < 1:2488571/147‑1 (MQ=0) cccctttttCGGGGTTTTGCTAAGAAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGAt < 1:3055602/144‑1 (MQ=1) cgtctCCGGGCTTTGCGACGCACAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGgcg < 2:4076996/146‑1 (MQ=0) gtgtCCGGGCTTTGCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGgcgc < 1:5375619/147‑1 (MQ=1) | GCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCA‑TTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGGCGC > NC_002947/5222945‑5223191 |
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: ‑ |