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 945,384 (T)5→4 coding (261/591 nt) ycaK → NAD(P)H oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273945,3800T.100.0% 88.8 / NA 21coding (257/591 nt)ycaKNAD(P)H oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (6/15);  total (6/15)

CGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGGCGCTGGTTG  >  NZ_CP009273/945256‑945511
                                                                                                                            |                                                                                                                                   
cGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCtggtg                                                                                                                      >  2:234164/1‑139 (MQ=255)
    aTTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTAc                                                                                                                  <  1:335846/139‑1 (MQ=255)
    aTTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTAc                                                                                                                  <  2:609371/139‑1 (MQ=255)
                          ggAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAg                                                                                            <  1:324242/139‑1 (MQ=255)
                                    ccAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGa                                                                                  <  1:234164/139‑1 (MQ=255)
                                      aGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACa                                                                                <  1:248601/139‑1 (MQ=255)
                                            gaagaaTATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTa                                                                           <  1:492848/138‑1 (MQ=255)
                                            gaagaaTATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTa                                                                           <  2:303439/138‑1 (MQ=255)
                                                       aTAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATggg                                                               <  2:645272/139‑1 (MQ=255)
                                                            cGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGggctggc                                                          <  1:340233/139‑1 (MQ=255)
                                                                         agGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATggg                                             <  1:605250/139‑1 (MQ=255)
                                                                              catcaGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATggg                                                               >  1:659516/1‑116 (MQ=255)
                                                                              catcaGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATggg                                                               <  2:659516/116‑1 (MQ=255)
                                                                                                   ctTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACa                                                                                <  1:63810/78‑1 (MQ=255)
                                                                                                   ctTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACa                                                                                >  2:63810/1‑78 (MQ=255)
                                                                                                             acacGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGgcg         <  2:31175/139‑1 (MQ=255)
                                                                                                                 gTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGgcg         <  2:416049/135‑1 (MQ=255)
                                                                                                                 gTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGgcg         >  1:416049/1‑135 (MQ=255)
                                                                                                                      tggtggTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATggg                                                               >  2:308046/1‑76 (MQ=255)
                                                                                                                      tggtggTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATggg                                                               <  1:308046/76‑1 (MQ=255)
                                                                                                                      tggtggTTTTCCTCTCTGGTGGGACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGGCGCTGGTTg  >  2:95154/1‑137 (MQ=255)
                                                                                                                            |                                                                                                                                   
CGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGGCGCTGGTTG  >  NZ_CP009273/945256‑945511

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ 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
AATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGGCGCTGGTTGGAGGAGA  >  NZ_CP009273/945259‑945518
                                                                                                                             |                                                                                                                                      
AATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGC                                                                                                              <  SRR3722133.340576/151‑1 (MQ=60)
                      CGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGAC                                                                                        <  SRR3722133.328725/151‑1 (MQ=60)
                                ACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGA                                                                              <  SRR3722133.237354/151‑1 (MQ=60)
                                  atacgagatctagtacggtctcgtgggctcggagatgtgtataagagacagATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAAT                                                                            <  SRR3722133.64628/100‑1 (MQ=60)
                                  CAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAAT                                                                            <  SRR3722133.252001/151‑1 (MQ=60)
                                        GGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGG                                                                       <  SRR3722133.499498/150‑1 (MQ=60)
                                                   gcatacgagatctagtacggtctcgtgggctcggagatgtgtataagagacagTGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGG                                                           <  SRR3722133.312279/98‑1 (MQ=60)
                                                        ACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATG                                                      <  SRR3722133.345011/151‑1 (MQ=60)
                                                             ATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACctgtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgccgtcttctgcttgaaaaa                                                  >  SRR3722133.98809/1‑79 (MQ=60)
                                                                CTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGctgtctcttatac                                              >  SRR3722133.668111/1‑138 (MQ=60)
                                                                     GAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACC                                         <  SRR3722133.613251/151‑1 (MQ=60)
                                                                                                   GAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGGCGCTGGT           >  SRR3722133.421797/1‑151 (MQ=60)
                                                                                                             CGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATctgtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgccgtcttctgcttgaaaaaaaa  >  SRR3722133.510119/1‑76 (MQ=60)
                                                                                                             CGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATGGGCACAAATTctgtctcttatacacatctccgagcccacgagaccgtact    >  SRR3722133.187371/1‑108 (MQ=60)
                                                                                                                             |                                                                                                                                      
AATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTCAATAAAGTTCGTTGGGTGGCGCTGGTTGGAGGAGA  >  NZ_CP009273/945259‑945518

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: