BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,679,306	T→C	S89G (AGC→GGC)	yhjJ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,679,306	0	T	C	97.0%	511.8 / 3.4	168	S89G (AGC→GGC)	yhjJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/4); new base C (91/72); total (92/76)
Fisher's exact test for biased strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 174 total aligned reads displayed.

CATAAGAGACAATTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAA  >  NC_000913/3679216‑3679395

  tAAGAGACAATTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGc                                                                                       >  1:4295464/1‑93 (MQ=255)

   aaGAGACAATTTCCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCg                                                                                      <  1:3006064/93‑1 (MQ=255)

   aaGAGACAATTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCg                                                                                      >  2:3844805/1‑93 (MQ=255)

   aaGAGACAATTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCg                                                                                      <  2:4325917/93‑1 (MQ=255)

     gagaCAATTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGtt                                                                                    <  1:696143/93‑1 (MQ=255)

      agaCAATTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTa                                                                                   <  2:2965440/93‑1 (MQ=255)

      agaCAATTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGGTa                                                                                   >  2:2248564/1‑93 (MQ=255)

        aCAATTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAgc                                                                                 >  1:3710317/1‑93 (MQ=255)

           aTTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCa                                                                              >  2:2321274/1‑93 (MQ=255)

           aTTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCa                                                                              >  1:2940575/1‑93 (MQ=255)

            ttACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCaa                                                                             <  1:398856/93‑1 (MQ=255)

              aCCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATa                                                                           <  1:4672453/93‑1 (MQ=255)

              aCCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATa                                                                           >  1:3037730/1‑93 (MQ=255)

               ccGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATAc                                                                          >  2:3852472/1‑93 (MQ=255)

                   cggcATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGACAAGGCCACCGCCTTGCGTTAGCGCAATACGAgg                                                                      <  2:365688/93‑1 (MQ=255)

                    ggcATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAggg                                                                     <  2:2823906/93‑1 (MQ=255)

                    ggcATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAggg                                                                     >  1:3382153/1‑93 (MQ=255)

                     gcATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGa                                                                    >  2:2850462/1‑93 (MQ=255)

                     gcATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGa                                                                    >  1:1487863/1‑93 (MQ=255)

                       aTCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATg                                                                  >  2:4260692/1‑93 (MQ=255)

                        tCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATgg                                                                 >  1:2614828/1‑93 (MQ=255)

                        tCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATgg                                                                 <  2:4013636/93‑1 (MQ=255)

                        tCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATgg                                                                 >  2:3584424/1‑93 (MQ=255)

                        tCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATgg                                                                 >  2:1740383/1‑93 (MQ=255)

                           ggCTTTTGGGGTCGACCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGt                                                              <  2:2118330/93‑1 (MQ=255)

                             cGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGa                                                            >  2:3542410/1‑93 (MQ=255)

                              gTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGAc                                                           >  2:3151454/1‑93 (MQ=255)

                               tttAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACt                                                          >  2:4177211/1‑93 (MQ=255)

                                 tAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGt                                                        <  2:1676688/93‑1 (MQ=255)

                                 tAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGt                                                        <  2:3074956/93‑1 (MQ=255)

                                    ggTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAAc                                                     <  1:1815151/93‑1 (MQ=255)

                                    ggTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAAc                                                     <  1:4656189/93‑1 (MQ=255)

                                      tCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCg                                                   >  2:1613330/1‑93 (MQ=255)

                                       cGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGc                                                  <  1:92882/93‑1 (MQ=255)

                                         aTCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGctt                                                >  2:1323157/1‑92 (MQ=255)

                                         aTCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGctc                                                <  2:346342/93‑1 (MQ=255)

                                         aTCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGctc                                                >  1:4392294/1‑93 (MQ=255)

                                         aTCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGctc                                                >  1:1614642/1‑93 (MQ=255)

                                          tCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGctct                                               >  2:1962485/1‑93 (MQ=255)

                                          tCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGctct                                               >  1:1257682/1‑93 (MQ=255)

                                          tCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGctct                                               >  2:1784441/1‑93 (MQ=255)

                                          tCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGctct                                               >  1:870742/1‑93 (MQ=255)

                                          tCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGctct                                               >  1:1258676/1‑93 (MQ=255)

                                           ccccTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCtg                                              >  1:16703/1‑93 (MQ=255)

                                             ccTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCtgtt                                            >  2:471164/1‑93 (MQ=255)

                                               tgctgcCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCtgttgt                                          <  2:2025326/93‑1 (MQ=255)

                                               cgctgcCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCtgttgt                                          <  1:1294336/92‑1 (MQ=255)

                                                 ctgcCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTgtgt                                        <  1:4523126/93‑1 (MQ=255)

                                                  tgcCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGCCTGTAACCGCTCTGTTGTGTa                                       >  2:1293391/1‑93 (MQ=255)

                                                  tgcCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTa                                       >  2:1332503/1‑93 (MQ=255)

                                                  tgcCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTa                                       >  2:4590064/1‑93 (MQ=255)

                                                   gcCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTAc                                      <  1:750347/93‑1 (MQ=255)

                                                    gcacaACGAACGCGCCTGTGCTGCTTCAAGGCCACCGCCTTGTGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGATCTGTTGTGGACt                                     <  2:3637673/92‑1 (MQ=255)

                                                    ccACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTTCt                                     >  2:2522637/1‑93 (MQ=255)

                                                     cacaATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACtt                                    <  2:375652/93‑1 (MQ=255)

                                                     aacaATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACtt                                    <  1:2783210/92‑1 (MQ=255)

                                                      acaATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACttt                                   >  1:4028414/1‑93 (MQ=255)

                                                       caATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTc                                  <  2:3833196/93‑1 (MQ=255)

                                                       caATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTc                                  <  1:101817/93‑1 (MQ=255)

                                                       aaaTGAACGCGCCTGTGCTGCGTCAAGGCCACAGCCTTGCGTTAGCGCAATACGAGGGATGGCGTAACTGTAACCGCTCTGTTGTGTACTTTc                                  <  2:4448925/92‑1 (MQ=255)

                                                        aaTGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCg                                 >  1:4609808/1‑93 (MQ=255)

                                                        aaTGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCg                                 >  1:3555988/1‑93 (MQ=255)

                                                        aaTGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCg                                 <  2:3765960/93‑1 (MQ=255)

                                                        aaTGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCg                                 <  2:333534/93‑1 (MQ=255)

                                                            aaCGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGgcga                             >  1:3443921/1‑93 (MQ=255)

                                                            aaCGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGgcga                             <  1:4076558/93‑1 (MQ=255)

                                                             aCGCTCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGgcgag                            <  2:1109041/93‑1 (MQ=255)

                                                             aCGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGgcgag                            >  1:3596638/1‑93 (MQ=255)

                                                             aCGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGgcgag                            <  2:4737718/93‑1 (MQ=255)

                                                             aCGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGgcgag                            <  2:329898/93‑1 (MQ=255)

                                                             aCGCGCCTCTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTAACTGTAACCGCCCTGTTGTGTACTTTCGgcgag                            <  2:1351871/93‑1 (MQ=255)

                                                              cgcgCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGgcgagc                           <  2:2156904/93‑1 (MQ=255)

                                                               gcgcCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGgcgagcg                          >  1:3229957/1‑93 (MQ=255)

                                                                  ccTGTGCTGCGTCACGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAAc                       <  1:3600455/93‑1 (MQ=255)

                                                                  ccTGTGCTGCGTCAAGGCCACCGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGgcgagcgatc                       <  2:693447/93‑3 (MQ=255)

                                                                    tgtgCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGTGCGAACCg                     >  1:4733341/1‑93 (MQ=255)

                                                                    tgtgCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCg                     >  1:4182072/1‑93 (MQ=255)

                                                                     gtgCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCgg                    <  2:579188/93‑1 (MQ=255)

                                                                       gctgcGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTa                  >  2:1694282/1‑93 (MQ=255)

                                                                       gctgcGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTa                  <  2:1161977/93‑1 (MQ=255)

                                                                         tgcGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTAtt                >  1:256881/1‑93 (MQ=255)

                                                                            gTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAAc             >  2:2865197/1‑93 (MQ=255)

                                                                            gTCAAGGCCACCGCCGTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAAc             <  2:3839311/93‑1 (MQ=255)

                                                                             tCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAAcc            >  2:2823316/1‑93 (MQ=255)

                                                                             tCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAAcc            <  1:331302/93‑1 (MQ=255)

                                                                              cAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACca           <  1:502847/93‑1 (MQ=255)

                                                                              cAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACca           >  2:44518/1‑93 (MQ=255)

                                                                              cAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACca           >  2:3430403/1‑93 (MQ=255)

                                                                              cAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACca           >  1:2384630/1‑93 (MQ=255)

                                                                              cAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACca           <  2:2358382/93‑1 (MQ=255)

                                                                                aGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACcagc         >  2:4414107/1‑93 (MQ=255)

                                                                                 ggCCACCGCCTTGCGTTAGCGCAATACGCGGGATGGCGTGACTGTCACCGCTCTGTTGTGTACTTTCGGCGAGCGAACAGGTTTTAACcagcc        <  1:4207385/93‑2 (MQ=255)

                                                                                 ggCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACcagca        >  1:4560680/1‑93 (MQ=255)

                                                                                 ggCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACcagca        <  1:620768/93‑1 (MQ=255)

                                                                                    cACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGc     >  1:3140517/1‑93 (MQ=255)

                                                                                    cACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGc     >  1:3971939/1‑93 (MQ=255)

                                                                                    cACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGc     >  1:1196331/1‑93 (MQ=255)

                                                                                      ccGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGa   <  2:4721663/93‑1 (MQ=255)

                                                                                       cGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGaa  <  2:3809696/93‑1 (MQ=255)

                                                                                       cGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGaa  >  1:1056094/1‑93 (MQ=255)

CATAAGAGACAATTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAA  >  NC_000913/3679216‑3679395

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