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 I227 R1 223 21.3 1167702 97.0% 1132670 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,848,051 T→G F257C (TTT→TGT)  emrD → multidrug efflux MFS transporter EmrD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,848,0510TG96.2% 80.4 / ‑3.4 26F257C (TTT→TGT) emrDmultidrug efflux MFS transporter EmrD
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (21/4);  total (22/4)
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.70e-01

CTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGC  >  NZ_CP009273/3847965‑3848136
                                                                                      |                                                                                     
cTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGtgttt                                                                                    <  1:132040/90‑1 (MQ=255)
  gCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTg                                                                                  >  2:497914/1‑90 (MQ=255)
                                   tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                 >  1:110080/1‑90 (MQ=255)
                                   tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                 >  2:400705/1‑90 (MQ=255)
                                   tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                 >  2:314090/1‑90 (MQ=255)
                                   tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                 >  2:268166/1‑90 (MQ=255)
                                   tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                 >  2:223910/1‑90 (MQ=255)
                                   tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                 >  1:583542/1‑90 (MQ=255)
                                   tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                 >  1:451407/1‑90 (MQ=255)
                                   tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                 >  1:383932/1‑90 (MQ=255)
                                   tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                 >  1:167500/1‑90 (MQ=255)
                                   tGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCt                                                 >  1:248692/1‑90 (MQ=255)
                                         gTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTtctc                                             >  2:90520/1‑88 (MQ=255)
                                         gTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTtctc                                             <  1:90520/88‑1 (MQ=255)
                                                   cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                 >  1:408224/1‑90 (MQ=255)
                                                   cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                 >  2:29290/1‑90 (MQ=255)
                                                   cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                 >  1:81026/1‑90 (MQ=255)
                                                   cAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAAtgtg                                 >  1:361410/1‑90 (MQ=255)
                                                         tttGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTc                           <  1:268166/90‑1 (MQ=255)
                                                         tttGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTc                           <  2:451407/90‑1 (MQ=255)
                                                               tATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTAt                     >  2:248766/1‑90 (MQ=255)
                                                               tATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTAt                     >  1:240560/1‑90 (MQ=255)
                                                                     gCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATctgctg               <  1:42593/90‑1 (MQ=255)
                                                                      ccGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCtgctgc              >  2:319041/1‑90 (MQ=255)
                                                                               gCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTgctggct     >  1:116779/1‑90 (MQ=255)
                                                                                 aGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGctggctgg   >  1:310138/1‑90 (MQ=255)
                                                                                  gCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCtggctggc  >  1:101181/1‑90 (MQ=255)
                                                                                      |                                                                                     
CTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGC  >  NZ_CP009273/3847965‑3848136

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ 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
CTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGG  >  NZ_CP009273/3847965‑3848148
                                                                                      |                                                                                                 
CTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGG                                                                                      <  SRR3722111.133699/100‑1 (MQ=60)
                         GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                             >  SRR3722111.111463/1‑100 (MQ=60)
                         GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                             >  SRR3722111.592693/1‑100 (MQ=60)
                         GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                             >  SRR3722111.169600/1‑100 (MQ=60)
                         GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                             >  SRR3722111.458255/1‑100 (MQ=60)
                         GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                             >  SRR3722111.389450/1‑100 (MQ=60)
                         GTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCT                                                             >  SRR3722111.251616/1‑100 (MQ=60)
                                       CAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATG                                               <  SRR3722111.91639/100‑1 (MQ=60)
                                         GTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTG                                             >  SRR3722111.366441/1‑100 (MQ=60)
                                         GTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTG                                             >  SRR3722111.414202/1‑100 (MQ=60)
                                         GTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTG                                             >  SRR3722111.82014/1‑100 (MQ=60)
                                                     GTATTTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTAT                                 >  SRR3722111.243428/1‑100 (MQ=60)
                                                         TTTGTTTATTCTGCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGC                             <  SRR3722111.271402/100‑1 (MQ=60)
                                                                     GCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCT                 <  SRR3722111.43114/100‑1 (MQ=60)
                                                                     GCCGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCT                 >  SRR3722111.118251/1‑100 (MQ=60)
                                                                       GGATTCCGGCAGCGTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGG               >  SRR3722111.314200/1‑100 (MQ=60)
                                                                        GATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGC              >  SRR3722111.102439/1‑100 (MQ=60)
                                                                                    GTGTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGG  >  SRR3722111.188559/1‑100 (MQ=60)
                                                                                      |                                                                                                 
CTGCTCCGGCGTGCTGATGGGCGCGGTGTTAGGGCTGAGCAGTATGACGGTCAGTATTTTGTTTATTCTGCCGATTCCGGCAGCGTTTTTTGGCGCATGGTTTGCCGGACGTCCCAATAAACGCTTCTCCACGTTAATGTGGCAGTCGGTTATCTGCTGCCTGCTGGCTGGCTTGCTGATGTGG  >  NZ_CP009273/3847965‑3848148

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: