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 1,812,807 (C)7→8 intergenic (‑50/+55) celF ← / ← chbR 6‑phospho‑beta‑glucosidase/transcriptional regulator ChbR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,812,8001.C100.0% 81.3 / NA 25intergenic (‑43/+62)celF/chbR6‑phospho‑beta‑glucosidase/transcriptional regulator ChbR
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (15/10);  total (15/10)

TAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAAT  >  NZ_CP009273/1812711‑1812883
                                                                                          |                                                                                    
taGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTg                                                                                       >  1:151054/1‑90 (MQ=255)
taGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTg                                                                                       >  2:532655/1‑90 (MQ=255)
     gctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑ccccc                                                                                >  1:136438/1‑90 (MQ=255)
     gctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑ccccc                                                                                >  2:11439/1‑90 (MQ=255)
     gctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑ccccc                                                                                >  1:230767/1‑90 (MQ=255)
        cccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑cccccccc                                                                             <  2:679901/90‑2 (MQ=255)
        cccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑cccccccc                                                                             <  1:699209/90‑2 (MQ=255)
          cGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑ccccccccgg                                                                           >  1:18201/1‑87 (MQ=255)
             caccaATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACa                                                                         >  2:140421/1‑90 (MQ=255)
                        acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCa                                                              >  2:84359/1‑90 (MQ=255)
                        acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCa                                                              <  2:624206/90‑1 (MQ=255)
                             ttttAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAAc                                                         <  2:64243/90‑1 (MQ=255)
                             ttttAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAAc                                                         <  2:122623/90‑1 (MQ=255)
                             ttttAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAAc                                                         <  2:286740/90‑1 (MQ=255)
                                 aaTTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                     <  2:367935/90‑1 (MQ=255)
                                  aTTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                    <  2:307372/90‑1 (MQ=255)
                                  aTTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                    <  1:521021/90‑1 (MQ=255)
                                              tAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAAt                                        <  2:352974/90‑1 (MQ=255)
                                                    ctcCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCa                                   >  2:626461/1‑90 (MQ=255)
                                                    ctcCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAg                                  >  2:422191/1‑90 (MQ=255)
                                                                  gAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTAc                    >  1:430912/1‑90 (MQ=255)
                                                                  gAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTAc                    >  1:359815/1‑90 (MQ=255)
                                                                    aTACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTg                  >  2:426771/1‑90 (MQ=255)
                                                                    aTACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTg                  >  1:350360/1‑90 (MQ=255)
                                                                       cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGat                                                >  2:156616/1‑57 (MQ=255)
                                                                       cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGat                                                <  1:156616/57‑1 (MQ=255)
                                                                       cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGAtt               >  1:619810/1‑90 (MQ=255)
                                                                            atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt          >  2:400666/1‑90 (MQ=255)
                                                                            atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt          >  2:76726/1‑90 (MQ=255)
                                                                             tCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAttt         >  2:75562/1‑89 (MQ=255)
                                                                                ggCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGt      >  1:358100/1‑90 (MQ=255)
                                                                                   atatCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCaa   >  2:18437/1‑90 (MQ=255)
                                                                                    tatCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAAt  <  2:95848/90‑1 (MQ=255)
                                                                                          |                                                                                    
TAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAAT  >  NZ_CP009273/1812711‑1812883

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ 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
TAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGG  >  NZ_CP009273/1812711‑1812903
                                                                                                 |                                                                                                
TAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCcgg                                                                                               >  SRR3722112.18430/1‑97 (MQ=60)
        CCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAAT                                                                                        <  SRR3722112.710135/100‑1 (MQ=60)
                                  ATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGA                                                              <  SRR3722112.528686/100‑1 (MQ=60)
                                      cgtcagatgtgtataagagacagGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTG                                                          <  SRR3722112.158534/77‑1 (MQ=60)
                                                        CTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTAC                                        >  SRR3722112.364410/1‑100 (MQ=60)
                                                        CTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTAC                                        >  SRR3722112.436761/1‑100 (MQ=60)
                                                          TCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTG                                      >  SRR3722112.354773/1‑100 (MQ=60)
                                                             GTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATT                                   >  SRR3722112.629377/1‑100 (MQ=60)
                                                                  aaaacctaaaATCTGCCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATT                              <  SRR3722112.435278/90‑1 (MQ=60)
                                                                      ACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGT                          >  SRR3722112.362672/1‑100 (MQ=60)
                                                                                             CCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAG   <  SRR3722112.427808/100‑1 (MQ=60)
                                                                                              CCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGG  <  SRR3722112.457503/100‑1 (MQ=60)
                                                                                              CCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGCATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGG  >  SRR3722112.835/1‑100 (MQ=60)
                                                                                                 |                                                                                                
TAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATAGCTCTTAGG  >  NZ_CP009273/1812711‑1812903

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