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 I228 R1 214 26.5 1448192 97.1% 1406194 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,785,687 G→C G75G (GGC→GGG kbl ← glycine C‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,785,6870GC100.0% 108.1 / NA 33G75G (GGC→GGGkblglycine C‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (21/12);  total (21/12)

TTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAA  >  NZ_CP009273/3785603‑3785771
                                                                                    |                                                                                    
ttCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGaa                                                                                 >  1:216890/1‑90 (MQ=255)
   cATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAcc                                                                              <  2:594393/90‑1 (MQ=255)
        ccAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgag                                                                         >  1:425181/1‑90 (MQ=255)
         cAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgaga                                                                        >  2:607046/1‑90 (MQ=255)
         cAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgaga                                                                        >  1:324619/1‑90 (MQ=255)
         cAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgaga                                                                        >  1:485088/1‑90 (MQ=255)
          aGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGaa                                                                       >  1:213467/1‑90 (MQ=255)
           ggAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAAt                                                                      >  2:229237/1‑90 (MQ=255)
                gccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCAtt                                                                 >  1:603447/1‑90 (MQ=255)
                   gccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTccc                                                              <  1:109126/90‑1 (MQ=255)
                              ttCAAGCTCTTTATGGCTGTCCGGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTgccgc                                                   <  1:499568/90‑1 (MQ=255)
                                              cTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGga                                   >  1:237183/1‑90 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  2:664283/1‑90 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  1:498917/1‑90 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  2:6025/1‑90 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  1:142658/1‑90 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  1:69503/1‑90 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCCTTGCCGCCGCAATCAGATCAGgatg                                 >  1:94882/1‑90 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACACACCGAAGCCTTCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  1:531794/1‑90 (MQ=255)
                                                     gAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCg                            >  1:200635/1‑90 (MQ=255)
                                                         gCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAg                        <  2:115992/90‑1 (MQ=255)
                                                          ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            >  1:124623/1‑69 (MQ=255)
                                                          ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            <  2:124623/69‑1 (MQ=255)
                                                            gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                     <  2:40094/90‑1 (MQ=255)
                                                            gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                     <  2:425181/90‑1 (MQ=255)
                                                            gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                     >  1:362874/1‑90 (MQ=255)
                                                            gCAAATAAAACACACCGAAGCCTTCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                     <  2:531794/90‑1 (MQ=255)
                                                               aaTAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            >  2:548075/1‑64 (MQ=255)
                                                               aaTAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            <  1:548075/64‑1 (MQ=255)
                                                                 tAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGAt                <  2:237183/90‑1 (MQ=255)
                                                                          ccGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTgg       >  1:678822/1‑90 (MQ=255)
                                                                               gCCATCCCGAACCCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaa  <  2:494233/90‑1 (MQ=255)
                                                                               gCCATCCCGAAACCGTGAGAACCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaa  <  1:639133/90‑1 (MQ=255)
                                                                                    |                                                                                    
TTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAA  >  NZ_CP009273/3785603‑3785771

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 29 ≤ 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
GTAGAGAATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGA  >  NZ_CP009273/3785588‑3785781
                                                                                                   |                                                                                              
GTAGAGAATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATC                                                                                                >  SRR3722112.123416/1‑100 (MQ=60)
GTAGAGAATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATC                                                                                                >  SRR3722112.673756/1‑100 (MQ=60)
     GAATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAA                                                                                           >  SRR3722112.219501/1‑100 (MQ=60)
       ATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAC                                                                                         <  SRR3722112.723536/100‑1 (MQ=60)
             TCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAG                                                                                   >  SRR3722112.430867/1‑100 (MQ=60)
              CTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGA                                                                                  >  SRR3722112.492061/1‑100 (MQ=60)
              CTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGA                                                                                  >  SRR3722112.328460/1‑100 (MQ=60)
               TTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAA                                                                                 >  SRR3722112.216043/1‑100 (MQ=60)
                     CCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATT                                                                           >  SRR3722112.612691/1‑100 (MQ=60)
                                  GCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCG                                                              <  SRR3722112.110482/100‑1 (MQ=60)
                                             TTCAAGCTCTTTATGGCTGTCCGGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGA                                                   <  SRR3722112.506807/100‑1 (MQ=60)
                                                   CTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGA                                             >  SRR3722112.240012/1‑100 (MQ=60)
                                                    tgtgtataagagacagGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGAT                                            <  SRR3722112.556337/84‑1 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                           >  SRR3722112.144442/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                           >  SRR3722112.70348/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                           >  SRR3722112.506148/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCCTTGCCGCCGCAATCAGATCAGGATG                                           >  SRR3722112.96065/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACACACCGAAGCCTTCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                           >  SRR3722112.539682/1‑100 (MQ=60)
                                                          TGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCG                                      >  SRR3722112.203071/1‑100 (MQ=60)
                                                               GTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATctgtctcttat                                 >  SRR3722112.126163/1‑89 (MQ=60)
                                                                 CCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCC                               >  SRR3722112.367524/1‑100 (MQ=60)
                                                                               ATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGG                 >  SRR3722112.689434/1‑100 (MQ=60)
                                                                                       CACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAG         >  SRR3722112.126274/1‑100 (MQ=60)
                                                                                       CACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAG         >  SRR3722112.119438/1‑100 (MQ=60)
                                                                                             AGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATG   >  SRR3722112.543481/1‑100 (MQ=60)
                                                                                              GCCATCCCGAAACCGTGAGAACCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGA  <  SRR3722112.649062/100‑1 (MQ=60)
                                                                                                   |                                                                                              
GTAGAGAATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGA  >  NZ_CP009273/3785588‑3785781

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ 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: