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 F1 I1 R1 744 338.8 17650688 97.8% 17262372 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,041,671 (T)6→7 intergenic (+46/‑108) dsbA → / → b3861 protein disulfide isomerase I/putative GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,041,6651.T100.0% 45.4 / NA 14intergenic (+40/‑114)dsbA/b3861protein disulfide isomerase I/putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (8/6);  total (8/6)

TTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAGTATTCACGGGCTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGATATTCAAAGTT  >  NC_000913/4041573‑4041752
                                                                                             |                                                                                       
tttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGc                                                                                          <  1:7599360/93‑1 (MQ=255)
tttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGc                                                                                          <  2:8472223/93‑1 (MQ=255)
tttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGc                                                                                          <  2:4374981/93‑1 (MQ=255)
tttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGc                                                                                          <  2:4374980/93‑1 (MQ=255)
 ttttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCt                                                                                        <  2:7868267/93‑1 (MQ=255)
  tttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCtt                                                                                       <  1:4297784/93‑1 (MQ=255)
  tttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCtt                                                                                       <  2:6497413/93‑1 (MQ=255)
  tttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCtt                                                                                       <  2:6108315/93‑1 (MQ=255)
  tttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCtt                                                                                       <  1:2293569/93‑1 (MQ=255)
  tttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCtt                                                                                       <  2:3656956/93‑1 (MQ=255)
  tttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCtt                                                                                       <  2:2515707/93‑1 (MQ=255)
   ttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCttt                                                                                      <  1:1112310/93‑1 (MQ=255)
   ttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCttt                                                                                      <  1:3479754/93‑1 (MQ=255)
   ttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCttt                                                                                      <  1:4189169/93‑1 (MQ=255)
   ttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCttt                                                                                      <  2:6167099/93‑1 (MQ=255)
   ttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCttt                                                                                      <  1:5352179/93‑1 (MQ=255)
   ttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCttt                                                                                      <  1:1140353/93‑1 (MQ=255)
   ttGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCttt                                                                                      <  1:6868474/93‑1 (MQ=255)
   ttGTTCAGCAGTAGGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAAGTCAATGTAAATTACATAAAGCCCGTGAATATTCCCGGGCttt                                                                                      <  1:3479753/93‑1 (MQ=255)
    tGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCtttt                                                                                     <  2:1391884/93‑1 (MQ=255)
     gTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCttttt                                                                                    <  1:4880886/93‑1 (MQ=255)
       tCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCttttttt                                                                                  <  2:1119032/93‑2 (MQ=255)
        cagcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCttttttta                                                                                 <  1:8138841/93‑3 (MQ=255)
        cagcagTATGATGATACAGTGAAATATCTGTCTGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCCCGGGCttttttta                                                                                 <  2:6367527/93‑3 (MQ=255)
          gcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCtttttttatt                                                                               <  1:6661332/93‑5 (MQ=255)
          gcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCtttttttatt                                                                               <  1:5649487/93‑5 (MQ=255)
          gcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCtttttttatt                                                                               <  2:2784316/93‑5 (MQ=255)
          gcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCtttttttatt                                                                               <  2:3404192/93‑5 (MQ=255)
          gcagTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCtttttttatt                                                                               <  2:364916/93‑5 (MQ=255)
             gTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTattatt                                                                             <  1:4726437/93‑1 (MQ=255)
             gTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTattatt                                                                             <  2:1023394/93‑1 (MQ=255)
              tATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTAttt                                                                            >  1:3442299/1‑93 (MQ=255)
               aTGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTta                                                                           <  1:6202015/93‑1 (MQ=255)
                  cTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTtaata                                                                        <  1:1056587/93‑1 (MQ=255)
                                        cGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAAt                                                  >  1:5810196/1‑93 (MQ=255)
                                            aaaaaaTAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCAt                                              >  1:5628394/1‑93 (MQ=255)
                                              aaaaTAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATcc                                            <  1:50484/93‑1 (MQ=255)
                                                    atTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGc                                      >  2:4989555/1‑93 (MQ=255)
                                                     ttCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCt                                     >  1:7386408/1‑93 (MQ=255)
                                                              aaTTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACattat                            <  1:4061095/93‑1 (MQ=255)
                                                                              tGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGat            >  2:4216208/1‑93 (MQ=255)
                                                                                  tATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGATATTc        >  2:2432782/1‑93 (MQ=255)
                                                                                        cGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGATATTCAAAGtt  >  1:8769106/1‑93 (MQ=255)
                                                                                             |                                                                                       
TTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAGTATTCACGGGCTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGATATTCAAAGTT  >  NC_000913/4041573‑4041752

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ 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
GTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAGTATTCACGGGCTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGATATTCAAAGTT  >  NC_000913/4041578‑4041752
                                                                                              |                                                                                 
GTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTT                                                                                   <  SRR4302541.4930699/94‑1 (MQ=60)
   CAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTtat                                                                                <  SRR4302541.8227898/94‑4 (MQ=60)
     GCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTA                                                                               <  SRR4302541.6729678/94‑1 (MQ=60)
     GCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTA                                                                               <  SRR4302541.5707208/94‑1 (MQ=60)
      CAGTATGCTGATACAGTGAAATATCTGTGCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTAT                                                                              <  SRR4302541.1029515/94‑1 (MQ=60)
      CAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTAT                                                                              <  SRR4302541.4197051/94‑1 (MQ=60)
      CAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTAT                                                                              <  SRR4302541.6148162/94‑1 (MQ=60)
        GTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTT                                                                            >  SRR4302541.3476697/1‑94 (MQ=60)
        GTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTT                                                                            <  SRR4302541.4774642/94‑1 (MQ=60)
          ATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAA                                                                          <  SRR4302541.6265391/94‑1 (MQ=60)
             CTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAA                                                                       <  SRR4302541.1066404/94‑1 (MQ=60)
                                  CCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAAT                                                  >  SRR4302541.5869227/1‑94 (MQ=60)
                                      GAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCAT                                              >  SRR4302541.5685972/1‑94 (MQ=60)
                                         AAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCT                                           <  SRR4302541.50872/94‑1 (MQ=60)
                                               ATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCT                                     >  SRR4302541.7465389/1‑94 (MQ=60)
                                                         AATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATT                           <  SRR4302541.4101167/94‑1 (MQ=60)
                                                                                  ACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGATATTCAAAGTT  >  SRR4302541.8862858/1‑94 (MQ=60)
                                                                                              |                                                                                 
GTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAGTATTCACGGGCTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCTGGACATTATTTCAGCATAAGTGGATATTCAAAGTT  >  NC_000913/4041578‑4041752

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