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 F3 I122 R1 193 20.9 1228102 94.7% 1163012 84.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,210,075 T→C F518L (TTT→CTT)  aceK → bifunctional isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,210,0750TC100.0% 104.5 / NA 29F518L (TTT→CTT) aceKbifunctional isocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/19);  total (10/19)

GACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGTTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCG  >  NZ_CP009273/4209988‑4210160
                                                                                       |                                                                                     
gACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGctt                                                                                     >  2:38890/1‑90 (MQ=255)
cccctCCCGCCGCCGCGCCCCCCGGAAGTTTACCTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGctt                                                                                     <  1:375015/86‑1 (MQ=255)
       cgccgccgcGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACt                                                                              <  1:133014/90‑1 (MQ=255)
       cgccgccgcGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACt                                                                              <  2:311141/90‑1 (MQ=255)
                  tATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGcc                                                                   <  1:156635/90‑1 (MQ=255)
                      cGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGAcc                                                               <  2:268200/90‑1 (MQ=255)
                      cGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGAcc                                                               <  1:513570/90‑1 (MQ=255)
                        gAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCcg                                                             <  2:89864/90‑1 (MQ=255)
                           gACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGt                                                          >  2:602039/1‑90 (MQ=255)
                             cGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACt                                                                              <  2:176235/68‑1 (MQ=255)
                             cGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACt                                                                              >  1:176235/1‑68 (MQ=255)
                             cGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGAcc                                                               <  1:519735/83‑1 (MQ=255)
                             cGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGAcc                                                               >  2:519735/1‑83 (MQ=255)
                             cGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTAt                                                        >  2:518112/1‑90 (MQ=255)
                                     ccAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGc                                                <  1:452640/90‑1 (MQ=255)
                                        gCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGt                                             >  2:399619/1‑90 (MQ=255)
                                                   tACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATg                                  <  1:62846/90‑1 (MQ=255)
                                                   tACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATg                                  <  2:421046/90‑1 (MQ=255)
                                                       gCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGc                                                <  2:267555/72‑1 (MQ=255)
                                                       gCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGc                                                >  1:267555/1‑72 (MQ=255)
                                                        cGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGc                             <  2:58811/90‑1 (MQ=255)
                                                        cGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGc                             <  2:182325/90‑1 (MQ=255)
                                                         gTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGcc                            >  1:116027/1‑90 (MQ=255)
                                                         gTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGcc                            >  2:603171/1‑90 (MQ=255)
                                                           ctcGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGAcc                                                               >  1:14414/1‑53 (MQ=255)
                                                           ctcGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGAcc                                                               <  2:14414/53‑1 (MQ=255)
                                                            tcGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGc                             <  2:289894/86‑1 (MQ=255)
                                                            tcGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGc                             >  1:289894/1‑86 (MQ=255)
                                                             cGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCATATTGGTCCGCTGTTTGAAGAGATGCACGCCGAcc                        <  1:496315/90‑1 (MQ=255)
                                                                       tGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCgcgc              <  1:602039/90‑1 (MQ=255)
                                                                                   agagCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGgcg  <  1:586724/90‑1 (MQ=255)
                                                                                       |                                                                                     
GACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGTTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCG  >  NZ_CP009273/4209988‑4210160

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ 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
TTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGTTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAAC  >  NZ_CP009273/4209982‑4210173
                                                                                             |                                                                                                  
agatctacactccttacgtcgtcggcagcgtcagatgtgtataagagacagCCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCC                                                                                              <  SRR3721989.398090/49‑1 (MQ=60)
  CCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCAC                                                                                            <  SRR3721989.425765/100‑1 (MQ=60)
      cccctcccgccgccgcgccccccggaagtttacCTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGC                                                                                        <  SRR3721989.382436/67‑1 (MQ=60)
             CGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGC                                                                                 <  SRR3721989.135445/100‑1 (MQ=60)
                        TATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTA                                                                      <  SRR3721989.159477/100‑1 (MQ=60)
                         ATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCtgtctcttata                                                                     >  SRR3721989.179459/1‑89 (MQ=60)
                            CGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGG                                                                  <  SRR3721989.524533/100‑1 (MQ=60)
                            CGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGG                                                                  <  SRR3721989.530841/100‑1 (MQ=60)
                                           CCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGA                                                   <  SRR3721989.462104/100‑1 (MQ=60)
                                                   CCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGActgtctct                                           >  SRR3721989.272307/1‑92 (MQ=60)
                                                     GTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCC                                         >  SRR3721989.118150/1‑100 (MQ=60)
                                                       GGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGctgtctcttatacacatctgacgctgc                                       >  SRR3721989.14678/1‑73 (MQ=60)
                                                        GTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCCGAC                                      >  SRR3721989.295249/1‑100 (MQ=60)
                                                         TACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACC                                     <  SRR3721989.63954/100‑1 (MQ=60)
                                                                   CGCCGGGCGATGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCATATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGC                           <  SRR3721989.506834/100‑1 (MQ=60)
                                                                             TGTTTTCCCGGAAGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGG                 <  SRR3721989.615120/100‑1 (MQ=60)
                                                                                         AGAGCTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAA     <  SRR3721989.599485/100‑1 (MQ=60)
                                                                                            GCTTCGCCACTGGCTATGCGCCGCCCCGCGTATTGGTCCGCTGTTTGAAGAGATGCCCGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAAC  <  SRR3721989.575532/100‑1 (MQ=60)
                                                                                             |                                                                                                  
TTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGTTTCGCCACTGGCTATGCGCCGACCCGCGTATTGGTCCGCTGTTTGAAGAGATGCACGCCGACCTGTTCCGCGCTGATTACTGGCGCGCACTACAAAAC  >  NZ_CP009273/4209982‑4210173

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