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 |
|
A1 F2 I208 R1
|
222 |
13.1 |
731276 |
96.6% |
706412 |
85.9 |
Breseq alignment
BRESEQ :: Evidence
|
| evidence |
seq id |
position |
mutation |
annotation |
gene |
description |
| RA |
NZ_CP009273 |
2,797,377 |
(T)5→4 |
coding (517/960 nt) |
nrdF → |
class 1b ribonucleoside‑diphosphate reductase subunit beta |
| |
seq id |
position |
ref |
new |
freq |
score (cons/poly) |
reads |
annotation |
genes |
product |
| * | NZ_CP009273 | 2,797,373 | 0 | T | . | 93.3%
| 53.2
/ ‑1.1
| 15 | coding (513/960 nt) | nrdF | class 1b ribonucleoside‑diphosphate reductase subunit beta |
| Reads supporting (aligned to +/- strand): ref base T (0/1); new base . (3/11); total (3/12) |
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 |
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01 |
ATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATA > NZ_CP009273/2797297‑2797457
|
atCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTA‑TTTTCCAGCCGCgg < 1:36679/90‑1 (MQ=255)
tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACcgc < 1:59931/90‑1 (MQ=255)
tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACcgc < 2:24801/90‑1 (MQ=255)
cGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCCGGCGGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACcgc < 2:329892/90‑1 (MQ=255)
aGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGAcc > 2:171455/1‑90 (MQ=255)
gtTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATc > 2:317777/1‑90 (MQ=255)
tttttGTTCTATTCCGGTTTCTGGCTGCCGATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATcc > 2:133715/1‑90 (MQ=255)
aTTCCGGTTTCTGGCTGCCGATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGc < 1:160675/90‑1 (MQ=255)
aTTCCGGTTTCTGGCTGCCGATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGc < 1:361780/90‑1 (MQ=255)
tttCTGGCTGCCGATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCAc < 1:262486/90‑1 (MQ=255)
ttCTGGCTGCCGATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACg < 1:238689/90‑1 (MQ=255)
tGGCTGCCGATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGtt < 1:133715/90‑1 (MQ=255)
tGGCTGCCGATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGtt < 1:171455/90‑1 (MQ=255)
gCTGCCGATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTtac < 2:46381/90‑1 (MQ=255)
gATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATa < 1:317777/90‑1 (MQ=255)
|
ATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATA > NZ_CP009273/2797297‑2797457
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
GATK/CNVnator alignment
BRESEQ :: bam2aln output
ATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAAT > NZ_CP009273/2797297‑2797467
|
ATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTA‑TTTTCCAGCCGCGGAAAGCTGACC < SRR3722090.37055/100‑1 (MQ=60)
TCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATC < SRR3722090.60578/100‑1 (MQ=60)
ATTCCGGTTTCTGGCTGCCGATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGT < SRR3722090.162296/100‑1 (MQ=60)
ATTCCGGTTTCTGGCTGCCGATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGT < SRR3722090.366506/100‑1 (MQ=60)
TTTCTGGCTGCCGATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACA < SRR3722090.265713/100‑1 (MQ=60)
TTCTGGCTGCCGATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACAT < SRR3722090.241527/100‑1 (MQ=60)
TGGCTGCCGATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGG < SRR3722090.135126/100‑1 (MQ=60)
TGGCTGCCGATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGG < SRR3722090.173248/100‑1 (MQ=60)
agacagaTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAA < SRR3722090.184355/93‑1 (MQ=60)
GATGTA‑TTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAAT < SRR3722090.321900/100‑1 (MQ=60)
|
ATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAAT > NZ_CP009273/2797297‑2797467
|
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |