Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,665,050 | C→G | 14.4% | Y46* (TAC→TAG) | nuoI → | NADH‑quinone oxidoreductase subunit I |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,665,050 | 0 | C | G | 14.4% | 80.3 / 9.1 | 35 | Y46* (TAC→TAG) | nuoI | NADH‑quinone oxidoreductase subunit I |
Reads supporting (aligned to +/- strand): ref base C (16/14); new base G (3/2); total (19/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-01 |
AGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGCTGGTACCCGGAGTTCTTCCGCA > NC_002947/4664901‑4665198 | aGGATTTGACCCATGTTCACGTATATCGGCGACATCGTTAAGAGCACCGGCACCCAGCGGCGGCGCCTGGCCCCGGTGTTCTCCTGCAGGCTCCGCAAGCGCGACACCCTGCAATACCCCGCAG‑ACCGGTAGTACCTGCTTCCGCGCTAc > 2:281831/1‑150 (MQ=255) aCCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGGTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCg > 1:99994/1‑150 (MQ=255) gCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGAcc > 1:275970/1‑150 (MQ=255) gCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGAcc > 1:565590/1‑150 (MQ=255) gggCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCgagga > 1:606990/1‑150 (MQ=255) gcaccggcaccCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGGTAGCGCGGCGGCag > 2:375033/1‑117 (MQ=255) gcaccggcaccCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGGTAGCGCGGCGGCag < 1:375033/118‑2 (MQ=255) gcaccggcaccCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAgc > 2:327208/1‑150 (MQ=255) ccggcaccCCGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCt < 1:694458/150‑1 (MQ=255) cccgCTGCGCAGCCTGGCATTTGTGTTCTCCCAGGGTTTTTGTAAGCGGGACACCCTGTAATACCCCGAAGAACCCGT‑GTACCTGCCGTCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGATCCCGATGGCGAGGAGCGCTGCGTAGc < 2:177579/147‑1 (MQ=255) ggTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGc < 1:412585/150‑1 (MQ=255) ctcCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATTGGCAGTAGCTCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGccc < 1:234007/150‑1 (MQ=255) ccACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGt > 2:507534/1‑150 (MQ=255) ccACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGt > 2:291750/1‑150 (MQ=255) cACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCcgc > 1:492642/1‑90 (MQ=255) cACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTg > 2:784564/1‑150 (MQ=255) cACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGT‑GTACATGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCcgc < 2:492642/90‑1 (MQ=255) gggTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGc < 1:817575/150‑1 (MQ=255) ttCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGAt > 1:843632/1‑93 (MQ=255) ttCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGAt < 2:843632/93‑1 (MQ=255) gcgACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGc < 2:521089/150‑1 (MQ=255) gACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACctct < 2:761728/112‑1 (MQ=255) gACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACctct > 1:761728/1‑112 (MQ=255) cacCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCgg < 1:797350/116‑1 (MQ=255) cacCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCgg > 2:797350/1‑116 (MQ=255) cccGT‑GTACCTTCCGCCGCGCTACCGCGGCCGCAGCGTCCTTACCAGCCGCCACCATCCCGTGGGGCGCTCCGTAGCGTACCGCCTCCGCTCGCTGCCCTGCCAGGTGGGTTGCATTTCGCTGCCGAAGGCCGAGGGAGAGGCCGGCCGc > 2:144070/1‑150 (MQ=255) ccGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGCt > 1:217591/1‑150 (MQ=255) tACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGCTGGTAc < 2:606990/150‑1 (MQ=255) gatgcggccgcGGTAGCGCGGCGGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGCTGGTAccc > 2:86503/7‑150 (MQ=255) gatGCGGCCGCGGTAGCGCGGCGGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGCTGGTAccc > 1:498689/3‑150 (MQ=255) cTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTTCAACCCCTGAACACCAGCTTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGGCGGTGGTACCCg > 1:488153/1‑150 (MQ=255) gcggccgcGGTAGCGCGGCGGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGCTGGTACCCGGa < 2:498689/147‑1 (MQ=255) ccgccgCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCACTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGCTGGTACCCGGAg < 1:281831/150‑1 (MQ=255) gccgcGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGCTGGTACCCGGAGtt < 1:784564/150‑1 (MQ=255) gcgcTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCACTTCGCTGCAGCAGGTAGAGACAGAGGGCGGCCGCTGGTATCCGGAGttctt > 2:520893/1‑150 (MQ=255) aCCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGCTGGTACCCGGAGTTCTTCCGCa < 2:696017/150‑1 (MQ=255) | AGGATTTGACCCATGTTCAAGTATATCGGCGACATCGTTAAGGGCACCGGCACCCAGCTGCGCAGCCTGGCAATGGTGTTCTCCCACGGGTTCCGCAAGCGCGACACCCTGCAATACCCCGAAGAACCCGT‑GTACCTGCCGCCGCGCTACCGCGGCCGCATCGTCCTCACCCGCGACCCCGATGGCGAGGAGCGCTGCGTAGCGTGCAACCTCTGCGCGGTGGCCTGCCCGGTGGGCTGCATTTCGCTGCAGAAGGCCGAGACAGAGGACGGCCGCTGGTACCCGGAGTTCTTCCGCA > NC_002947/4664901‑4665198 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |