BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,030,706	(C)_5→6	coding (396/696 nt)	yedJ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,030,701	1	.	C	97.2%	541.7 / 7.9	177	coding (401/696 nt)	yedJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (3/2); new base C (93/79); total (96/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.07e-01

Only 100 of 219 total aligned reads displayed.

AAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑G‑CCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTGCGCGCTGAAACTGTGAGCGG  >  NC_000913/2030610‑2030788

aaGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGGAAACACGCGCGCCAGGCCAATC‑G‑c                                                                                        >  1:3926283/1‑93 (MQ=255)

aaGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑G‑c                                                                                        >  1:3941687/1‑93 (MQ=255)

aaGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑G‑c                                                                                        >  2:642756/1‑93 (MQ=255)

aaGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCACGCCAATC‑G‑c                                                                                        >  1:3827265/1‑93 (MQ=255)

 aGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑G‑cc                                                                                       >  1:1056735/1‑93 (MQ=255)

 aGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑G‑cc                                                                                       >  2:1253644/1‑93 (MQ=255)

    gggCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑G‑ccccc                                                                                    <  1:5029123/93‑1 (MQ=255)

    gggCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑G‑ccccc                                                                                    >  2:161444/1‑93 (MQ=255)

      gCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑G‑cccccca                                                                                  <  2:4528603/93‑3 (MQ=255)

       cGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑G‑ccccccag                                                                                 <  2:101578/93‑4 (MQ=255)

       cGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑G‑ccccccag                                                                                 <  2:254915/93‑4 (MQ=255)

       cGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑G‑ccccccag                                                                                 <  1:1458926/93‑4 (MQ=255)

        gATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑G‑ccccccagg                                                                                <  2:4353165/93‑5 (MQ=255)

           gCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCt                                                                              <  2:797636/93‑1 (MQ=255)

            cTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCtt                                                                             >  1:6043419/1‑93 (MQ=255)

              gtgCGAACGGATCTTCACCATCAAGCAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTcc                                                                           <  2:3427823/93‑1 (MQ=255)

                 cGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGc                                                                        >  2:206745/1‑93 (MQ=255)

                  gAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGcc                                                                       >  1:41950/1‑93 (MQ=255)

                  gAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGcc                                                                       >  1:2689095/1‑93 (MQ=255)

                   aaCGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCg                                                                      >  1:5169476/1‑93 (MQ=255)

                   aaCGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCg                                                                      >  2:3355380/1‑93 (MQ=255)

                   aaCGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCg                                                                      >  1:2070962/1‑93 (MQ=255)

                   aaCGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCg                                                                      >  1:5400108/1‑93 (MQ=255)

                    aCGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGa                                                                     <  2:5690620/93‑1 (MQ=255)

                     cagATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGAt                                                                    >  2:2847944/3‑93 (MQ=255)

                     cGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGAt                                                                    >  2:4753191/1‑93 (MQ=255)

                      ggATCTTCACCATCAAACAGCGCCACGCCCAATTCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATc                                                                   <  2:248535/93‑1 (MQ=255)

                      ggATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATc                                                                   <  1:469267/93‑1 (MQ=255)

                      ggATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATc                                                                   >  1:156517/1‑93 (MQ=255)

                      cgATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATc                                                                   <  1:5232586/92‑1 (MQ=255)

                        aTCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCg                                                                 >  2:5076260/1‑93 (MQ=255)

                         tCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGc                                                                >  2:1343416/1‑93 (MQ=255)

                         tCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGc                                                                >  2:5233485/1‑93 (MQ=255)

                         tCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGc                                                                >  2:61740/1‑93 (MQ=255)

                         tCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGc                                                                >  1:1213231/1‑93 (MQ=255)

                          cTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCa                                                               >  2:5112594/1‑93 (MQ=255)

                             cACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATcc                                                            <  1:3105974/93‑1 (MQ=255)

                               ccATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTg                                                          <  2:4301735/93‑1 (MQ=255)

                               ccATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTg                                                          <  1:5633204/93‑1 (MQ=255)

                                cagcAAGCAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTgg                                                         <  2:1991613/90‑1 (MQ=255)

                                cATCAAACATCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCATGCCAATC‑GCCCCCCAGCGCTTCCAGCCGATCCGCATGCTgg                                                         >  1:1311258/1‑93 (MQ=255)

                                cATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGACAATC‑GCCCCCCAGGGCTTCCAGACGATTCGCATCCTgg                                                         >  1:840888/1‑93 (MQ=255)

                                 aTCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGa                                                        <  2:3083364/93‑1 (MQ=255)

                                  tCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGAc                                                       >  1:2616190/1‑93 (MQ=255)

                                       cAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAAttt                                                  <  2:5556780/93‑1 (MQ=255)

                                       cAGCGCCACGCCCAAGGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAAttt                                                  <  1:187026/93‑1 (MQ=255)

                                        aGCGCCACGCCCAATCCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAAtttt                                                 <  1:357855/93‑1 (MQ=255)

                                           gcCACTCCCAATTCCCCTGAAACGGCAACCACGCGCGCCAGGCCAATCAC‑CCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGc                                              <  1:4871463/93‑1 (MQ=255)

                                           gcCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGc                                              >  2:4021042/1‑93 (MQ=255)

                                           gcCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGc                                              >  1:2808839/1‑93 (MQ=255)

                                           gcCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGc                                              <  2:4259787/93‑1 (MQ=255)

                                              ccGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTc                                           <  1:2983668/92‑1 (MQ=255)

                                               cGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGACAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTcc                                          >  1:717576/1‑93 (MQ=255)

                                                   cAATGCCCCTGAAACGGCAAACACGCGCGCCAGGGCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCCTCCTGGACAATTTTAGCTTCCgttg                                      >  2:3130178/1‑93 (MQ=255)

                                                   cAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCgttg                                      >  1:1755221/1‑93 (MQ=255)

                                                    aaTGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCgttgt                                     <  2:839709/93‑1 (MQ=255)

                                                     aTGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCgttgtt                                    >  1:4468603/1‑93 (MQ=255)

                                                      tGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTa                                   >  2:3250183/1‑93 (MQ=255)

                                                       gCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTaa                                  <  1:3372080/93‑1 (MQ=255)

                                                       gCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTaa                                  <  2:5986368/93‑1 (MQ=255)

                                                         cccTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAgg                                <  1:5334187/93‑1 (MQ=255)

                                                         cccTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGACGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAgg                                >  1:1422449/1‑93 (MQ=255)

                                                          ccTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAggg                               >  2:602806/1‑93 (MQ=255)

                                                          ccTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAggg                               <  1:2519461/93‑1 (MQ=255)

                                                           cTGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGATTCCCGCCGCTCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAgggg                              >  2:2611311/1‑93 (MQ=255)

                                                            tGAAACGGCAACCACGCGCGCCAGCCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAggggg                             <  1:4103929/93‑1 (MQ=255)

                                                            tGAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAggggg                             <  2:584129/93‑1 (MQ=255)

                                                             gAAACGGCAACCCCGCGCGCCAGCCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGc                            <  2:2239808/93‑1 (MQ=255)

                                                             gAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGc                            <  2:4046238/93‑1 (MQ=255)

                                                             gAAACGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGc                            <  1:4003437/93‑1 (MQ=255)

                                                              aaaCGGCAAACACGCGCGCCATGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCg                           >  1:4951253/1‑93 (MQ=255)

                                                              aaaCGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCg                           >  2:5312880/1‑93 (MQ=255)

                                                                aCGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGAt                         <  1:3577152/93‑1 (MQ=255)

                                                                aCGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGAt                         <  2:1803504/93‑1 (MQ=255)

                                                                aCGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGAt                         <  2:4891977/93‑1 (MQ=255)

                                                                aCGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGAt                         <  1:946623/93‑1 (MQ=255)

                                                                 cGGCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGAtt                        <  1:6219149/93‑1 (MQ=255)

                                                                  ggCAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGAttt                       >  1:1146747/1‑93 (MQ=255)

                                                                    cAAACACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTgc                     <  1:5908333/93‑1 (MQ=255)

                                                                     aaaCACGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTgcg                    >  2:4791940/1‑93 (MQ=255)

                                                                        cacGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTgcgcgc                 >  2:397877/1‑93 (MQ=255)

                                                                         acGCGCGCCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTGCGCGCt                <  1:439833/93‑1 (MQ=255)

                                                                           ccgcccCAGGCCAACC‑GCCCCCCAGTGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTGCGCGCTGa              <  2:3703209/92‑1 (MQ=255)

                                                                             gcgcCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTGCGCGCTGaaa            <  2:6100089/93‑1 (MQ=255)

                                                                              cgcCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTGCGCGCTGAAAc           <  2:1949505/93‑1 (MQ=255)

                                                                              cgcCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTGCGCGCTGAAAc           <  1:2141066/93‑1 (MQ=255)

                                                                              cgcCAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTGCGCGCTGAAAc           <  1:2050508/93‑1 (MQ=255)

                                                                                ccAGGCGATTC‑GCCCCCCAGGCCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTGCGCGCTGAAACtg         <  2:1574122/93‑1 (MQ=255)

                                                                                ccAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTGCGCGCTGAAACtg         >  2:3157616/1‑93 (MQ=255)

                                                                                 cAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTGCGCGCTGAAACtgt        >  2:92289/1‑93 (MQ=255)

                                                                                 cAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTGCGCGCTGAAACtgt        <  2:277175/93‑1 (MQ=255)

                                                                                 cAGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTAGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTGCGCGCTGAAACtgt        >  1:3382725/1‑93 (MQ=255)

                                                                                  aGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTGCGCGCTGAAACtgtg       >  2:4619118/1‑93 (MQ=255)

                                                                                  aGGCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTGCGCGCTGAAACtgtg       >  2:151174/1‑93 (MQ=255)

                                                                                   ggCCAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTGCGCGCTGAAACTGTGa      >  1:2373228/1‑93 (MQ=255)

                                                                                     ccAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTGCGCGCTGAAACTGTGAGc    <  2:3833999/93‑1 (MQ=255)

                                                                                     ccAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTGCGCGCTGAAACTGTGAGc    <  1:3333699/93‑1 (MQ=255)

                                                                                      cAATC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTGCGCGCTGAAACTGTGAGCg   >  2:4655299/1‑93 (MQ=255)

                                                                                       aaTC‑GCCCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTGCGCGCTGAAACTGTGAGCgg  <  2:3062408/93‑1 (MQ=255)

                                                                                       aaTC‑GCCCCCCAGGGCCTCCAGCCGATCAGAATCCTGGACAATTTTAGCCTCCGTTGTTTAGGGGGCGAGCTGCGAGCTGAGACTATGAGCgg  >  2:1733870/1‑93 (MQ=255)

AAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCGCCAGGCCAATC‑G‑CCCCCAGGGCTTCCAGCCGATCCGCATCCTGGACAATTTTAGCTTCCGTTGTTAAGGGGGCGATTTGCGCGCTGAAACTGTGAGCGG  >  NC_000913/2030610‑2030788

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 16 ≤` `ATCG/ATCG` `< 24 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`