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 I176 R1	180	42.9	1710034	96.5%	1650182	119.3

Breseq alignment

BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	449,442	A→G	V61A (GTC→GCC)	yajG ←	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	449,442	0	A	G	100.0%	145.3 / NA	39	V61A (GTC→GCC)	yajG	lipoprotein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (14/25); total (14/25)

ATGATTTGCAGATTAACCGGGCCATTCGGACCAACCATGTAACCACGCGCGGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATCCTGCTGTGGCAGTGTAATCGTCGGGGAAACTTCAATAGTTGTTGGCGGTTTTGCGCATCCTG  >  NZ_CP009273/449307‑449577

atgatTTGCAGATTAACCGGGCCATTCGGACCAACCATGTAACCACGCGCGGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCtt                                                                                                                                      >  1:526671/1‑139 (MQ=255)

      tGCAGATTAACCGGGCCATTCGGACCAACCATGTAACCACGCGCGGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGcc                                                                                                                                 <  2:93442/138‑1 (MQ=255)

            ttAACCGGGCCATTCGGACCAACCATGTAACCACGCGCGGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCt                                                                                                                          <  2:620722/139‑1 (MQ=255)

               aCCGGGCCATTCGGACCAACCATGTAACCACGCGCGGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGAt                                                                                                                       <  2:787205/139‑1 (MQ=255)

               aCCGGGCCATTCGGACCAACCATGTAACCACGCGCGGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGAt                                                                                                                       <  2:823978/139‑1 (MQ=255)

                      cATTCGGACCAACCATGTAACCACGCGCGGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACg                                                                                                                <  1:351405/139‑1 (MQ=255)

                         tCGGACCAACCATGTAACCACGCGCGGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGCTTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCt                                                                                                              >  2:683282/1‑138 (MQ=255)

                               caaccaTGTAACCACGCGCGGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGc                                                                                                       <  2:526671/139‑1 (MQ=255)

                                  ccaTGTAACCACGCGCGGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCAc                                                                                                     >  1:595642/1‑138 (MQ=255)

                                           cACGCGCGGTCATCTGTTTTTCCAGCCCTTCTTGCAGCAGGAAACGCCGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCt                                                                                           <  1:395268/139‑1 (MQ=255)

                                            aCGCGCGGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCt                                                                                           >  2:364051/1‑138 (MQ=255)

                                             cgcgcgGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTg                                                                                                             >  1:301571/1‑119 (MQ=255)

                                             cgcgcgGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTg                                                                                                             <  2:301571/119‑1 (MQ=255)

                                             cgcgcgGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTa                                                                                         >  1:291351/1‑139 (MQ=255)

                                             cgcgcgGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTa                                                                                         >  1:437885/1‑139 (MQ=255)

                                                   gTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGa                                                                                   <  2:765115/139‑1 (MQ=255)

                                                           tttttCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCAt                                                                            <  2:80854/138‑1 (MQ=255)

                                                           tttttCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCAt                                                                            <  2:753318/138‑1 (MQ=255)

                                                           tttttCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCAt                                                                            <  2:486838/138‑1 (MQ=255)

                                                           tttttCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCAt                                                                            <  2:621605/138‑1 (MQ=255)

                                                                ccAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCt                                                                      >  1:295826/1‑139 (MQ=255)

                                                                   gCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTg                                                                    >  2:543723/1‑138 (MQ=255)

                                                                     aCTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGac                                                                 <  1:543723/139‑2 (MQ=255)

                                                                     aCTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGac                                                                 <  1:364051/139‑2 (MQ=255)

                                                                     aCTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGAt                                                                 <  2:437885/139‑1 (MQ=255)

                                                                                   gAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATCCTGCTGTGGCAgt                                                   >  2:744491/1‑139 (MQ=255)

                                                                                       cGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATCCTGCTGTGGCAGTGTaa                                               <  1:683282/139‑1 (MQ=255)

                                                                                           gATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATCCTGCTGTGGCAGTGTAAtcgt                                           <  2:362085/139‑1 (MQ=255)

                                                                                             tCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGACCCTGCTGTGGCAGTGTAAtcgtcg                                         >  1:662493/1‑139 (MQ=255)

                                                                                                         gTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGACCCTGCTGTGGCAGTGTAATCGTCGGGGAAACTTCaa                             <  2:662493/139‑1 (MQ=255)

                                                                                                                   cGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATCCTGCTGTGGCAGTGTAATCGTCGGGGAAACTTCAATAGTTGTTgg                   >  1:625381/1‑139 (MQ=255)

                                                                                                                        tGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAgg                                                                        >  1:174850/1‑81 (MQ=255)

                                                                                                                        tGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAgg                                                                        <  2:174850/81‑1 (MQ=255)

                                                                                                                        tGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGc                                                                       >  2:96652/1‑82 (MQ=255)

                                                                                                                        tGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGc                                                                       <  1:96652/82‑1 (MQ=255)

                                                                                                                              tCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATCCTGCTGTGGCAGTGTAATCGTCGGGGAAACTTCAATAGTTGTTGGCGGTTTTgcg         <  1:461901/138‑1 (MQ=255)

                                                                                                                                    gTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATCCTGCTGTGGCAGTGTAATCGTCGGGGAAACTTCAATAGTTGTTGGCGGTTTTGCGCATCCTg  <  2:154567/139‑1 (MQ=255)

                                                                                                                                    gTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATCCTGCTGTGGCAGTGTAATCGTCGGGGAAACTTCAATAGTTGTTGGCGGTTTTGCGCATCCTg  <  1:48837/139‑1 (MQ=255)

                                                                                                                                    gTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGACCCTGCTGTGGCAGTGTAATCGTCGGGGAAACTTCAATAGTTGTTGGCGGTTTTGCGCATCCTg  <  2:295826/139‑1 (MQ=255)

