BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,579,447	T→C	S598P (TCG→CCG)	tktB →	transketolase 2 isozyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,579,447	0	T	C	99.6%	774.4 / ‑4.9	244	S598P (TCG→CCG)	tktB	transketolase 2 isozyme
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (125/118); total (125/119)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 250 total aligned reads displayed.

GAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGC  >  NC_000913/2579357‑2579535

gagaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCg                                                                                         >  1:143272/1‑93 (MQ=255)

gagaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCg                                                                                         >  2:958120/1‑93 (MQ=255)

gagaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCg                                                                                         >  2:8670476/1‑93 (MQ=255)

  gaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgt                                                                                       >  1:3566441/1‑93 (MQ=255)

   aaaTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgtg                                                                                      <  2:2323658/93‑1 (MQ=255)

    aaTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgtgt                                                                                     <  2:1700546/93‑1 (MQ=255)

    aaTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgtgt                                                                                     >  2:8749649/1‑93 (MQ=255)

           aGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTc                                                                              <  1:5876132/93‑1 (MQ=255)

            ggAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCt                                                                             >  2:3224506/1‑93 (MQ=255)

            ggAGAAGGACGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGGCTATCGGGAGCCGGTGTTGGCTTCt                                                                             >  2:8244473/1‑93 (MQ=255)

             gagaAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGGCTTCTa                                                                            >  2:7995788/1‑93 (MQ=255)

             gagaAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTa                                                                            >  1:1701603/1‑93 (MQ=255)

              agaAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTaa                                                                           >  2:6640578/1‑93 (MQ=255)

              agaAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTaa                                                                           <  1:2027166/93‑1 (MQ=255)

              agaAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTaa                                                                           <  2:5272154/93‑1 (MQ=255)

              agaAGGTAGCAATGTACGAGTAGTTTCCCTGCCCTCGACAGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTaa                                                                           <  2:838642/93‑1 (MQ=255)

                aaGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACg                                                                         <  2:4133456/93‑1 (MQ=255)

                 aGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGt                                                                        >  2:802050/1‑93 (MQ=255)

                 aGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGt                                                                        >  1:5707012/1‑93 (MQ=255)

                  ggTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGtt                                                                       >  1:3006596/1‑93 (MQ=255)

                   gTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTg                                                                      >  2:9908278/1‑93 (MQ=255)

                       cAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGc                                                                  >  2:2527709/1‑93 (MQ=255)

                        aaTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCt                                                                 >  2:3618676/1‑93 (MQ=255)

                        aaTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCt                                                                 <  1:9490409/93‑1 (MQ=255)

                         aTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTc                                                                <  1:3072330/93‑1 (MQ=255)

                         aTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTc                                                                <  2:2256/93‑1 (MQ=255)

                          tGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTcc                                                               >  2:6143298/1‑92 (MQ=255)

                           gTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGCATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTcgg                                                              >  1:1911803/1‑92 (MQ=255)

                            tACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTTCGGCTcgc‑g                                                            >  1:702170/1‑93 (MQ=255)

                            tACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTcgc‑g                                                            >  2:3002965/1‑93 (MQ=255)

                             aCGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑Gt                                                           <  1:7873039/93‑1 (MQ=255)

                              cgcgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTg                                                          <  2:1964637/93‑1 (MQ=255)

                              cgcgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTg                                                          <  1:9729763/93‑1 (MQ=255)

                               gcgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑CTgg                                                         >  1:9612204/1‑93 (MQ=255)

                                  tAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCgg                                                      >  2:5062007/1‑93 (MQ=255)

                                  tAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCgg                                                      <  1:1012380/93‑1 (MQ=255)

                                     tttCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGgggg                                                   >  1:430450/1‑93 (MQ=255)

                                     tttCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCCC‑GTGGCGGTgg                                                   >  1:4087892/1‑93 (MQ=255)

                                       tCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGaa                                                 <  1:6402720/93‑1 (MQ=255)

                                       tCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGaa                                                 <  2:6943205/93‑1 (MQ=255)

                                        cccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGagg                                                >  1:3191240/1‑91 (MQ=255)

                                        cccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAg                                                <  1:7188520/93‑1 (MQ=255)

                                          cttcccTCGACCGATATTTTCCACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCa                                              <  1:2811522/90‑1 (MQ=255)

                                          cTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCa                                              <  1:7297711/93‑1 (MQ=255)

                                           tGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAg                                             <  2:5780955/93‑1 (MQ=255)

                                           tGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAg                                             <  1:4604304/93‑1 (MQ=255)

                                           tGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGAGTTGCCTTCTAACGTTGCGGCTCGC‑GCGGGGGGGGGAGCAg                                             >  1:9618108/1‑93 (MQ=255)

                                            gCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGGAGCAgg                                            >  2:5938687/1‑93 (MQ=255)

                                             cccTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGt                                           <  1:4184491/93‑1 (MQ=255)

                                              ccTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTa                                          <  1:1051363/93‑1 (MQ=255)

                                              ccTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTa                                          <  2:3781069/93‑1 (MQ=255)

                                               cTCGACCGTTATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCTCGGTGGGGGT‑GAAGCAGCTAt                                         >  2:8630047/1‑93 (MQ=255)

                                               cTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGGGGGAGCAGGGAt                                         >  1:4729622/1‑93 (MQ=255)

                                                tCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGGAtt                                        >  1:2316811/1‑93 (MQ=255)

                                                tCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGGGGGAGGAGGGAtt                                        >  2:9765928/1‑93 (MQ=255)

                                                   aCCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGGAGCAGGTATTGcc                                     >  2:6411830/1‑93 (MQ=255)

                                                   aCCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGGGGGAGGAGGTATTGcc                                     >  2:5417669/1‑93 (MQ=255)

                                                   aCCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGGGGGAGCAGGGATTGcc                                     >  2:430294/1‑93 (MQ=255)

                                                    ccGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGGGGGAGCAGGGATTGCCg                                    >  1:3246184/1‑93 (MQ=255)

                                                     cGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCgg                                   >  1:525991/1‑92 (MQ=255)

                                                      gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGGAGCAGGTATTGCCGAt                                  >  1:8412336/1‑93 (MQ=255)

                                                       atatTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGGGGGAGCAGGTATTGCCGAtt                                 >  2:7405299/1‑93 (MQ=255)

                                                        tatTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTa                                <  2:2871580/93‑1 (MQ=255)

                                                        tatTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTa                                <  1:9269911/93‑1 (MQ=255)

                                                        tatTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGGGGAAGCAGGTATTGCCGATTa                                >  1:9138076/1‑93 (MQ=255)

                                                          ttttCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGGAGCAGGTATTGCCGATTACt                              >  1:1441257/1‑93 (MQ=255)

                                                          ttttCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGAGGAAGCAGGTATTGCCGATTACt                              >  2:4974659/1‑93 (MQ=255)

                                                             tCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGt                           <  2:7685930/93‑1 (MQ=255)

                                                              cGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTa                          <  1:6006197/93‑1 (MQ=255)

                                                                aCGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACa                        <  2:6210538/93‑1 (MQ=255)

                                                                  gCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAg                      <  1:7280608/93‑1 (MQ=255)

                                                                  gCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAg                      <  2:118475/93‑1 (MQ=255)

                                                                   cccAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGt                     >  2:1152556/1‑93 (MQ=255)

                                                                    ccAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTa                    <  2:1468444/93‑1 (MQ=255)

                                                                    ccAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTa                    <  1:857090/93‑1 (MQ=255)

                                                                     cAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGGAt                   >  1:3812446/1‑93 (MQ=255)

                                                                     cAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑CGGGCGGTGGGAGCAGGCGTTGCCGATTTCTGGAACCAGTAt                   >  1:2566924/1‑93 (MQ=255)

                                                                       ggATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGGATTGCCGATTACTGGTACAAGTATGt                 >  2:7692933/1‑93 (MQ=255)

                                                                       ggATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGGATTGCCGATTACTGGTACAAGTATGt                 >  1:6784554/1‑93 (MQ=255)

                                                                       ggATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGGATTGCCGATTACTGGGACAAGTATGt                 >  1:9378196/1‑93 (MQ=255)

                                                                        gATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGtt                <  1:6297646/93‑1 (MQ=255)

                                                                        gATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGtt                <  2:2498664/93‑1 (MQ=255)

                                                                         aTGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTg               <  1:4671762/93‑1 (MQ=255)

                                                                         aTGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGGATTGCCGATTACTGGTACAAGTATGTTg               >  1:9849898/1‑93 (MQ=255)

                                                                          tGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTgg              <  1:3960260/93‑1 (MQ=255)

                                                                          tGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGGATTGCCGATTACTGGTACAAGTATGTTgg              >  1:2518005/1‑93 (MQ=255)

                                                                           gAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGt             <  1:8882530/93‑1 (MQ=255)

                                                                            aGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTc            >  1:9033535/1‑93 (MQ=255)

                                                                               aaTATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTTGTCTGa         >  2:3156177/1‑93 (MQ=255)

                                                                               aaTATCGGGAGCCGGTGGTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGa         <  1:7074950/93‑1 (MQ=255)

                                                                                atatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaa        <  1:4512826/93‑1 (MQ=255)

                                                                                atatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaa        <  1:8279138/93‑1 (MQ=255)

                                                                                 tatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaaa       >  2:7033245/1‑93 (MQ=255)

                                                                                 tatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaaa       >  1:2158620/1‑93 (MQ=255)

                                                                                 tatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGGATGTTGGGCTGaaa       >  2:1376617/1‑93 (MQ=255)

                                                                                 tatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGGATTGCCGATTACTGGTACAAGGATGTTGGTCTGaaa       >  1:59485/1‑93 (MQ=255)

                                                                                 tatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGGGGGAGCAGGTATTGCCGGTTACTTGTTCCAGCATGGTGGGCTGaaa       >  2:7479234/1‑93 (MQ=255)

                                                                                  atCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAg      <  2:8744229/93‑1 (MQ=255)

                                                                                  atCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAg      <  2:9320714/93‑1 (MQ=255)

                                                                                      ggAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGc  <  1:5075063/93‑1 (MQ=255)

GAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGC‑GTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGC  >  NC_000913/2579357‑2579535

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