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 I96 R1 172 34.7 2005794 96.2% 1929573 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,662,649 T→G H58P (CAC→CCC)  mlc ← DNA‑binding transcriptional repressor Mlc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,662,6490TG100.0% 104.3 / NA 31H58P (CAC→CCC) mlcDNA‑binding transcriptional repressor Mlc
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (14/17);  total (14/17)

GTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGTGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATACGCGAGACTGGA  >  NZ_CP009273/1662564‑1662735
                                                                                     |                                                                                      
gTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGc                                                                                    >  1:24975/1‑90 (MQ=255)
gTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGc                                                                                    >  1:548153/1‑90 (MQ=255)
  gCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCtt                                                                                  >  2:666216/1‑90 (MQ=255)
  gCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCtt                                                                                  >  1:693530/1‑90 (MQ=255)
  gCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCtt                                                                                  <  2:465759/90‑1 (MQ=255)
     aGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGa                                                                               >  2:78299/1‑90 (MQ=255)
      ggCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAg                                                                              >  2:319541/1‑90 (MQ=255)
        cTTCAGTTTCAACCCCCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCa                                                                            <  1:989840/90‑1 (MQ=255)
        cTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCa                                                                            <  2:455404/90‑1 (MQ=255)
          tCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCAtc                                                                          >  2:320027/1‑90 (MQ=255)
                  aaCCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGAc                                                                  <  2:71914/90‑1 (MQ=255)
                       ccaGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAAttt                                                             >  1:961770/1‑90 (MQ=255)
                          gCCCCCCCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAg                                                          <  2:17725/90‑1 (MQ=255)
                             ccACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGa                                                       <  1:200484/90‑1 (MQ=255)
                              cACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGa                                                                   <  1:701833/77‑1 (MQ=255)
                              cACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGa                                                                   >  2:701833/1‑77 (MQ=255)
                              cACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGAt                                                      >  2:44596/1‑90 (MQ=255)
                                                gTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTg                                    >  2:991542/1‑90 (MQ=255)
                                                    cccGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGcc                                >  2:805921/1‑90 (MQ=255)
                                                        cttcttTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGAc                            <  1:862752/90‑1 (MQ=255)
                                                        cttcttTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGAc                            <  2:803604/90‑1 (MQ=255)
                                                                  ttCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATc                  <  2:974650/90‑1 (MQ=255)
                                                                    ccAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGa                >  2:584487/1‑90 (MQ=255)
                                                                    ccAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGa                <  1:584487/90‑1 (MQ=255)
                                                                         tctTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCATGAGCCAGTTGCGCCAGACGGGAAAGATCGATAcgc           <  2:467151/90‑1 (MQ=255)
                                                                         tctTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATAcgc           <  2:548153/90‑1 (MQ=255)
                                                                         tctTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATAcgc           <  1:304831/90‑1 (MQ=255)
                                                                         tctTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATAcgc           <  2:355663/90‑1 (MQ=255)
                                                                         tctTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATAcgc           <  2:305798/90‑1 (MQ=255)
                                                                          ctTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATAcgcg          <  1:991542/90‑1 (MQ=255)
                                                                                  aGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGACCGTTGGGCCAGACGGGCAAGATCGCAACGCGCGACTGGa  >  1:994568/1‑90 (MQ=255)
                                                                                     |                                                                                      
GTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGTGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATACGCGAGACTGGA  >  NZ_CP009273/1662564‑1662735

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ 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
GCAGAGAAAGATAGTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGTGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATACGCGAGACTGGACCAAGCTGATCAA  >  NZ_CP009273/1662551‑1662748
                                                                                                  |                                                                                                   
gtataagagacagGTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGG                                                                                                    <  SRR3722236.768932/87‑1 (MQ=60)
cagGAGAAAGATAGTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGG                                                                                                    <  SRR3722236.673654/97‑1 (MQ=60)
GCAGAGAAAGATAGTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGG                                                                                                    <  SRR3722236.972624/100‑1 (MQ=60)
   GAGAAAGATAGTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGC                                                                                                 >  SRR3722236.25335/1‑100 (MQ=60)
   GAGAAAGATAGTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGC                                                                                                 >  SRR3722236.556603/1‑100 (MQ=60)
   GAGAAAGATAGTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGC                                                                                                 >  SRR3722236.492189/1‑100 (MQ=60)
     GAAAGATAGTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTT                                                                                               >  SRR3722236.705214/1‑100 (MQ=60)
     GAAAGATAGTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTT                                                                                               <  SRR3722236.404579/100‑1 (MQ=60)
                     CTTCAGTTTCAACCCCCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGAC                                                                               <  SRR3722236.1007279/100‑1 (MQ=60)
                          GTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTT                                                                          >  SRR3722236.978740/1‑100 (MQ=60)
                              CAGCCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGT                                                                      <  SRR3722236.713699/100‑1 (MQ=60)
                                          CCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGG                                                          <  SRR3722236.203472/100‑1 (MQ=60)
                                                                     CTTCTTTGATTTCCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATC                               <  SRR3722236.877883/100‑1 (MQ=60)
                                                                                 CCAGCTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATACGCGAGAC                   <  SRR3722236.593591/100‑1 (MQ=60)
                                                                                     CTCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGACCGTTGGGCCAGACGGGCAAGATCGCAACGCGCGACTGGA               >  SRR3722236.1012080/1‑100 (MQ=60)
                                                                                      TCTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATACGCGAGACTGGAC              <  SRR3722236.309226/100‑1 (MQ=60)
                                                                                       CTTGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATACGCGAGACTGGACC             <  SRR3722236.1009005/100‑1 (MQ=60)
                                                                                         TGCACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATACGCGAGACTGGACCAA           >  SRR3722236.26586/1‑100 (MQ=60)
                                                                                           CACCAGGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATACGCGAGACTGGACCAAGC         >  SRR3722236.944624/1‑100 (MQ=60)
                                                                                                GGGGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATACGCGAGACTGGACCAAGCTGATC    >  SRR3722236.914262/1‑100 (MQ=60)
                                                                                                  GGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATACGCGAGACTGGACCAAGCTGATCAA  <  SRR3722236.420397/100‑1 (MQ=60)
                                                                                                  |                                                                                                   
GCAGAGAAAGATAGTGCCAGGCTTCAGTTTCAACCACCAGCCCCACCGCCGGACGGCCACGGTTCCCCGCTTCTTTGATTTCCAGCTCTTGCACCAGGTGTGCTTCGAGCATCTCACGGACAATTTTAGTGATACTGGCAGGAGCCAGTTGCGCCAGACGGGAAAGATCGATACGCGAGACTGGACCAAGCTGATCAA  >  NZ_CP009273/1662551‑1662748

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