| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 1,433,718 | T→C | intergenic (‑44/+273) | pinR ← / ← ynaE | Rac prophage; putative site‑specific recombinase/cold shock protein, Rac prophage |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,433,718 | 0 | T | C | 96.9% | 62.5 / ‑4.1 | 33 | intergenic (‑44/+273) | pinR/ynaE | Rac prophage; putative site‑specific recombinase/cold shock protein, Rac prophage |
| Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (18/13); total (20/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01 | |||||||||||
TCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTC > NC_000913/1433622‑1433766 | tCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAg < 2:1353747/100‑1 (MQ=17) ggtggtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAACTTGGTGTTCATATCATAATGCATCCCAg > 2:1532674/1‑98 (MQ=12) cTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGGATCCCAGAAaataat > 1:1553158/1‑100 (MQ=16) cTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAaataat < 2:620378/100‑1 (MQ=18) ngATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAAtt > 1:956255/2‑100 (MQ=17) gATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTa > 1:748440/1‑100 (MQ=18) aTCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTAt > 2:446075/1‑100 (MQ=18) cAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGat > 1:949753/1‑100 (MQ=18) aGCGTTGATATTCGACGGTAAGCAGGAATTCGGGACATGATGAGATTCTATACGAAATTGTGGTTCATATCATAATGCATCCCAGAAAATATTTATGa < 2:627929/98‑1 (MQ=11) aGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGatt > 2:1523846/1‑100 (MQ=18) aGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGatt > 1:84316/1‑100 (MQ=18) gCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGatta < 2:8071/100‑1 (MQ=18) gTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTAGTGTTCATATCATAATGCATCCCAGAAAATAATTATGattatt < 2:1086060/100‑1 (MQ=255) tGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTAtttt < 1:209406/100‑1 (MQ=17) gATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTAttttt > 2:867899/1‑100 (MQ=18) gATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTAttttt > 1:416499/1‑100 (MQ=18) atatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtg < 1:231298/100‑1 (MQ=18) ccGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCa > 1:1159757/1‑100 (MQ=18) ccGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATAATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCa > 2:1314463/1‑100 (MQ=18) ccGACAGTAAGCAAAAATTCGAGACAAAGTAAGACAATATACGAAATAGATGTTCATATCATAATGCATCACATATAGTAATTATAATTAtttttataca > 2:1340593/1‑95 (MQ=2) gTAAACAAAAATTCGAGACATAGTGAGACTCTATACGGAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtg < 2:294864/100‑1 (MQ=255) tAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgt > 2:1231349/1‑100 (MQ=18) aaGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta < 1:600748/100‑1 (MQ=18) aaGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta < 2:802003/100‑1 (MQ=17) aaGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta < 1:209008/100‑1 (MQ=18) aaGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta < 1:368874/100‑1 (MQ=18) aGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtat > 1:111151/1‑100 (MQ=18) gCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtatg > 2:163825/1‑100 (MQ=18) cAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtatgt < 1:409047/100‑1 (MQ=17) aaaaaTTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtat > 1:424088/1‑97 (MQ=17) ttCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGt > 1:1172072/1‑100 (MQ=17) cGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTc > 1:951130/1‑100 (MQ=11) cGAGACATAGTGAGACTCTATACGAAATTGGCGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTc > 2:1336208/1‑100 (MQ=2) | TCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTC > NC_000913/1433622‑1433766 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |