Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,857,969 | A→C | 9.7% | T4T (ACT→ACG) ‡ | PP_4272 ← | membrane protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,857,969 | 0 | A | C | 9.7% | 122.3 / 2.9 | 45 | T4T (ACT→ACG) ‡ | PP_4272 | membrane protein |
Reads supporting (aligned to +/- strand): ref base A (19/21); new base C (3/2); total (22/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.65e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AACCGCCTCAGCAGCAACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGA > NC_002947/4857820‑4858116 | aaCCGCCTCAGCAGCAACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGaa < 2:129107/150‑1 (MQ=255) cTCAGCAGCAACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCaaa > 1:917378/1‑149 (MQ=255) cagcaACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCATTCTGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTtca‑a < 1:660051/150‑1 (MQ=255) agcaACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTtca‑ag < 1:19358/150‑1 (MQ=255) aCCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGaa > 2:855627/1‑134 (MQ=255) aCCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGaa < 1:855627/134‑1 (MQ=255) ccACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGc > 2:946408/1‑150 (MQ=255) ccACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGc > 1:377580/1‑150 (MQ=255) cacaGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCAAAAAAGTTGAAGTCAACTACA‑AGTCAAGGGcc > 2:308441/1‑150 (MQ=255) tacATCAGTAAAGTGACACCAAAGCAAATCGGCGATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAATCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAAt < 2:232218/148‑1 (MQ=255) cATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATgg < 1:1292023/150‑1 (MQ=255) cATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATgg < 1:789316/150‑1 (MQ=255) aTCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCGTGCCGCGCAGCATAGAAGTTGAAGTCAGCTTCA‑AGTCAAGGGCCCGAAATGGa > 2:1157412/1‑150 (MQ=255) cACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACa < 1:308441/150‑1 (MQ=255) acctggCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGc > 2:895134/6‑150 (MQ=255) cAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATgg > 2:1125554/1‑133 (MQ=255) cAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATgg < 1:1125554/133‑1 (MQ=255) aaaTGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCCCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGcc < 2:568306/149‑1 (MQ=255) aaTCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCaaa < 2:1205114/146‑1 (MQ=255) aaTCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCaaa > 1:1205114/1‑146 (MQ=255) cGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTg < 1:946408/150‑1 (MQ=255) tGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTAc > 2:403666/1‑150 (MQ=255) cTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCg > 1:1211173/1‑121 (MQ=255) cTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCg < 2:1211173/121‑1 (MQ=255) tAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTtca‑agtcaa > 2:564219/1‑98 (MQ=255) tAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTtca‑agtcaa < 1:564219/98‑1 (MQ=255) cATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGacac < 2:775838/150‑1 (MQ=255) tCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTg > 2:774272/1‑150 (MQ=255) aGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGtt > 2:1244146/1‑150 (MQ=255) gAACAAGCCCTCGAAGGTTGACGGTCATTACCGTGCCGCGCAGCATAGAAGTTGAAGTCAGCTTCACAGT‑AAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGaaa < 1:1157412/150‑1 (MQ=255) gAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGACTTCAACTtctat < 1:1053355/68‑4 (MQ=38) gAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGACTTCAACTtctat > 2:1053355/1‑65 (MQ=38) cAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGc < 1:1244146/150‑1 (MQ=255) cAAGCCCGCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGc < 1:567834/150‑1 (MQ=255) gCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTg > 2:361757/1‑150 (MQ=255) cTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTc < 2:377580/150‑1 (MQ=255) tCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCa > 2:236083/1‑150 (MQ=255) gAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTAAAGGGCCCGAAATGGACATTGCCGACGTCTCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGACGAAAGGCTTGCTCAAg > 2:1366650/1‑150 (MQ=255) tCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCggg > 1:817944/1‑150 (MQ=255) ccGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGCGCCCGGCCaagg < 1:774272/150‑4 (MQ=255) gcgcAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTc < 2:817944/150‑1 (MQ=255) ttgAAGTTGCCTTCAACTTCA‑AGTCAAGGGCCCGAGATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGttt > 2:59822/3‑150 (MQ=255) tgAAGTTGACTTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTg < 2:324015/85‑1 (MQ=255) tgAAGTTGACTTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTg > 1:324015/2‑86 (MQ=255) aGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTg > 1:45114/1‑150 (MQ=255) gTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGcc < 2:450076/121‑1 (MQ=255) gTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGcc > 1:450076/1‑121 (MQ=255) gAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGa < 1:361757/150‑1 (MQ=255) | AACCGCCTCAGCAGCAACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCATAGAAGTTGAAGTCAACTTCA‑AGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGA > NC_002947/4857820‑4858116 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |