Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,497,350 | T→C | 10.5% | E215G (GAA→GGA) | PP_5651 ← | hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,497,350 | 0 | T | C | 10.5% | 123.9 / 3.8 | 46 | E215G (GAA→GGA) | PP_5651 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (18/23); new base C (3/2); total (21/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.48e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01 |
TAAGAGTTTCTTGACTTAACTTCATTGGATGAACGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGTATAAAGACTATTCCATGTCCCGGGCTTATATCTTGCAAAGGAAGCACACAGGTCAAT > NC_002947/4497204‑4497492 | tAAGAGTTTCTTGACTTAACTTCATTGGATGAACGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCt > 2:1396411/1‑150 (MQ=255) aaGAGTTTCTTGACTTAACTTCATTGGATGAACGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCtt < 1:1396411/150‑1 (MQ=255) gagTTTCTTGACTTAACTTCATTGGATGAACGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTAt < 1:1073227/150‑1 (MQ=255) agTTTCTTGACTTAACTTCATTGGATGAACGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTAtt < 2:1071245/150‑1 (MQ=255) tttCTTGACTTAACTTCATTGGATGAACGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTgg > 1:891554/1‑150 (MQ=255) aCTTAACTTCATTGGATGAACGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTCCCTTATTGGTCCGGGt > 1:1049163/1‑150 (MQ=255) aCTTAACTTCATTGGATGAACGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTCCCTTATTGGTCCGGGt > 1:339272/1‑150 (MQ=255) aTTGGATGAACGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGaaa > 2:772625/1‑150 (MQ=255) aaCGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCa > 1:468558/1‑150 (MQ=255) aCGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCTGGTCgat > 1:1354326/1‑133 (MQ=255) aCGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCTGGTCgat < 2:1354326/135‑3 (MQ=255) cGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGCCTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAg < 1:549943/150‑1 (MQ=255) tCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTCCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGa < 2:1049163/150‑1 (MQ=255) tCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTCCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGa < 2:339272/150‑1 (MQ=255) acacCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGaaaa < 1:1412285/127‑1 (MQ=255) acacCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGaaaa > 2:1412285/1‑127 (MQ=255) ccATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTAt < 1:895050/102‑1 (MQ=255) ccATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTAt > 2:895050/1‑102 (MQ=255) ccATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTCCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAgg > 2:853469/1‑150 (MQ=255) aCGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTg < 1:227278/150‑1 (MQ=255) ccAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCggg < 2:367418/103‑1 (MQ=255) ccAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCggg > 1:367418/1‑103 (MQ=255) tCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGca < 1:1208749/147‑1 (MQ=255) tCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGca > 2:1208749/1‑147 (MQ=255) tCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACaa > 1:23896/1‑150 (MQ=255) ttCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCt < 2:468558/150‑1 (MQ=255) tCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACa < 2:1226052/141‑1 (MQ=255) tCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACa > 1:1226052/1‑141 (MQ=255) aacaaATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAAtt > 2:431784/1‑150 (MQ=255) aTCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGt < 2:891554/150‑1 (MQ=255) cAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTcc < 2:1322547/47‑1 (MQ=255) cAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTcc > 1:1322547/1‑47 (MQ=255) aTTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGTATAAAGACTATTCCATGTCCCGGGCt < 1:908539/150‑1 (MQ=255) ttCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGTATAAAGGCTATTCCATGTCCCGGGCTta < 1:853469/150‑1 (MQ=255) tCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGTATAAAGACTATTCCATGTCCCGGGCTtat < 1:431784/150‑1 (MQ=255) aTTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGTATAAAGACTATTCCATGTCCCGGGCTTatat > 2:1384114/1‑150 (MQ=255) aTTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGTATAAAGACTATTCCATGTCCCGGGCTTatat > 2:752313/1‑150 (MQ=255) aaCTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGTATAAAGACTATTCCATGTCCCGGGCTTATATCtt < 1:1384114/150‑1 (MQ=255) tGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGt > 2:598964/1‑88 (MQ=255) tGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGt < 1:598964/88‑1 (MQ=255) tGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGtata > 1:1009303/1‑116 (MQ=255) tGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGtata < 2:1009303/116‑1 (MQ=255) aaaaaaGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGTATAAAGACTATTCCATGTCCCGGGCTTATATCTTGCAAAg < 1:772625/150‑1 (MQ=255) aaaaaGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCa > 2:289390/1‑107 (MQ=255) aaaaaGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCa < 1:289390/107‑1 (MQ=255) gACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGTATAAAGACTATTCCATGTCCCGGGCTTATATCTTGCAAAGGAAGCACACAGGTCAAt < 2:23896/150‑1 (MQ=255) | TAAGAGTTTCTTGACTTAACTTCATTGGATGAACGAGGGCGGAATCCACACCCATACGCCCAGTCATGAATTCGATCAACAAATCTGGCGATACGGCCTCACCAACCGATTTCATTAACTGTAAAAAAGCATCAGCATCGACTTTTTCCTTATTGGTCCGGGTCGTTATGAAAAACAACCAAGGCTCCGTCCGACAAAAGGTCTGGTGCACAAACAGTCTGTGAATTGCAGTATAAAGACTATTCCATGTCCCGGGCTTATATCTTGCAAAGGAAGCACACAGGTCAAT > NC_002947/4497204‑4497492 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |