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 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,475,701 T→C S352P (TCG→CCG)  yheS → ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,475,7010TC100.0% 85.8 / NA 26S352P (TCG→CCG) yheSABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/12);  total (14/12)

ATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGC  >  NZ_CP009273/3475618‑3475781
                                                                                   |                                                                                
aTCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAt                                                                            >  1:1055273/1‑90 (MQ=255)
 tCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAtt                                                                           <  1:1107370/90‑1 (MQ=255)
 tCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAtt                                                                           <  2:135055/90‑1 (MQ=255)
       tattCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCaa                                                                     <  2:760573/90‑1 (MQ=255)
       tattCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCaa                                                                     <  1:1205947/90‑1 (MQ=255)
          tctcGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACt                                                                  >  1:614914/1‑90 (MQ=255)
                              aCCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTgcgc                                              >  1:448693/1‑90 (MQ=255)
                              aCCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTgcgc                                              >  2:607550/1‑90 (MQ=255)
                                  ggTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAAACAAACTGTTAGCCGGGGAACGTGCGCCAGt                                          >  1:10019/1‑90 (MQ=255)
                                   gTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTc                                         <  2:695283/90‑1 (MQ=255)
                                   gTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTc                                         >  2:63557/1‑90 (MQ=255)
                                     gCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAg                                       <  2:606127/90‑1 (MQ=255)
                                        cGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCgg                                    >  1:174776/1‑90 (MQ=255)
                                         ggCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGt                                   <  1:463130/90‑1 (MQ=255)
                                         ggCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGt                                   <  2:1213114/90‑1 (MQ=255)
                                         ggCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGt                                   <  2:1171612/90‑1 (MQ=255)
                                              gcgTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAAt                              <  2:30262/90‑1 (MQ=255)
                                                                  gCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAAt                              <  2:1138798/70‑1 (MQ=255)
                                                                  gCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAAt                              >  1:1138798/1‑70 (MQ=255)
                                                                  gCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGCCTGGCGAAAGGGATCa          >  1:207114/1‑90 (MQ=255)
                                                                    aaTGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCaaa        >  1:711009/1‑90 (MQ=255)
                                                                    aaTGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCaaa        >  1:1058443/1‑90 (MQ=255)
                                                                    aaTGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCaaa        >  2:948292/1‑90 (MQ=255)
                                                                        gcgcGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCg    >  1:674094/1‑90 (MQ=255)
                                                                          gcgGGTAAACGGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGc  >  2:946434/1‑90 (MQ=255)
                                                                          gcgGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGc  <  2:674094/90‑1 (MQ=255)
                                                                                   |                                                                                
ATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGC  >  NZ_CP009273/3475618‑3475781

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 31 ≤ 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
GGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCG  >  NZ_CP009273/3475608‑3475779
                                                                                             |                                                                              
GGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAT                                                                          >  SRR3722089.1071111/1‑100 (MQ=60)
     TGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATC                                                                     <  SRR3722089.1247447/100‑1 (MQ=60)
          ATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACT                                                                >  SRR3722089.622593/1‑100 (MQ=60)
           TCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTG                                                               <  SRR3722089.1124114/100‑1 (MQ=60)
                 agatctacactaagatggtcgtcggcagcgtcagatgtgtataagagacagGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCC                                                         <  SRR3722089.1695/49‑1 (MQ=60)
                 TATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCC                                                         <  SRR3722089.1224284/100‑1 (MQ=60)
                              ATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGC                                            >  SRR3722089.454077/1‑100 (MQ=60)
                                  AACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAAACAAACTGTTAGCCGGGGAACGTGCGCCAGT                                        >  SRR3722089.10138/1‑100 (MQ=60)
                                        ACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGG                                  >  SRR3722089.176962/1‑100 (MQ=60)
                                                   GGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTC                       <  SRR3722089.468657/100‑1 (MQ=60)
                                                                  CTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGctgtctctta        >  SRR3722089.1156121/1‑90 (MQ=60)
                                                                  CTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGCCTGGCGAAAGGGATCA        >  SRR3722089.209718/1‑100 (MQ=60)
                                                                    GTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAA      >  SRR3722089.1074326/1‑100 (MQ=60)
                                                                    GTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAA      >  SRR3722089.720263/1‑100 (MQ=60)
                                                                        GGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCG  >  SRR3722089.682826/1‑100 (MQ=60)
                                                                                             |                                                                              
GGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCG  >  NZ_CP009273/3475608‑3475779

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