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 I175 R1	328	64.7	2582366	96.2%	2484236	120.9

Breseq alignment

BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	493,309	(C)_5→4	intergenic (+34/‑202)	adk → / → hemH	adenylate kinase/ferrochelatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	493,305	0	C	.	97.7%	187.4 / ‑1.9	44	intergenic (+30/‑206)	adk/hemH	adenylate kinase/ferrochelatase
Reads supporting (aligned to +/- strand): ref base C (1/0); new base . (14/29); total (15/29)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

GGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACC  >  NZ_CP009273/493168‑493436

ggCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGGAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGc                                                                                                                                     <  1:543016/138‑1 (MQ=255)

ggCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGcc                                                                                                                                    <  2:569458/139‑1 (MQ=255)

ggCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGc                                                                                                                                     <  2:295752/138‑1 (MQ=255)

ggCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGc                                                                                                                                     <  1:999818/138‑1 (MQ=255)

  ctactaCTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGcccc                                                                                                                                  <  2:1277583/139‑1 (MQ=255)

  ctactaCTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGcccc                                                                                                                                  >  1:980681/1‑139 (MQ=255)

           cAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGga                                                                                                                          <  2:693667/137‑1 (MQ=255)

           cAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGga                                                                                                                          <  2:336276/137‑1 (MQ=255)

           cAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGga                                                                                                                          <  2:633241/137‑1 (MQ=255)

           cAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGga                                                                                                                          <  1:1277058/137‑1 (MQ=255)

           cAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGga                                                                                                                          <  1:974637/137‑1 (MQ=255)

           cAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGga                                                                                                                          >  1:56229/1‑137 (MQ=255)

           cAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGga                                                                                                                          <  1:755667/137‑1 (MQ=255)

           cAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGga                                                                                                                          <  1:627833/137‑1 (MQ=255)

                  gcagaagcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGga                                                                                                                          >  1:1119705/1‑130 (MQ=255)

                  gcagaagcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGga                                                                                                                          <  2:1119705/130‑1 (MQ=255)

                   cagaagcGGGTAATACCGAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTa                                                                                                                  >  1:468676/1‑137 (MQ=255)

                   cagaagcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTa                                                                                                                  <  2:468676/137‑1 (MQ=255)

                   cagaagcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAgg                                                                                                                <  1:1264202/139‑1 (MQ=255)

                     gaagcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGt                                                                                                              <  1:4021/139‑1 (MQ=255)

                        gcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAg                                                                                                           <  2:828198/139‑1 (MQ=255)

                        gcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAg                                                                                                           <  2:18669/139‑1 (MQ=255)

                             tAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                      <  1:699624/139‑1 (MQ=255)

                                     aaTACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGga                                                                                                                          <  2:509006/111‑1 (MQ=255)

                                     aaTACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGga                                                                                                                          >  1:509006/1‑111 (MQ=255)

                                     aaTACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAgg                                                                                                                <  1:597028/121‑1 (MQ=255)

                                     aaTACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAgg                                                                                                                >  2:597028/1‑121 (MQ=255)

                                     aaTACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCg                                                                                              <  2:980681/139‑1 (MQ=255)

                                     aaTACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAAGCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGAAGAGGGTTAGGGTGAGGGGAACAGGCCCg                                                                                              <  2:326475/139‑1 (MQ=255)

                                                   gACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACtt                                                                                 <  2:1205364/138‑1 (MQ=255)

                                                       gCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATGCGAAAGCGCGCACGGACAGCCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAg                                                                            <  1:1272814/139‑1 (MQ=255)

                                                         aCCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGc                                                                           >  1:51066/1‑138 (MQ=255)

                                                          ccAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCaa                                                                         >  2:921043/1‑139 (MQ=255)

                                                                 ggTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAgg                                                                  <  2:237673/139‑1 (MQ=255)

                                                                            ttCGCTCTGATCTTGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGAGAGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATattca                                                       >  2:1287555/1‑139 (MQ=38)

                                                                              cgcgCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCATCG‑CCCCTCGGGGAGAGGGCTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATattcatt                                                     >  2:961555/1‑139 (MQ=255)

                                                                                 gcTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAg                                                                                                           >  2:1189450/1‑82 (MQ=255)

                                                                                 gcTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAg                                                                                                           <  1:1189450/82‑1 (MQ=255)

                                                                                               aaTCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGgaga                                                                                                                        <  2:1242467/55‑1 (MQ=255)

                                                                                               aaTCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGgaga                                                                                                                        >  1:1242467/1‑55 (MQ=255)

                                                                                               aaTCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCa                                                                                      >  2:71312/1‑89 (MQ=255)

                                                                                               aaTCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCa                                                                                      <  1:71312/89‑1 (MQ=255)

                                                                                                        cTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTg                           <  2:490532/139‑1 (MQ=255)

                                                                                                                 aaGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTcc                  <  2:794740/139‑1 (MQ=255)

                                                                                                                      cgcACGGGCAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCttttt             <  1:1097042/139‑1 (MQ=255)

                                                                                                                      cgcACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCtttt              >  2:361156/1‑139 (MQ=255)

                                                                                                                      cgcACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCttttt             >  1:426776/1‑139 (MQ=255)

                                                                                                                      cgcACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCttttt             <  1:563110/139‑1 (MQ=255)

                                                                                                                              cAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTg     <  2:56229/139‑1 (MQ=255)

                                                                                                                                 tCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGAcc  >  2:1218947/1‑139 (MQ=255)

GGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACC  >  NZ_CP009273/493168‑493436

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

CCGCTGATCGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATT  >  NZ_CP009273/493159‑493449

CCGCTGATCGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCT                                                                                                                                              >  SRR3722129.994553/1‑151 (MQ=60)

        CGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGGAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAG                                                                                                                                      <  SRR3722129.550977/150‑1 (MQ=60)

        CGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAG                                                                                                                                      <  SRR3722129.1013912/150‑1 (MQ=60)

         GGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAG                                                                                                                                      >  SRR3722129.56963/1‑149 (MQ=60)

                ACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGG                                                                                                                             >  SRR3722129.1135325/1‑151 (MQ=60)

                 CTCCAAAGAAGCAGAAGCGGGTAATACCGAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGG                                                                                                                            >  SRR3722129.475489/1‑151 (MQ=60)

                   CCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGG                                                                                                                            <  SRR3722129.1294793/149‑1 (MQ=60)

                   CCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGG                                                                                                                            <  SRR3722129.988435/149‑1 (MQ=60)

                   CCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGG                                                                                                                            <  SRR3722129.766686/149‑1 (MQ=60)

                   CCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGG                                                                                                                            <  SRR3722129.637167/149‑1 (MQ=60)

                           agagacagGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAAC                                                                                                                  <  SRR3722129.605941/143‑1 (MQ=60)

                           GCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAAC                                                                                                                  <  SRR3722129.1281778/151‑1 (MQ=60)

                             AGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAG                                                                                                                <  SRR3722129.4075/151‑1 (MQ=60)

                                gagattcgccttagtctcgtgggctcggagatgtgtataagagacaGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCC                                                                                                             <  SRR3722129.1205999/105‑1 (MQ=60)

                                   GGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGctgtctcttatacacatc                                                                                                          >  SRR3722129.516461/1‑133 (MQ=60)

                                     GTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCAC                                                                                                        <  SRR3722129.709946/151‑1 (MQ=60)

                                                     gccttagtctcgtgggctcggagatgtgtataagagacagGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGC                                                                                        <  SRR3722129.72247/111‑1 (MQ=60)

                                                       AAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCA                                                                                       >  SRR3722129.51753/1‑150 (MQ=60)

                                                               GGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATGCGAAAGCGCGCACGGACAGCCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGC                                                                              <  SRR3722129.1290487/151‑1 (MQ=60)

                                                                                             GATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaa                                                 >  SRR3722129.1259712/1‑77 (MQ=60)

                                                                                                   GAACAAATGCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCCCG‑CCCCGCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATC                                          <  SRR3722129.466918/151‑1 (MQ=60)

                                                                                                                    ATTCGAAAGCGCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTC                         >  SRR3722129.432981/1‑151 (MQ=60)

                                                                                                                              GCGCACGGGCAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACC               <  SRR3722129.1112365/151‑1 (MQ=60)

                                                                                                                              GCGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACC               <  SRR3722129.571465/151‑1 (MQ=60)

                                                                                                                               CGCACGGACAGTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaaaaaaacggatgacgac               >  SRR3722129.577787/1‑61 (MQ=34)

                                                                                                                                         GTCCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCA     >  SRR3722129.193314/1‑150 (MQ=60)

                                                                                                                                           CCCCTCG‑CCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATT  >  SRR3722129.1061322/1‑151 (MQ=60)

CCGCTGATCGGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATT  >  NZ_CP009273/493159‑493449

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