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 I225 R1 227 21.4 1179702 97.1% 1145490 86.5

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% 86.0 / NA 27F186V (TTT→GTT) fabBbeta‑ketoacyl‑ACP synthase I
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (10/17);  total (10/17)

ATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGATGGAG  >  NZ_CP009273/2434454‑2434614
                                                                           |                                                                                     
aTTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGtt                                                                         >  2:483101/1‑90 (MQ=255)
   tagtagACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGc                                                                      >  1:491439/1‑90 (MQ=255)
    agtagACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGcc                                                                     <  2:457748/90‑1 (MQ=255)
        gACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTc                                                                              <  1:134985/77‑1 (MQ=255)
        gACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTc                                                                              >  2:134985/1‑77 (MQ=255)
          cAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTg                                                               >  1:526102/1‑90 (MQ=255)
            gcgcACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGa                                                             >  1:425705/1‑90 (MQ=255)
            gcgcACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGa                                                             <  2:395234/90‑1 (MQ=255)
                 cccATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGc                                                        <  2:151818/90‑1 (MQ=255)
                           cGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGt                                              >  2:55241/1‑90 (MQ=255)
                            gAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGtt                                             <  2:526102/90‑1 (MQ=255)
                            gAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGtt                                             <  2:491439/90‑1 (MQ=255)
                            gAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGtt                                             <  1:478301/90‑1 (MQ=255)
                            gAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGtt                                             <  2:110364/90‑1 (MQ=255)
                            gAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGtt                                             <  2:218627/90‑1 (MQ=255)
                                    aaGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATAc                                     <  1:363263/90‑1 (MQ=255)
                                    aaGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATAc                                     <  2:319974/90‑1 (MQ=255)
                                      gCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAg                                   <  1:343923/90‑1 (MQ=255)
                                      gCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAg                                   <  1:301497/90‑1 (MQ=255)
                                                      gCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGcc                                                                     <  2:559389/40‑1 (MQ=255)
                                                      gCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGcc                                                                     >  1:559389/1‑40 (MQ=38)
                                                      gCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGcac                   >  1:531498/1‑90 (MQ=255)
                                                      gCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGcac                   >  1:361064/1‑90 (MQ=255)
                                                           tCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACAcgcg              <  1:564273/90‑1 (MQ=255)
                                                           tCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACAcgcg              >  1:380591/1‑90 (MQ=255)
                                                                     cAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGATgg    <  1:202779/90‑1 (MQ=255)
                                                                       gCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGATGGag  <  1:124752/90‑1 (MQ=255)
                                                                           |                                                                                     
ATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGATGGAG  >  NZ_CP009273/2434454‑2434614

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ 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
AGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGATGGAGTAGTTAACGC  >  NZ_CP009273/2434430‑2434624
                                                                                                   |                                                                                               
AGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAC                                                                                                 <  SRR3722109.27073/100‑1 (MQ=60)
 GGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACC                                                                                                >  SRR3722109.569762/1‑100 (MQ=60)
                 TCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGC                                                                                >  SRR3722109.498849/1‑100 (MQ=60)
                   cagGTATTTAGTAGACAGCGCCCCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCC                                                                              <  SRR3722109.136652/97‑1 (MQ=60)
                        ATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTG                                                                         >  SRR3722109.534118/1‑100 (MQ=60)
                          TTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGA                                                                       >  SRR3722109.431965/1‑100 (MQ=60)
                                  CAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCctgtctcttatacacatctgacgc                                                               >  SRR3722109.369963/1‑76 (MQ=60)
                                                    GAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAG                                             <  SRR3722109.485509/100‑1 (MQ=60)
                                                            AAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGA                                     <  SRR3722109.368311/100‑1 (MQ=60)
                                                              GCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAG                                   <  SRR3722109.305407/100‑1 (MQ=60)
                                                              GCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAG                                   <  SRR3722109.348509/100‑1 (MQ=60)
                                                                    TCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCctgtctcttatacacatctgacgctgccgacgaatgtcaa                             >  SRR3722109.567870/1‑60 (MQ=60)
                                                                    TCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCAC                             >  SRR3722109.366054/1‑100 (MQ=60)
                                                                    TCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCAC                             >  SRR3722109.539591/1‑100 (MQ=60)
                                                                         GCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCG                        >  SRR3722109.385998/1‑100 (MQ=60)
                                                                                   TCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGATGG              <  SRR3722109.572837/100‑1 (MQ=60)
                                                                                             CAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGATGGAGTAGTTAAC    <  SRR3722109.205242/100‑1 (MQ=60)
                                                                                               GCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGATGGAGTAGTTAACGC  <  SRR3722109.126289/100‑1 (MQ=60)
                                                                                                   |                                                                                               
AGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGATGGAGTAGTTAACGC  >  NZ_CP009273/2434430‑2434624

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