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 F4 I96 R1	227	23.3	882352	98.0%	864704	126.9

Breseq alignment

BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	453,912	A→G	E344G (GAA→GGA)	clpX →	ATP‑dependent protease ATP‑binding subunit ClpX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	453,912	0	A	G	96.9%	109.7 / ‑3.2	32	E344G (GAA→GGA)	clpX	ATP‑dependent protease ATP‑binding subunit ClpX
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (15/16); total (15/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

TCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGAAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATAC  >  NZ_CP009273/453799‑454029

tcgtcTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGAc                                                                                              >  1:4357/1‑139 (MQ=255)

  gtcTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACg                                                                                             <  2:401914/138‑1 (MQ=255)

          gTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTg                                                                                     >  2:105510/1‑138 (MQ=255)

             gTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATg                                                                                 >  2:333126/1‑139 (MQ=255)

                 cAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTAt                                                                             >  2:37695/1‑139 (MQ=255)

                 cAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTAt                                                                             >  1:246555/1‑139 (MQ=255)

                 cAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTAt                                                                             <  1:41552/139‑1 (MQ=255)

                        gaatgaaCTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTaa                                                                       <  1:105510/138‑1 (MQ=255)

                                        gaagCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTa                                                      <  1:37695/139‑1 (MQ=255)

                                            ctctGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAAc                                                  >  2:51463/1‑139 (MQ=255)

                                              ctGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCg                                                >  2:283404/1‑139 (MQ=255)

                                                     aGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGccc                                          <  2:257982/138‑1 (MQ=255)

                                                      gATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGt                                        >  2:238106/1‑139 (MQ=255)

                                                           ccAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCt                                   >  1:322253/2‑139 (MQ=255)

                                                                gagCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGtt                              <  2:424488/139‑1 (MQ=255)

                                                                 agCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTc                             >  1:87075/1‑139 (MQ=255)

                                                                   ccGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGgcg                                                                                        <  1:378301/78‑1 (MQ=255)

                                                                   ccGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGgcg                                                                                        >  2:378301/1‑78 (MQ=255)

                                                                               cTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCg               <  1:238106/139‑1 (MQ=255)

                                                                                  aCCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCa             <  2:246555/138‑1 (MQ=255)

                                                                                  aCCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCa             <  1:51463/138‑1 (MQ=255)

                                                                                  aCCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCa             <  1:333126/138‑1 (MQ=255)

                                                                                  aCCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCa             <  1:398248/138‑1 (MQ=255)

                                                                                  aCCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCa             <  2:382140/138‑1 (MQ=255)

                                                                                       gCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTc       >  2:52517/1‑139 (MQ=255)

                                                                                           tATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATa   <  2:87075/139‑1 (MQ=255)

                                                                                           tATCAGGCGCTGTTTAATCTGGAAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATa   <  1:314926/139‑1 (MQ=255)

                                                                                            aTCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGGGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATAc  >  1:115228/1‑139 (MQ=255)

                                                                                                     tGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCa             >  2:181540/1‑119 (MQ=255)

                                                                                                     tGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCa             <  1:181540/119‑1 (MQ=255)

                                                                                                     tGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATa   >  1:434576/1‑129 (MQ=255)

                                                                                                     tGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATa   <  2:434576/129‑1 (MQ=255)

TCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGAAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATAC  >  NZ_CP009273/453799‑454029

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

GAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGAAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATACCATGTACGATCTGCCGTCCATGGAAG  >  NZ_CP009273/453788‑454055

GAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACG                                                                                                                       >  SRR3722119.4397/1‑151 (MQ=60)

                atacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagTCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTAT                                                                                                       <  SRR3722119.382110/100‑1 (MQ=60)

                 GCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATC                                                                                                      >  SRR3722119.249116/1‑151 (MQ=60)

                           GCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAG                                                                                            <  SRR3722119.41923/151‑1 (MQ=60)

                                  TGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGG                                                                                      <  SRR3722119.106547/150‑1 (MQ=60)

                                                  AGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCC                                                                     <  SRR3722119.38027/151‑1 (MQ=60)

                                                           GATTCAGATCCCCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTG                                                            >  SRR3722119.325544/1‑151 (MQ=60)

                                                                 GATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCC                                                      >  SRR3722119.87930/1‑151 (MQ=60)

                                                                                         CCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGAT                              <  SRR3722119.240593/151‑1 (MQ=60)

                                                                                            aagagacagGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATAC                            <  SRR3722119.183417/141‑1 (MQ=60)

                                                                                            GACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGGGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATACC                           >  SRR3722119.116356/1‑151 (MQ=60)

                                                                                            GACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATAC                            <  SRR3722119.402214/150‑1 (MQ=60)

                                                                                            GACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATAC                            <  SRR3722119.336511/150‑1 (MQ=60)

                                                                                            GACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATAC                            <  SRR3722119.51913/150‑1 (MQ=60)

                                                                                                     GTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATACCATGTACGAT                  >  SRR3722119.438926/1‑151 (MQ=60)

                                                                                                     GTATCAGGCGCTGTTTAATCTGGAAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATACCATGTACGAT                  <  SRR3722119.318136/151‑1 (MQ=60)

                                                                                                                     AATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATACCATGTACctgtctcttatacacatct  >  SRR3722119.134775/1‑132 (MQ=60)

GAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGAAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATACCATGTACGATCTGCCGTCCATGGAAG  >  NZ_CP009273/453788‑454055

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