Sample Resequencing Stats
Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate |
Predicted Mutations |
Mean Coverage |
Total Reads |
Percent Mapped |
Mapped Reads |
Average Read Length |
A15 F177 I0 R1
|
51 |
26.9 |
575130 |
96.6% |
555575 |
140.9 |
Breseq alignment
BRESEQ :: Evidence
|
evidence |
seq id |
position |
mutation |
freq |
annotation |
gene |
description |
RA |
CP000730 |
1,318,237 |
T→A |
100% |
F124Y (TTT→TAT) |
miaB → |
2‑methylthioadenine synthetase |
|
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
* | CP000730 | 1,318,237 | 0 | T | A | 99.5%
| 73.3
/ ‑6.9
| 27 | F124Y (TTT→TAT) | miaB | 2‑methylthioadenine synthetase |
Reads supporting (aligned to +/- strand): ref base T (0/0); major base A (18/8); minor base G (1/0); total (19/8) |
Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
Rejected as polymorphism: E-value score below prediction cutoff. |
Rejected as polymorphism: Frequency below/above cutoff threshold. |
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
TGAATGCACATGACACTGAGGTCATTGCTGGTATACTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTTTAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTAAAATCGTATCAAAATGTAGATATGATATTT > CP000730/1318102‑1318373
|
tgaatgCACATGACACTGAGGTCATTGCTGGTATACTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTg > 2:228737/1‑141 (MQ=255)
atgCACATGACACTGAGGTCATTGCTGGTATACTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGaaa < 1:102244/141‑1 (MQ=255)
atgCACATGACACTGAGGTCATTGCTGGTATACTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGaaa < 1:31088/141‑1 (MQ=255)
gTCATTGCTGGTATACTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCAt > 1:80523/1‑141 (MQ=255)
ataCTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGAGTAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGaa > 2:5825/1‑141 (MQ=255)
ataCTTAAAGCATTAGGCTATCAAGCAACGACAGAAATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGAGAAATAAGTAATTTGAAGCATTTGAGAAAAGaa > 2:219656/1‑141 (MQ=255)
aCTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAACTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACg < 1:5229/141‑1 (MQ=255)
ttGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGGGATTAGAGATAATGCCGAGTACACAGTGTATAGGGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGAc > 2:5874/1‑141 (MQ=255)
tCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAAGTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCgg < 1:196134/141‑1 (MQ=255)
aGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGgtgt < 1:259065/141‑1 (MQ=255)
aaCGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTtg > 2:35854/1‑141 (MQ=255)
aCGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTtgt < 1:161565/141‑1 (MQ=255)
aCGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTtgt > 2:14517/1‑141 (MQ=255)
gATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGt > 1:19288/1‑141 (MQ=255)
ttAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTGGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGt > 2:57467/1‑141 (MQ=255)
aGATGTTATTTTAATCAATACATGTGCGATTAGGGAAAATGCCGAGAACAAAGGGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAagaaga < 2:175430/141‑1 (MQ=255)
gATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGGATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAagag > 1:105458/1‑141 (MQ=255)
taattaatACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTgc > 1:10378/1‑140 (MQ=255)
aattaatACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGaa > 2:18109/1‑141 (MQ=255)
aattaatACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGaa > 2:156545/1‑141 (MQ=255)
attaatACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGCAt > 2:212879/1‑141 (MQ=255)
cATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGATATATGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAAttt > 1:151066/1‑141 (MQ=255)
tgtgCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTa > 2:92360/1‑141 (MQ=255)
agagAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTaaactagaag > 2:140676/1‑134 (MQ=255)
gagaAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTCGTGAATAAAATTTTAAAATCGTATc > 2:98980/1‑141 (MQ=255)
gaAAGTGCGCAGGACAAAGTGTATAGTGAAATAGGTAATTGGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTTTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTAAAATCGTATCaa < 2:12970/141‑1 (MQ=255)
gtgtATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTAAAATCGTATCAAAATGTAGATATGATAttt > 1:78712/1‑141 (MQ=255)
|
TGAATGCACATGACACTGAGGTCATTGCTGGTATACTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTTTAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTAAAATCGTATCAAAATGTAGATATGATATTT > CP000730/1318102‑1318373
|
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
N/A