Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 4,480,143 | T→C | 21.9% | Y9H (TAT→CAT) | yjgL → | hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 4,480,143 | 0 | T | C | 21.9% | 59.3 / 15.1 | 32 | Y9H (TAT→CAT) | yjgL | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (11/14); new base C (5/2); total (16/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.94e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01 |
AAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCG > W3110S.gb/4480074‑4480207 | aaaTTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAatta > 1:144754/1‑71 (MQ=255) tAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCaa < 1:1511811/71‑1 (MQ=255) aGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACacatta < 1:816673/70‑1 (MQ=255) aCGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACacattacat < 1:638640/71‑1 (MQ=255) tCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACacattaca > 1:1064560/1‑66 (MQ=255) tCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACacattaca > 1:154756/1‑66 (MQ=255) tCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGc > 1:262953/1‑71 (MQ=255) tATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGc > 1:45526/1‑68 (MQ=255) tATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGa > 1:1162908/1‑71 (MQ=255) tCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACaa > 1:553174/1‑71 (MQ=255) cACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACaa < 1:59113/70‑1 (MQ=255) cACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACaa < 1:1118207/70‑1 (MQ=255) cACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACaa < 1:164502/70‑1 (MQ=255) cTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATtctc < 1:814456/46‑1 (MQ=255) tCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAatcatt > 1:1206834/1‑42 (MQ=255) gAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGc < 1:390090/56‑1 (MQ=255) gAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGc < 1:979816/56‑1 (MQ=255) gAGTTCGTCATGAGCAAAATATCAGATTTAAATCATTCTCAACACATTACATTAg > 1:1635480/1‑55 (MQ=255) aGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACaa < 1:800544/71‑1 (MQ=255) cGTCATGAGCAAAATATCAGATTTAAATCATTCTCAAcaca > 1:452237/1‑41 (MQ=255) aTGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACa > 1:1728283/1‑62 (MQ=255) aTGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTg > 1:1476327/1‑70 (MQ=255) aTGAGCAAAATATCAGATTTAAATCATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTg > 1:1124716/1‑70 (MQ=255) aTGAGCAAAATATCAGATTTAAATCATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTg > 1:1695191/1‑70 (MQ=255) gCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACa < 1:1020114/58‑1 (MQ=255) cAAAATATCAGATTTAAATCATTCTCAACACATTACATTAGCCGACAATTTTaaa > 1:1326713/1‑55 (MQ=255) atatCAGATTTAAATCATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTaaa < 1:229144/71‑1 (MQ=255) tatCAGATTTAAATCATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAAt < 1:1322298/71‑1 (MQ=255) tCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATAc < 1:680510/71‑1 (MQ=255) tCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATAc < 1:1084862/71‑1 (MQ=255) aGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCt < 1:1737632/71‑1 (MQ=255) ttAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTa > 1:368921/1‑60 (MQ=255) ttAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCg > 1:675489/1‑71 (MQ=255) aaaTTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCt > 1:456271/1‑65 (MQ=255) attattCTCAACACATTACATTAGCCGACAATTTTaa < 1:297541/37‑1 (MQ=255) | AAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCG > W3110S.gb/4480074‑4480207 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |