| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,669,686 | A→G | 12.2% | M38V (ATG→GTG) ‡ | metF → | 5,10‑methylenetetrahydrofolate reductase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 5,669,686 | 0 | A | G | 12.2% | 103.8 / 6.0 | 41 | M38V (ATG→GTG) ‡ | metF | 5,10‑methylenetetrahydrofolate reductase |
| Reads supporting (aligned to +/- strand): ref base A (17/19); new base G (2/3); total (19/22) | |||||||||||
| 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.77e-01 | |||||||||||
GCCCAAGCCAGACCGGCCATCACCGATCTGGCTTTTACTTGCAGCTTGCCGCCAGACGCTTGTTGCCAGAGGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGCCGGTGGGTCGACCCGCGACCGCACGCTCAATACCGTGCTGCAGCTGGAAAGCGAAGTGAAGATACCTG > NC_002947/5669537‑5669815 | gCCCAAGCCAGACCGGCCATCACCGATCTGGCTTTTACTTGCAGCTTGCCGCCAGACGCTTGTTGCCAGAGGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGa < 1:834985/150‑1 (MQ=255) cccAAGCCAGACCGGCCATCACCGATCTGGCTTTTACTTGCAGCTTGCCGCCAGACGCTTGTTGCCAGAGGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGAt < 1:585492/150‑1 (MQ=255) ccAAGCCAGACCGGCCATCACCGATCTGGCTTTTACTTGCAGCTTGCCGCCAGACGCTTGTTGCCAGAGGAACGCCTGATGTCACGAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCCGTACCGATGCCGGTCACGAAAAGCTGATg > 1:217856/1‑150 (MQ=255) ccGGCCATCACCGATCTGGCTTTTACTTGCAGCTTGCCGCCAGACGCTTGTTGCCAGAGGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTcgccc < 2:54912/150‑1 (MQ=255) cATCACCGATCTGGCTTTTACTTGCAGCTTGCCGCCAGACGCTTGTTGCCAGAGGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGGCGGGCGAcgcccatcag < 1:76151/150‑6 (MQ=255) cATCACCGATCTGGCTTTTACTTGCAGCTTGCCGCCAGACGCTTGTTGCCAGAGGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGGCGGGCGAcgcccatcag > 2:76151/1‑145 (MQ=255) cATCACCGATCTGGCTTTTACTTGCAGCTTGCCGCCAGACGCTTGTTGCCAGAGGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGGCGGGCGAcgcccatcag < 2:661814/150‑6 (MQ=255) tCACCGATCTGGCTTTTACTTGCAGCTTGCCGCCAGACGCTTGTTGCCAGAGGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCt < 2:513301/150‑1 (MQ=255) aCCGATCTGGCTTTTACTTGCAGCTTGCCGCCAGACGCTTGTTGCCAGAGGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTgg > 1:782241/1‑150 (MQ=255) ccGATCTGGCTTTTACTTGCAGCTTGCCGCCAGACGCTTGTTGCCAGAGGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGc < 1:537458/150‑1 (MQ=255) gCTTTTACTTGCAGCTTGCCGCCAGACGCTTGTTGCCAGAGGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCtata > 1:100925/1‑150 (MQ=255) cgccAGACGCTTGTTGCCAGAGGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTg > 2:814476/1‑150 (MQ=255) aCGCTTGTTGCCAGAGGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTa > 1:422543/1‑150 (MQ=255) gCTTGTTGCCAGAGGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACg > 2:782364/1‑150 (MQ=255) gCTTGTTGCCAGAGGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACg > 2:455402/1‑150 (MQ=255) gCTTGTTGCCAGAGGAACACTTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACg < 1:683394/150‑1 (MQ=255) ttgCCAGAGGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGGCGGGCGAcgcccatca > 2:607709/1‑100 (MQ=255) ttgCCAGAGGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGGCGGGCGAcgcccatca < 1:607709/104‑5 (MQ=255) cAGAGGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGCCGGTg > 1:786701/1‑150 (MQ=255) cAGAGGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACCAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCTGACTTCTTCTCCTGCACCTACGGTGCCGGTg > 2:289529/1‑150 (MQ=255) agagGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGCCGGTgg < 1:191053/150‑1 (MQ=255) gagGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGAACGATCCCGGTCACGTAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGCCGGTggg < 1:103553/150‑1 (MQ=255) aCGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGCCGGTGGGTCGAc < 2:786701/150‑1 (MQ=255) cTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTACCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGCCGGTGGGTCGACCcgc < 1:696893/150‑1 (MQ=255) acaAGAACGCCGCTGCAGTTTCGAGTTTTTCCGGACCAAGAGCGAAGCCTGTCACGAAAAGCTGATGGGCGTGGACCGCGAGCTGGCCACCTATAAGACGGACTTCTTCGTCGGAACCTATGTGGCTGGTGGGTGGACCAGTGAcaggag > 2:339445/1‑145 (MQ=255) acaAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGCCGGTGGGTCGACCCGCGACCGCAc > 1:156010/1‑150 (MQ=255) gAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGc < 1:780908/122‑1 (MQ=255) gAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGc > 2:780908/1‑122 (MQ=255) cTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGCCGGTGGGTCGACCCGCGACCGCACGCTCAATACCGt < 1:782364/150‑1 (MQ=255) cTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGCCGGTGGGTCGACCCGCGACCGCACGCTCAATACCGt > 2:647245/1‑150 (MQ=255) gTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGc > 2:276810/1‑102 (MQ=255) gTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGc < 1:276810/102‑1 (MQ=255) tcttcCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGCCGGTGGGTCGACCCGCGACCGCACGCTCAATACCGTGCTGCAGCTGGaaa < 2:782241/150‑1 (MQ=255) tcttcCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGCCGGTGGGTCGACCCGCGACAGCACGCTCCATACCGTTCTGCAGCTGGaaa < 2:100925/150‑1 (MQ=255) ccGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTGTTCTCCTGCACCTACGGAGCCGGTGGCAGGAGTCGGGACCGCACGCTCAATACCGTGCTGCAGCTGGAACGCGaa > 1:828645/1‑150 (MQ=255) aGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGCCGGTGGGTCGACCCGCGACCGCACGCTCAATACCGTGCTGCAGCTGGAAAGCGAAGTGAAGa < 1:455402/150‑1 (MQ=255) gACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCt > 2:590013/1‑74 (MQ=255) gACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCt < 1:590013/74‑1 (MQ=255) aTGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGCCGGTGGGTCGACCCGCGACCGCACGCTCAATACCGTGCTGCAGCTGGAAAGCGAAGTGAAGATACCTg < 1:289529/150‑1 (MQ=255) aTGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGCCGGTGGGTCGACCCGCGACCGCACGCTCAATACCGTGCTGCAGCTGGAAAGCGAAGTGAAGATACCTg < 2:422543/150‑1 (MQ=255) aCGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGCCGGTGGGTCGACCCGCGACCGCACGCTCAATACCGtg < 1:360138/110‑1 (MQ=255) aCGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGCCGGTGGGTCGACCCGCGACCGCACGCTCAATACCGtg > 2:360138/1‑110 (MQ=255) | GCCCAAGCCAGACCGGCCATCACCGATCTGGCTTTTACTTGCAGCTTGCCGCCAGACGCTTGTTGCCAGAGGAACGCCTGATGTCACAAGAACGCCGCTACAGTTTCGAGTTCTTCCCGACCAAGACCGATGCCGGTCACGAAAAGCTGATGGGCGTCGCCCGCCAGCTGGCCACCTATAACCCGGACTTCTTCTCCTGCACCTACGGTGCCGGTGGGTCGACCCGCGACCGCACGCTCAATACCGTGCTGCAGCTGGAAAGCGAAGTGAAGATACCTG > NC_002947/5669537‑5669815 |
| 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: ‑ |