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 I206 R1 206 23.7 1314362 96.8% 1272302 86.2

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,434,529 A→C F186V (TTT→GTT)  fabB ← beta‑ketoacyl‑ACP synthase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,434,5290AC100.0% 78.6 / NA 25F186V (TTT→GTT) fabBbeta‑ketoacyl‑ACP synthase I
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (10/15);  total (10/15)

CGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGATGGAGT  >  NZ_CP009273/2434448‑2434615
                                                                                 |                                                                                      
cGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGt                                                                                <  1:291661/90‑1 (MQ=255)
           gtagACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGccc                                                                     >  1:488603/1‑90 (MQ=255)
               aCAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGtt                                                                 <  2:512487/90‑1 (MQ=255)
                  gcgcACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGa                                                              >  2:87279/1‑90 (MQ=255)
                  gcgcACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGa                                                              >  1:264730/1‑90 (MQ=255)
                    gcACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATc                                                            <  1:19531/90‑1 (MQ=255)
                    gcACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATc                                                            <  1:206414/90‑1 (MQ=255)
                             gCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACt                                                   >  1:329586/1‑90 (MQ=255)
                               gTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATc                                                            <  1:500727/79‑1 (MQ=255)
                               gTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATc                                                            >  2:500727/1‑79 (MQ=255)
                                tCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCg                                                <  2:19705/90‑1 (MQ=255)
                                 cGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGt                                               >  1:469208/1‑90 (MQ=255)
                                  gAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGtt                                              <  1:586332/90‑1 (MQ=255)
                                          aaGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATAc                                      <  1:151251/90‑1 (MQ=255)
                                          aaGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATAc                                      <  2:488603/90‑1 (MQ=255)
                                            gCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAg                                    <  1:318265/90‑1 (MQ=255)
                                            gCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAg                                    <  2:532190/90‑1 (MQ=255)
                                                           aGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGca                     <  2:286164/90‑1 (MQ=255)
                                                            gCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCGGTTGGATCTGCTCTACCGCGTTACCGACACAGTGTGCGGAAGTCCcac                    >  2:598170/1‑90 (MQ=255)
                                                            gCTCTTCGCCGCCGCCAGCAACCACGAAGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGcac                    >  2:258461/1‑90 (MQ=255)
                                                                      gccgccAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGtt                                              <  2:643492/54‑1 (MQ=255)
                                                                      gccgccAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGtt                                              >  1:643492/1‑54 (MQ=255)
                                                                      gccgccAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGtt                                              <  2:19504/54‑1 (MQ=255)
                                                                      gccgccAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGtt                                              >  1:19504/1‑54 (MQ=255)
                                                                              cAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGATGGagt  <  1:92295/90‑1 (MQ=255)
                                                                                 |                                                                                      
CGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGATGGAGT  >  NZ_CP009273/2434448‑2434615

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ 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
CGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGATGGAGTAGTTAACGCC  >  NZ_CP009273/2434448‑2434625
                                                                                 |                                                                                                
CGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCC                                                                                <  SRR3722088.294450/100‑1 (MQ=60)
 GTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCC                                                                               >  SRR3722088.494424/1‑100 (MQ=60)
        TTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGA                                                                        >  SRR3722088.267295/1‑100 (MQ=60)
                   CGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACT                                                             >  SRR3722088.332987/1‑100 (MQ=60)
                    GCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTG                                                            <  SRR3722088.19738/100‑1 (MQ=60)
                    GCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTG                                                            <  SRR3722088.208478/100‑1 (MQ=60)
                    GCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTG                                                            <  SRR3722088.506735/100‑1 (MQ=60)
                       CCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGT                                                         >  SRR3722088.474728/1‑100 (MQ=60)
                                  GAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAG                                              <  SRR3722088.593518/100‑1 (MQ=60)
                                          AAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGA                                      <  SRR3722088.152787/100‑1 (MQ=60)
                                            GCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAG                                    <  SRR3722088.321490/100‑1 (MQ=60)
                                                            GCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGctgtctcttatacacatctccgagcc                    >  SRR3722088.19710/1‑74 (MQ=60)
                                                            GCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGctgtctcttatacacatctccgagcc                    >  SRR3722088.651462/1‑74 (MQ=60)
                                                                              CAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGATGGAGTAGTTAACGCC  <  SRR3722088.93229/100‑1 (MQ=60)
                                                                                 |                                                                                                
CGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGATGGAGTAGTTAACGCC  >  NZ_CP009273/2434448‑2434625

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: