Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,785,984 | 0 | G | C | 36.0% | 58.3 / 55.8 | 50 | R242P (CGT→CCT) | gpmM | 2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase |
Reads supporting (aligned to +/- strand): ref base G (16/16); new base C (9/9); total (25/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGT > NC_000913/3785906‑3786063 | ccGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCCTg > 2:744151/1‑81 (MQ=255) gATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCt > 1:744133/1‑81 (MQ=255) gATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCCTGCt > 1:744119/1‑81 (MQ=255) tACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGa < 1:744057/81‑1 (MQ=255) ccgTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGt > 1:744130/1‑80 (MQ=255) tGCTGGTTTGCCGTCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGTAAGCGACCGTTATCCGTGCTGAAGGTCAGc > 2:744116/1‑80 (MQ=255) gCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGATAGCGACCGTTATCCGTGCTGAAGGTCAGc > 2:744110/1‑79 (MQ=255) gCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGcc < 1:744058/80‑1 (MQ=255) gCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCCTGCTGAAGGTCAGcc < 2:744059/80‑1 (MQ=255) cTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAg < 1:744060/81‑1 (MQ=255) gTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCCTGCTGAAGGTCAGCCAGAcg < 1:744061/81‑1 (MQ=255) tttGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAcgc < 2:744062/81‑1 (MQ=255) ttGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGAcgcg > 2:744131/1‑81 (MQ=255) aGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCCTGCTGAAGGTCAGCCAGACGCGGCaa < 2:744063/81‑1 (MQ=255) aGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCCTGCTGAAGGTCAGCCAGACGCGGCaa > 2:744147/1‑81 (MQ=255) gctTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATgg > 1:744146/1‑79 (MQ=255) ctTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGaa < 2:744065/80‑1 (MQ=255) ctTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCCTGCTGAAGGTCAGCCAGACGCGGCAATGGAAg < 2:744064/81‑1 (MQ=255) tATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCCTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGAc > 2:744136/1‑81 (MQ=255) aTGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCCTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACg > 1:744156/1‑81 (MQ=255) cTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCg < 2:744066/81‑1 (MQ=255) cTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCCTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGc < 2:744067/80‑1 (MQ=255) gcgACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATg < 1:744068/81‑1 (MQ=255) gacgaAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCa > 1:744108/1‑81 (MQ=255) acgaAAATGATGAATTCGTGAAAGCGACCGTTATCCCTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCAc < 1:744069/81‑1 (MQ=255) gaAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTg < 2:744070/81‑1 (MQ=255) aaaaTGATGAATTCGTGAAAGCGACCGTTATCCCTGCTGAAGGTCAGCCAGACGCAGCAATGGAAGACGGCGATGCACTGa > 1:744139/1‑81 (MQ=255) aaaTGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGAt < 1:744071/81‑1 (MQ=255) aaTGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCGGACGCGGCAATGGAAGACGGCGATGCACTGAtt < 1:744072/81‑1 (MQ=255) aaTGATGAATTCGTGAAAGCGACCGTTATCCCTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGAtt > 1:744157/1‑81 (MQ=255) atgatgAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGAttt < 1:744073/81‑1 (MQ=255) gatgaATTCGTGAAAGCGACCGTTATCCCTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGAtttt > 2:744138/1‑80 (MQ=255) atgaATTCGTGAAAGCGACCGTTATCCCTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTc < 1:744074/80‑1 (MQ=255) gAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCg > 1:744144/1‑81 (MQ=255) aaaGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCg < 2:744075/80‑1 (MQ=255) aGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGc > 2:744174/1‑81 (MQ=255) gCGACCGTTATCCCTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCt < 2:744076/81‑1 (MQ=255) gACCGTTATCCCTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGa > 2:744150/1‑81 (MQ=255) ccGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGAcc < 2:744077/81‑1 (MQ=255) ccGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGAcc < 2:744078/81‑1 (MQ=255) gTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcgc > 2:744175/1‑81 (MQ=255) gTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcgc > 2:744159/1‑81 (MQ=255) gTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcgc > 1:744171/1‑81 (MQ=255) gTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcgc > 1:744126/1‑81 (MQ=255) ttATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcgcg < 1:744080/81‑1 (MQ=255) ttATCCCTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcgcg < 2:744079/81‑1 (MQ=255) tATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcgcgc > 1:744143/1‑81 (MQ=255) tCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcgcgcgc > 1:744166/1‑81 (MQ=255) tCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACcgcgcgc < 1:744081/81‑1 (MQ=255) cGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGt > 1:744141/1‑81 (MQ=255) gTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGt > 2:744142/1‑80 (MQ=255) | CCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGT > NC_000913/3785906‑3786063 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |