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
A2 F5 I179 R1 182 40.3 1547118 96.8% 1497610 120.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 634,370 A→C intergenic (‑341/‑31) dsbG ← / → ahpC thiol:disulfide interchange protein DsbG/alkyl hydroperoxide reductase subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273634,3700AC100.0% 104.6 / NA 27intergenic (‑341/‑31)dsbG/ahpCthiol:disulfide interchange protein DsbG/alkyl hydroperoxide reductase subunit C
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (12/15);  total (12/15)

AGGGGGCCATTTTACTTTTTATCGCCGCTGGCGGTGCAAAGTTCACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTACATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGAATTCATCGAAATCACCGAAAAAGATACCGAAGGCCGCTGGAGCG  >  NZ_CP009273/634237‑634503
                                                                                                                                     |                                                                                                                                     
aGGGGGCCATTTTACTTTTTATCGCCGCTGGCGGTGCAAAGTTCACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCAtct                                                                                                                                  >  2:237654/1‑139 (MQ=255)
  ggggCCATTTTACTTTTTATCGCCGCTGGCGGTGCAAAGTTCACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCa                                                                                                                                >  2:572228/1‑139 (MQ=255)
    ggCCATTTTACTTTTTATCGCCGCTGGCGGTGCAAAGTTCACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATc                                                                                                                              >  2:478413/1‑139 (MQ=255)
              cTTTTTATCGCCGCTGGCGGTGCAAAGTTCACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACgg                                                                                                                    <  2:704609/139‑1 (MQ=255)
                 tttATCGCCGCTGGCGGTGCAAAGTTCACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACgg                                                                                                                    <  2:468323/136‑1 (MQ=255)
                  ttATCGCCGCTGGCGGTGCAAAGTTCACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACggagg                                                                                                                 <  1:433143/138‑1 (MQ=255)
                  ttATCGCCGCTGGCGGTGCAAAGTTCACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACggagg                                                                                                                 <  1:692991/138‑1 (MQ=255)
                        ccgcTGGCGGTGCAAAGTTCACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAGta                                                                                                            >  1:369106/1‑137 (MQ=255)
                          gcTGGCGGTGCAAAGTTCACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGa                                                                                                        <  2:769554/139‑1 (MQ=255)
                                           cacaAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAacac                                                                                        >  2:744010/1‑138 (MQ=255)
                                                               gTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGaaaa                                                                                                                         <  1:123466/85‑1 (MQ=255)
                                                               gTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGaaaa                                                                                                                         >  2:123466/1‑85 (MQ=255)
                                                                  gTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAAcc                                                                 <  2:92969/138‑1 (MQ=255)
                                                                                 tCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTa                                                                                <  1:742345/108‑1 (MQ=255)
                                                                                 tCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTa                                                                                >  2:742345/1‑108 (MQ=255)
                                                                                      ttGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACggagg                                                                                                                 >  1:648635/1‑70 (MQ=255)
                                                                                      ttGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACggagg                                                                                                                 <  2:648635/70‑1 (MQ=255)
                                                                                        gATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTTCCTTCCATCTCATCGAAAACACGGAGGaa                                                                                                               <  2:438223/70‑1 (MQ=255)
                                                                                        gATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGaa                                                                                                               >  1:438223/1‑70 (MQ=255)
                                                                                             tGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGAATTCATCg                                      >  1:560672/1‑138 (MQ=255)
                                                                                                aaaCGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGAATTCATCGAAAt                                  <  2:560672/139‑1 (MQ=255)
                                                                                                       ttAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGAATTCATCGAAATCACCGaa                           <  2:224385/139‑1 (MQ=255)
                                                                                                        tAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGAATTCATCGAAATCACCGaa                           >  1:312777/1‑138 (MQ=255)
                                                                                                                   gCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGAATTCATCGAAATCACCGAAAAAGATACCGaa               <  2:6070/139‑1 (MQ=255)
                                                                                                                    cAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGAATTCATCGAAATCACCGAAAAAGATACCGAAg              >  2:410005/1‑139 (MQ=255)
                                                                                                                             gTTACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGAATTCATCGAAATCACCGAAAAAGATACCGAAGGCCGCTGGa     <  1:572228/139‑1 (MQ=255)
                                                                                                                                 ccTTCCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGAATTCATCGAAATCACCGAAAAAGATACCGAAGGCCGCTGGAGCg  <  1:478413/138‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                     
AGGGGGCCATTTTACTTTTTATCGCCGCTGGCGGTGCAAAGTTCACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTACATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGAATTCATCGAAATCACCGAAAAAGATACCGAAGGCCGCTGGAGCG  >  NZ_CP009273/634237‑634503

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGGTCACCTGGAAAGGGGGCCATTTTACTTTTTATCGCCGCTGGCGGTGCAAAGTTCACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTACATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGAATTCATCGAAATCACCGAAAAAGATACCGAAGGCCGCTGGAGCGTCTTCTTCTTC  >  NZ_CP009273/634224‑634514
                                                                                                                                                  |                                                                                                                                                
agacagCCTGGAAAGGGGGCCATTTTACTTTTTATCGCCGCTGGCGGTGCAAAGTTCACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATC                                                                                                                                              <  SRR3722133.223367/145‑1 (MQ=60)
                         agtacggtctcgtgggctcggagatgtgtataagagacagGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAG                                                                                                                         <  SRR3722133.125096/107‑1 (MQ=60)
                          ACTTTTTATCGCCGCTGGCGGTGCAAAGTTCACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAGTAT                                                                                                                      >  SRR3722133.374277/1‑149 (MQ=60)
                              TTTATCGCCGCTGGCGGTGCAAAGTTCACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGATG                                                                                                                 <  SRR3722133.439065/150‑1 (MQ=60)
                              TTTATCGCCGCTGGCGGTGCAAAGTTCACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGATG                                                                                                                 <  SRR3722133.701985/150‑1 (MQ=60)
                                                              ggagatgtgtataagagacagGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAA                                                                                <  SRR3722133.751941/130‑1 (MQ=60)
                                                                                        AACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGATGctgtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgccgtctt                                                      >  SRR3722133.657101/1‑92 (MQ=60)
                                                                                          CTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTgtctcttatacacatctccgagcccacgagaccgtactagatatcgtatgccgaatt                                                    >  SRR3722133.444196/1‑94 (MQ=60)
                                                                                               CGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGAATTCATCGA                                                >  SRR3722133.568128/1‑150 (MQ=60)
                                                                                                          TGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGAATTCATCGAAATCACCGAAA                                     >  SRR3722133.317076/1‑150 (MQ=60)
                                                                                                                                         GGTTACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGAATTCATCGAAATCACCGAAAAAGATACCGAAGGCCGCTGGAGCGTCTTCTTC     <  SRR3722133.579803/151‑1 (MQ=60)
                                                                                                                                             ACCTTCCATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGAATTCATCGAAATCACCGAAAAAGATACCGAAGGCCGCTGGAGCGTCTTCTTCTTC  <  SRR3722133.484890/150‑1 (MQ=60)
                                                                                                                                                  |                                                                                                                                                
CGGTCACCTGGAAAGGGGGCCATTTTACTTTTTATCGCCGCTGGCGGTGCAAAGTTCACAAAGTTGTCTTACGAAGGTTGTAAGGTAAAACTTATCGATTTGATAATGGAAACGCATTAGCCGAATCGGCAAAAATTGGTTACCTTACATCTCATCGAAAACACGGAGGAAGTATAGATGTCCTTGATTAACACCAAAATTAAACCTTTTAAAAACCAGGCATTCAAAAACGGCGAATTCATCGAAATCACCGAAAAAGATACCGAAGGCCGCTGGAGCGTCTTCTTCTTC  >  NZ_CP009273/634224‑634514

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