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 I181 R1	159	32.7	1279082	96.9%	1239430	121.0

Breseq alignment

BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	788,865	+G	coding (380/1473 nt)	modF ←	molybdate ABC transporter ATP‑binding protein ModF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	788,862	1	.	G	100.0%	101.0 / NA	29	L128P (CTC→CCC)	modF	molybdate ABC transporter ATP‑binding protein ModF
Reads supporting (aligned to +/- strand): ref base . (0/0); new base G (14/15); total (14/15)

CTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGA‑GGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATAT  >  NZ_CP009273/788729‑788979

ctgACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGA‑ggg                                                                                                                    <  1:623333/137‑1 (MQ=255)

              cATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGGTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAAct                                                                                                     >  1:365746/1‑139 (MQ=255)

               aTCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAActg                                                                                                    <  2:453641/139‑1 (MQ=255)

                         tCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGt                                                                                          >  1:25732/1‑139 (MQ=255)

                         tCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGATTTTCGCGTCTCGCCAGTGGAAAGGAATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGt                                                                                          >  1:476202/1‑139 (MQ=255)

                          cGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCCGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGtt                                                                                         >  2:209413/1‑139 (MQ=255)

                          cGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGtt                                                                                         >  1:527531/1‑139 (MQ=255)

                                              aacaaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATcc                                                                     >  1:52300/1‑139 (MQ=255)

                                                    tCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACt                                                               <  1:590173/139‑1 (MQ=255)

                                                        gCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCAt                                                           <  2:527531/139‑1 (MQ=255)

                                                             gACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGa                                                      <  2:9491/139‑1 (MQ=255)

                                                                  cAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGAt                                                 >  1:164289/1‑139 (MQ=255)

                                                                    gcgcCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGAtct                                               <  2:52300/139‑1 (MQ=255)

                                                                    gcgcCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGAtct                                               <  2:164289/139‑1 (MQ=255)

                                                                      gcCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGtt                                                                                         <  1:472810/95‑1 (MQ=255)

                                                                      gcCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGtt                                                                                         >  2:472810/1‑95 (MQ=255)

                                                                               agcagGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAActgc                                                                                                   <  2:389670/76‑1 (MQ=255)

                                                                               agcagGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAActgc                                                                                                   >  1:389670/1‑76 (MQ=255)

                                                                                         ttCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGgtgt                          >  1:404424/1‑139 (MQ=255)

                                                                                             cgTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATc                      <  1:26716/139‑1 (MQ=255)

                                                                                                                  tttAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCatat  <  2:25732/138‑1 (MQ=255)

                                                                                                                  tttAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCatat  <  2:404424/138‑1 (MQ=255)

                                                                                                                    tAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCat    <  1:638329/134‑1 (MQ=255)

                                                                                                                    tAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCat    >  2:638329/1‑134 (MQ=255)

                                                                                                                        gCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTtgcatgca                                                                                  >  1:488560/1‑52 (MQ=255)

                                                                                                                        gCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTtgcatgca                                                                                  <  2:488560/52‑1 (MQ=255)

                                                                                                                          gtcggtcgAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGAt                                                 <  2:126127/83‑1 (MQ=255)

                                                                                                                          gtcggtcgAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGAt                                                 >  1:126127/1‑83 (MQ=255)

                                                                                                                          gtcggtcgAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAg                                      >  2:196926/1‑94 (MQ=255)

                                                                                                                          gtcggtcgAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAg                                      <  1:196926/94‑1 (MQ=255)

CTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGA‑GGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATAT  >  NZ_CP009273/788729‑788979

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output

GCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGA‑GGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCG  >  NZ_CP009273/788728‑789002

GCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACC                                                                                                                                 <  SRR3722135.631475/149‑1 (MQ=60)

    ACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGGTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTG                                                                                                                           >  SRR3722135.370898/1‑151 (MQ=60)

               CATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTT                                                                                                                >  SRR3722135.26071/1‑151 (MQ=60)

               CATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGATTTTCGCGTCTCGCCAGTGGAAAGGAATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTT                                                                                                                >  SRR3722135.482594/1‑151 (MQ=60)

                ATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTG                                                                                                               >  SRR3722135.534525/1‑151 (MQ=60)

                          gtctcgtgggctcggagatgtgtataagagacagCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGC                                                                                                     <  SRR3722135.479169/117‑1 (MQ=60)

                                    CATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCT                                                                                           >  SRR3722135.53004/1‑151 (MQ=60)

                                                    GTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAAT                                                                           <  SRR3722135.597895/151‑1 (MQ=60)

                                                        GGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATC                                                                       >  SRR3722135.166553/1‑151 (MQ=60)

                                                         GCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACctgtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgccg                                                                      >  SRR3722135.444167/1‑96 (MQ=60)

                                                                     GCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCtgtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgc                                                          >  SRR3722135.395103/1‑99 (MQ=60)

                                                                             ggtctcgtgggctcggagatgtgtataagagacagggATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTG                                                  <  SRR3722135.199707/114‑1 (MQ=60)

                                                                               CAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTC                                                >  SRR3722135.410009/1‑151 (MQ=60)

                                                                                             GCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGC                                  <  SRR3722135.27071/151‑1 (MQ=60)

                                                                                                              AGGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCtgtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgccgtcttctgcttgaaaaaaaaaa                  >  SRR3722135.495121/1‑75 (MQ=60)

                                                                                                               GGTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTAT                <  SRR3722135.646615/151‑1 (MQ=60)

                                                                                                                GTATTTAAAGCGTCGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGctgtctcttatacacatctccgagcccacgagaccgtactagatct               >  SRR3722135.127899/1‑105 (MQ=60)

                                                                                                                             GGGTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCG  >  SRR3722135.25635/1‑151 (MQ=60)

                                                                                                                               GTCGAGGAGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCctgtctcttatacacatctccgagcccacgagaccgtact  >  SRR3722135.437521/1‑109 (MQ=60)

GCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGA‑GGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCG  >  NZ_CP009273/788728‑789002

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