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	4,589,054	T→G	T69P (ACC→CCC)	opgB ←	phosphatidylglycerol‑‑membrane‑oligosaccharide glycerophosphotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,589,054	0	T	G	97.9%	167.0 / ‑2.9	47	coding (205/2292 nt)	opgB	phosphatidylglycerol‑‑membrane‑oligosaccharide glycerophosphotransferase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base G (19/27); minor base . (0/1); total (19/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

CTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAG  >  NZ_CP009273/4588927‑4589184

cTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTgg                                                                                                                         <  1:406910/139‑1 (MQ=255)

   tAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGtt                                                                                                                       <  2:482487/138‑1 (MQ=255)

   tAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTa                                                                                                                      <  1:289654/139‑1 (MQ=255)

   tAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTa                                                                                                                      <  1:451392/139‑1 (MQ=255)

   tAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTa                                                                                                                      <  2:702561/139‑1 (MQ=255)

   tAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTa                                                                                                                      <  1:616960/139‑1 (MQ=255)

     aCCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAg                                                                                                                    <  1:498534/139‑1 (MQ=255)

        aaaaGGGTGCGGAGGATGGCGACGACGGCGCAGGATCCAGCCCGGCGCGCCGAGCACCGCTGTCGGCCCCAGCGAAATCCGGATGCCCGGCAAAATTTATTTGCTGGCGCCAGCCGCGGGTGGGCTTTTGGTAAAGGta                                                                                                                 <  2:250480/139‑1 (MQ=255)

               tGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGaa                                                                                                          >  2:588618/1‑139 (MQ=255)

                    atgatgGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACcgcg                                                                                                     <  1:310599/139‑1 (MQ=255)

                              cgacgGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATAc                                                                                           >  1:736390/1‑139 (MQ=255)

                              cgacgGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATAc                                                                                           >  1:4864/1‑139 (MQ=255)

                                                cccagcGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGt                                                                         <  2:532009/139‑1 (MQ=255)

                                                  cagcGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCg                                                                       <  1:621202/139‑1 (MQ=255)

                                                     cgcACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAg                                                                                                                    >  2:303220/1‑91 (MQ=255)

                                                     cgcACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAg                                                                                                                    <  1:303220/91‑1 (MQ=255)

                                                     cgcACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTg                                                                    >  2:632825/1‑139 (MQ=255)

                                                     cgcACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTg                                                                    >  1:86248/1‑139 (MQ=255)

                                                          cGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCaaa                                                               <  2:4864/139‑1 (MQ=255)

                                                                  gCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACcgc                                                                                                      <  1:228880/92‑1 (MQ=255)

                                                                  gCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACcgc                                                                                                      >  2:228880/1‑92 (MQ=255)

                                                                        aGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATAttt                                                  >  2:36981/1‑138 (MQ=255)

                                                                              aGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGa                                           <  1:632825/139‑1 (MQ=255)

                                                                                   aaTGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAAt                                      <  1:696870/139‑1 (MQ=255)

                                                                                    aTGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACaaa                                       >  2:696506/1‑137 (MQ=255)

                                                                                    aTGCCGATACCCGACAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAg                                                                                                                    >  2:342560/1‑60 (MQ=255)

                                                                                    aTGCCGATACCCGACAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAg                                                                                                                    <  1:342560/60‑1 (MQ=255)

                                                                                     tGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGcc                                                                                    >  2:488237/1‑91 (MQ=255)

                                                                                     tGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGcc                                                                                    <  1:488237/91‑1 (MQ=255)

                                                                                     tGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTg                                                        >  2:160988/1‑119 (MQ=255)

                                                                                     tGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTg                                                        <  1:160988/119‑1 (MQ=255)

                                                                                     tGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAg                                    >  2:148680/1‑139 (MQ=255)

                                                                                     tGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAg                                    <  1:612216/139‑1 (MQ=255)

                                                                                       ccGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCAGCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGcc                                  <  2:728736/139‑1 (MQ=255)

                                                                                               cGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAg                                                                                                                    >  1:220533/1‑49 (MQ=255)

                                                                                               cGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAg                                                                                                                    <  2:220533/49‑1 (MQ=255)

                                                                                                 gCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCt                                                                     <  2:249932/94‑1 (MQ=255)

                                                                                                 gCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCt                                                                     >  1:249932/1‑94 (MQ=255)

                                                                                                 gCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGt                        >  2:520954/1‑139 (MQ=255)

                                                                                                 gCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGt                        >  1:519235/1‑139 (MQ=255)

                                                                                                           tttGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGc               <  2:426187/138‑1 (MQ=255)

                                                                                                               cTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCa           <  1:696506/138‑1 (MQ=255)

                                                                                                                   cGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAacca      >  2:506345/1‑139 (MQ=255)

                                                                                                                   cGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAacca      >  1:464422/1‑139 (MQ=255)

                                                                                                                     ccAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAaccacc    <  2:384530/139‑1 (MQ=255)

                                                                                                                     ccAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAaccacc    <  2:103696/139‑1 (MQ=255)

                                                                                                                      aggCACCGG‑CAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAg  <  2:678032/137‑1 (MQ=255)

CTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAG  >  NZ_CP009273/4588927‑4589184

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

GCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACG  >  NZ_CP009273/4588926‑4589192

GCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAG                                                                                                                      <  SRR3722133.412537/151‑1 (MQ=60)

   GTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGA                                                                                                                   <  SRR3722133.625077/151‑1 (MQ=60)

   GTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGA                                                                                                                   <  SRR3722133.293592/151‑1 (MQ=60)

   GTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGA                                                                                                                   <  SRR3722133.457525/151‑1 (MQ=60)

     AACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAAC                                                                                                                 <  SRR3722133.505264/151‑1 (MQ=60)

      tccaagcagaagacggcatacgagatctagtacggtctcgggggctcggagatgtgtataagagacagGCCCCAGCACAATGCCGATACCCGACAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAAC                                                                                                                 <  SRR3722133.347359/82‑1 (MQ=60)

      acggtctcgtgggctcggagatgtgtataagagacaGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAAC                                                                                                                 <  SRR3722133.307382/114‑1 (MQ=60)

                   cagcgtggtcgggatggccttttttttttttcaagcagaagacggcatacgatatctagtacggtctcgtgggctcggagatgtgtataagagacagCGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCTTCGTTAATAC                                                                                                   <  SRR3722133.479717/54‑1 (MQ=60)

                   acggtctcgtgggctcggagatgtgtataagagacaGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATAC                                                                                                   <  SRR3722133.232003/115‑1 (MQ=60)

                    GATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACC                                                                                                  >  SRR3722133.745910/1‑151 (MQ=60)

                    GATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACC                                                                                                  <  SRR3722133.314866/151‑1 (MQ=60)

                    GATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACC                                                                                                  >  SRR3722133.4931/1‑151 (MQ=60)

                                      acggtctcgtgggctcggagatgtgtataagagaCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGT                                                                                 <  SRR3722133.494821/116‑1 (MQ=60)

                                           ATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGG                                                                           >  SRR3722133.87369/1‑151 (MQ=60)

                                                  CCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAG                                                                    <  SRR3722133.629352/151‑1 (MQ=60)

                                                                  aagagaCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAAC                                                     <  SRR3722133.163089/144‑1 (MQ=60)

                                                                              CAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCC                                        <  SRR3722133.641116/151‑1 (MQ=60)

                                                                                   CAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCAC                                   <  SRR3722133.705900/151‑1 (MQ=60)

                                                                                     ATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCtgtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgccgtcttctgcttgaaaaaaaaaaaat                                  >  SRR3722133.223544/1‑72 (MQ=60)

                                                                                     ATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCG                                 <  SRR3722133.620295/151‑1 (MQ=60)

                                                                                       GCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGctgtctcttatacacatctccgagcccacgagacc                               >  SRR3722133.253343/1‑116 (MQ=60)

                                                                                       GCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTT                               >  SRR3722133.526236/1‑151 (MQ=60)

                                                                                                         GTATTTGCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCAC             >  SRR3722133.470723/1‑151 (MQ=60)

                                                                                                               GCTGACGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGT        <  SRR3722133.705532/150‑1 (MQ=60)

                                                                                                                    CGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACG  >  SRR3722133.323755/1‑151 (MQ=60)

                                                                                                                    CGCCAGCACCGGGCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACG  >  SRR3722133.298461/1‑151 (MQ=60)

GCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACG  >  NZ_CP009273/4588926‑4589192

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