ATGATTTGCAGATTAACCGGGCCATTCGGACCAACCATGTAACCACGCGCGGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATCCTGCTGTGGCAGTGTAATCGTCGGGGAAACTTCAATAGTTGTTGGCGGTTTTGCGCATCCTG  >  NZ_CP009273/449307‑449577

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

GTTGGCTAACGATGATTTGCAGATTAACCGGGCCATTCGGACCAACCATGTAACCACGCGCGGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATCCTGCTGTGGCAGTGTAATCGTCGGGGAAACTTCAATAGTTGTTGGCGGTTTTGCGCATCCTGCAAGCATAAAC  >  NZ_CP009273/449296‑449588

GTTGGCTAACGATGATTTGCAGATTAACCGGGCCATTCGGACCAACCATGTAACCACGCGCGGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTT                                                                                                                                                >  SRR3722130.534727/1‑151 (MQ=60)

                                CCATTCGGACCAACCATGTAACCACGCGCGGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCAC                                                                                                                <  SRR3722130.356865/151‑1 (MQ=60)

                                  ATTCGGACCAACCATGTAACCACGCGCGGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACC                                                                                                               >  SRR3722130.604602/1‑150 (MQ=60)

                                             CCATGTAACCACGCGCGGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATctgtctctt                                                                                                    >  SRR3722130.306243/1‑141 (MQ=60)

                                             CCATGTAACCACGCGCGGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTAC                                                                                                   >  SRR3722130.295886/1‑151 (MQ=60)

                                             CCATGTAACCACGCGCGGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTAC                                                                                                   >  SRR3722130.444762/1‑151 (MQ=60)

                                                     CCACGCGCGGTCATCTGTTTTTCCAGCCCTTCTTGCAGCAGGAAACGCCGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGC                                                                                           <  SRR3722130.401552/151‑1 (MQ=60)

                                                                CATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTT                                                                                >  SRR3722130.300418/1‑151 (MQ=60)

                                                                         gatattcgccttagtctcgtgggctcggagatgtgtataagagacagGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGACCCTGC                                                                       <  SRR3722130.97986/104‑1 (MQ=60)

                                                                               CACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGACCCTGCTGTGGC                                                                 <  SRR3722130.551980/151‑1 (MQ=60)

                                                                               CACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGACCCTGCTGTGGC                                                                 <  SRR3722130.369752/151‑1 (MQ=60)

                                                                                             GGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGACCCTGCTGTGGCAGTGTAATCGTCGG                                                   >  SRR3722130.672377/1‑151 (MQ=60)

                                                                                                 ACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATCCTGCTGTGGCAGTGTAATCGTCGGGGAA                                               <  SRR3722130.693429/151‑1 (MQ=60)

                                                                                                                   GGTCAGGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATCCTGCTGTGGCAGTGTAATCGTCGGGGAAACTTCAATAGTTGTTGGC                             >  SRR3722130.634754/1‑151 (MQ=60)

                                                                                                                        GGGTAACGATTTGATTATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATCCTGCTGTctcttatacacatctccgagcccacgagactaaggcgaatctcg                        >  SRR3722130.177479/1‑107 (MQ=60)

                                                                                                                                        ATCGCGGGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATCCTGCTGTGGCAGTGTAATCGTCGGGGAAACTTCAATAGTTGTTGGCGGTTTTGCGCATCCTGCAAG         <  SRR3722130.469107/150‑1 (MQ=60)

                                                                                                                                              GGTGGCTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATCCTGCTGTGGCAGTGTAATCGTCGGGGAAACTTCAATAGTTGTTGGCGGTTTTGCGCATCCTGCAAGCATAAAC  <  SRR3722130.49503/151‑1 (MQ=60)

GTTGGCTAACGATGATTTGCAGATTAACCGGGCCATTCGGACCAACCATGTAACCACGCGCGGTCATCTGTTTTTCCAGCACTTCTTGCAGCAGGAAACGCAGATCGCGGGAGGCGGTCAGGGTAACGATTTGATTATCGCGGGTGACTTTTGCCAGCGCCTGATCGGTACGCTGATCGGCACCATTAATGCTTACGGTGACGCCCATCAGGCTTGGATCCTGCTGTGGCAGTGTAATCGTCGGGGAAACTTCAATAGTTGTTGGCGGTTTTGCGCATCCTGCAAGCATAAAC  >  NZ_CP009273/449296‑449588

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