BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	308,935	G→A	R106C (CGC→TGC)	yagY ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	308,935	0	G	A	99.4%	610.7 / NA	180	R106C (CGC→TGC)	yagY	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (83/95); minor base C (0/1); total (83/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 181 total aligned reads displayed.

GCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGC  >  NC_000913/308847‑309023

gcggAACCAGTGTCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCActct                                                                                      <  2:1658218/93‑1 (MQ=255)

 cggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCtt                                                                                     >  1:2434321/1‑93 (MQ=255)

  ggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCttt                                                                                    >  2:1155631/1‑93 (MQ=255)

  ggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCttt                                                                                    >  1:1956984/1‑93 (MQ=255)

   gAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCtttt                                                                                   <  1:452567/93‑1 (MQ=255)

   gAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCtttt                                                                                   <  2:3961630/93‑1 (MQ=255)

   gAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCtttt                                                                                   >  1:1600882/1‑93 (MQ=255)

   gAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCtttt                                                                                   >  2:2740561/1‑93 (MQ=255)

     aGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTc                                                                                 >  1:496300/1‑93 (MQ=255)

     aGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTc                                                                                 >  1:3255360/1‑93 (MQ=255)

     aGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTc                                                                                 <  2:3220816/93‑1 (MQ=255)

      gCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACGCTTTTTCa                                                                                <  2:4581264/93‑1 (MQ=255)

        aGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATc                                                                              <  1:2285575/93‑1 (MQ=255)

        aGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATc                                                                              <  1:2285576/93‑1 (MQ=255)

        aGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATc                                                                              <  1:2524889/93‑1 (MQ=255)

         gTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCt                                                                             >  2:2934572/1‑93 (MQ=255)

          tGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTg                                                                            <  1:217708/93‑1 (MQ=255)

           ggCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGc                                                                           <  2:3531534/93‑1 (MQ=255)

            gCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCg                                                                          <  1:849298/93‑1 (MQ=255)

             ccACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCgg                                                                         <  1:2812985/93‑1 (MQ=255)

              cacaGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCggg                                                                        <  1:651662/93‑1 (MQ=255)

               acaGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCgggc                                                                       >  1:3155325/1‑92 (MQ=255)

               acaGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGt                                                                       <  2:3435087/93‑1 (MQ=255)

               acaGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGt                                                                       >  2:3980645/1‑93 (MQ=255)

                caGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTc                                                                      >  1:3532647/1‑93 (MQ=255)

                    gcTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCtt                                                                  >  1:3923064/1‑93 (MQ=255)

                     cTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTta                                                                 <  2:3680396/93‑1 (MQ=255)

                     cTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTta                                                                 <  1:4272931/93‑1 (MQ=255)

                     cTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTta                                                                 >  1:4111432/1‑93 (MQ=255)

                      tGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTtat                                                                >  2:2608279/1‑93 (MQ=255)

                       gCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTtata                                                               <  2:2261137/93‑1 (MQ=255)

                       gCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTtata                                                               <  2:1409305/93‑1 (MQ=255)

                        cGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATag                                                              <  1:974052/93‑1 (MQ=255)

                          gTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATagaa                                                            <  1:1742792/93‑1 (MQ=255)

                          gTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATagaa                                                            <  2:1126242/93‑1 (MQ=255)

                          gTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATagaa                                                            >  2:2045361/1‑93 (MQ=255)

                           ttGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATagaag                                                           <  2:1629186/93‑1 (MQ=255)

                           ttGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATagaag                                                           >  2:4077086/1‑93 (MQ=255)

                           ttGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATagaag                                                           <  1:1015645/93‑1 (MQ=255)

                             ggCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAgaagaa                                                         <  2:4133084/93‑1 (MQ=255)

                             ggCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAgaagaa                                                         <  2:3662095/93‑1 (MQ=255)

                              gCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGaagaag                                                        <  1:1094117/93‑1 (MQ=255)

                              gCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGaagaag                                                        >  1:3336393/1‑93 (MQ=255)

                               cATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGc                                                       <  2:2361128/93‑1 (MQ=255)

                                aTTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCg                                                      <  1:4608973/93‑1 (MQ=255)

                                 ttATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCgg                                                     >  2:2383971/1‑93 (MQ=255)

                                 ttATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCgg                                                     <  1:3399423/93‑1 (MQ=255)

                                   aTCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGAt                                                   <  2:1773063/93‑1 (MQ=255)

                                     cgcgCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCa                                                 <  2:2755981/93‑1 (MQ=255)

                                     cgcgCTGCGCATCACGGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCCCTCTTTTTCATCTTCGGGTCCCTTATAGAAGAAGCGGATCa                                                 <  2:3984043/93‑1 (MQ=255)

                                         cTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTc                                             <  2:1967605/93‑1 (MQ=255)

                                           tcgcATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGc                                           <  2:2338453/92‑1 (MQ=255)

                                           gcgcATCACTGAGGGCCTGATCAAACGAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGc                                           <  2:3793600/93‑1 (MQ=255)

                                            cgcATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCt                                          >  2:1503930/1‑93 (MQ=255)

                                            cgcATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCt                                          <  2:65041/93‑1 (MQ=255)

                                              cATCACTGAGGGCCTGAACAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTgg                                        >  2:3218100/1‑93 (MQ=255)

                                               aTCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGc                                       >  1:34386/1‑93 (MQ=255)

                                               aTCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGc                                       >  1:691566/1‑93 (MQ=255)

                                                tCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCt                                      >  2:3701733/1‑93 (MQ=255)

                                                  aCTGAGGTCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTg                                    <  2:535060/93‑1 (MQ=255)

                                                  aCTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTg                                    >  1:1907063/1‑93 (MQ=255)

                                                   cTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTgg                                   <  2:765498/93‑1 (MQ=255)

                                                    tGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTggg                                  >  2:438267/1‑93 (MQ=255)

                                                    tGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTggg                                  >  2:2993733/1‑93 (MQ=255)

                                                       gggCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCgg                               <  1:3125787/93‑1 (MQ=255)

                                                        ggCCTGATCAAACCAGACAATGCGGTAGTAGAACTCTTTTCCATTTGCGGGTCCCTTATAGACGAACCGGATCACTTCGCTGGCTTGGGCggg                              >  1:1505284/1‑93 (MQ=255)

                                                         gCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGt                             >  1:1692772/1‑93 (MQ=255)

                                                          ccTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTa                            <  2:4326269/93‑1 (MQ=255)

                                                          ccTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTa                            <  1:450561/93‑1 (MQ=255)

                                                           cTGATCAAACCAGTCAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTag                           >  2:4603431/1‑93 (MQ=255)

                                                            tGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagc                          <  2:4055080/93‑1 (MQ=255)

                                                            tGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGGagc                          >  2:4522821/1‑93 (MQ=255)

                                                             gATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagca                         >  1:4521708/1‑93 (MQ=255)

                                                             gATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagca                         <  2:1044207/93‑1 (MQ=255)

                                                                cAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAgcagca                      >  2:3075674/1‑93 (MQ=255)

                                                                 aaaCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCaa                     <  2:2687399/93‑1 (MQ=255)

                                                                 aaaCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCaa                     <  1:4537362/93‑1 (MQ=255)

                                                                 aaaCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGAAGCaa                     >  1:1127123/1‑93 (MQ=255)

                                                                  aaCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAg                    >  2:2614436/1‑93 (MQ=255)

                                                                   aCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCTGGTAGCAGGAAGc                   >  1:1454567/1‑93 (MQ=255)

                                                                   aCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGc                   <  2:3751898/93‑1 (MQ=255)

                                                                   aCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGc                   >  1:5393/1‑93 (MQ=255)

                                                                   aCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGc                   <  2:3546998/93‑1 (MQ=255)

                                                                    ccAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCt                  >  2:1676493/1‑93 (MQ=255)

                                                                     cAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTg                 <  1:170100/93‑1 (MQ=255)

                                                                      aGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTgg                <  2:136585/93‑1 (MQ=255)

                                                                        aCAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCg              >  1:2688862/1‑93 (MQ=255)

                                                                         cAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCgg             >  1:1518638/1‑93 (MQ=255)

                                                                         cAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCgg             >  2:4480145/1‑93 (MQ=255)

                                                                          aaTGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCggg            >  1:3894783/1‑93 (MQ=255)

                                                                          aaTGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCggg            >  1:1860685/1‑93 (MQ=255)

                                                                          aaTGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCggg            <  2:1930910/93‑1 (MQ=255)

                                                                          aaTGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCACTGGCTTGGGCGGGGAGCAGCAAGCTGGCggg            >  2:1984059/1‑93 (MQ=255)

                                                                           aTGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGga           >  1:3881113/1‑93 (MQ=255)

                                                                           aTGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGga           >  2:2018698/1‑93 (MQ=255)

                                                                                gtagtaGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagtgag      <  1:3227201/93‑1 (MQ=255)

                                                                                   gtagCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAg   >  1:1014078/1‑93 (MQ=255)

                                                                                   gtagCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAg   <  1:3805294/93‑1 (MQ=255)

                                                                                   gtagCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAg   >  2:3825209/1‑93 (MQ=255)

                                                                                    tagCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGc  <  2:4159245/93‑1 (MQ=255)

GCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGC  >  NC_000913/308847‑309023

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