Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,606,449 | A→T | 16.1% | intergenic (‑157/+54) | PP_1407 ← / ← PP_t64 | sulfate transporter/tRNA‑Gly |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,606,449 | 0 | A | T | 16.1% | 76.3 / 8.8 | 31 | intergenic (‑157/+54) | PP_1407/PP_t64 | sulfate transporter/tRNA‑Gly |
Reads supporting (aligned to +/- strand): ref base A (12/14); new base T (3/2); total (15/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.54e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01 |
TTTCGAATAGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGACT > NC_002947/1606302‑1606590 | tttCGAATAGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAg > 2:171744/1‑150 (MQ=255) tCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCTTGGCTCCtt < 2:340530/145‑1 (MQ=255) tCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCTTGGCTCCtt > 1:340530/1‑145 (MQ=255) gTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGTCAAGGCTCCTTTTGACTGCGATGACGt > 1:64406/1‑150 (MQ=255) tatCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGa > 1:291450/1‑150 (MQ=255) atCCGCAACAAGCAGTTTCGCTTGCGCCCACAAAGGCAGGACGCGCAGTGAACCGGAGCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGAAGGAAAAAAAAAGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGaa < 2:55751/150‑1 (MQ=255) tCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAAc < 1:171744/150‑1 (MQ=255) tCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAAc < 2:291450/150‑1 (MQ=255) tCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGGGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGAAGCCAAGGCGCCTTTTGACTGCGATGACGTTGCTGGAAc > 1:360289/1‑150 (MQ=255) gcaAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCggg > 1:189327/1‑150 (MQ=255) gCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAgg < 2:466259/150‑1 (MQ=255) cGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGaa < 2:64406/150‑1 (MQ=255) gcATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAAc > 2:219897/1‑150 (MQ=255) aTTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTgg > 2:283805/1‑102 (MQ=255) aTTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTgg < 1:283805/102‑1 (MQ=255) aaGCGGATCGGGGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCTGGTAGAGGGAATCGAACCCTCaa < 2:181045/150‑1 (MQ=255) ggCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTc < 2:189327/150‑1 (MQ=255) gCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTaga > 1:325250/1‑112 (MQ=255) gCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTaga < 2:325250/112‑1 (MQ=255) gCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGaa > 1:434756/1‑110 (MQ=255) gCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGaa < 2:434756/110‑1 (MQ=255) aaGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACt > 1:165753/1‑150 (MQ=255) cAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATAc < 2:165753/150‑1 (MQ=255) aaaaCGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCt < 1:219897/150‑1 (MQ=255) gAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCt > 2:73950/1‑150 (MQ=255) gAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCt < 1:73950/150‑1 (MQ=255) aaaaGGAGCCTTGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCt > 2:479945/1‑142 (MQ=255) aaaaGGAGCCTTGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCt < 1:479945/142‑1 (MQ=255) aaaaGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGc < 2:278071/146‑1 (MQ=255) aaaaGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGc > 1:278071/1‑146 (MQ=255) aaGGAGCCTTGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGACt > 2:406469/1‑150 (MQ=255) | TTTCGAATAGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGACT > NC_002947/1606302‑1606590 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |