BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,567,293	A→G	L554L (TTA→CTA)	glgX ←	GlgX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,567,293	0	A	G	97.5%	724.2 / 4.6	239	L554L (TTA→CTA)	glgX	GlgX
Reads supporting (aligned to +/- strand): ref base A (1/5); new base G (117/116); total (118/121)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 250 total aligned reads displayed.

ATATCGATTTAGCCAACGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAACGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCTGACAGTAGGCATTGTTATTGC  >  NC_000913/3567204‑3567382

 tatCGATTTAGCCAACGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGcgcg                                                                                       >  1:3682150/1‑93 (MQ=255)

   tCGATTTAGCCAACGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCgcggc                                                                                     <  1:2120744/93‑1 (MQ=255)

     gATTTAGCCAACGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGcggcgg                                                                                   >  2:1104605/1‑93 (MQ=255)

      aTTTAGCCAACGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGt                                                                                  >  2:3109414/1‑93 (MQ=255)

      aTTTAGCCAACGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGt                                                                                  <  1:2099221/93‑1 (MQ=255)

      aTTTAGCCAACGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGAGCGGCGGt                                                                                  >  1:3524620/1‑93 (MQ=255)

        ttAGCCAACGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAACGCGGCGGTaa                                                                                <  1:159649/93‑1 (MQ=255)

         tAGCCAACGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTaaa                                                                               <  2:3871027/93‑1 (MQ=255)

         tAGCCAACGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTaaa                                                                               <  2:4565034/93‑1 (MQ=255)

          aGCCAACGGACATTGCCGTCGCCTTCTTCCCACCAGTGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAAt                                                                              <  2:4171637/93‑1 (MQ=255)

            ccAACGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGc                                                                            <  2:3345141/93‑1 (MQ=255)

            ccAACGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGc                                                                            >  2:3264032/1‑93 (MQ=255)

             cAACGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGCTTAGCGCGGCGGTAAATGCg                                                                           <  2:3762269/93‑1 (MQ=255)

             cAACGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAACTGCg                                                                           >  2:3466686/1‑93 (MQ=255)

             cAACGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCg                                                                           <  2:1541533/93‑1 (MQ=255)

              aaCGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCgg                                                                          >  2:197289/1‑93 (MQ=255)

              aaCGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCgg                                                                          >  2:290683/1‑93 (MQ=255)

               gcGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGt                                                                         <  2:4436224/92‑1 (MQ=255)

               aCGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGt                                                                         >  2:2362540/1‑93 (MQ=255)

               aCGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGt                                                                         <  1:3045947/93‑1 (MQ=255)

               aCGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGt                                                                         >  1:1639534/1‑93 (MQ=255)

                  gACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTaa                                                                      >  1:361097/1‑93 (MQ=255)

                   aCATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTaaa                                                                     <  1:1439213/93‑1 (MQ=255)

                    cATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAAc                                                                    >  1:1805683/1‑93 (MQ=255)

                    cATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAAc                                                                    >  2:2575684/1‑93 (MQ=255)

                     aTTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAAcc                                                                   <  2:2739824/93‑1 (MQ=255)

                        gCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAGGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACt                                                                <  1:4021815/93‑1 (MQ=255)

                        gCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACt                                                                <  2:4093487/93‑1 (MQ=255)

                        gCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAAACACt                                                                >  2:806935/1‑93 (MQ=255)

                         ccGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTg                                                               <  2:3727411/93‑1 (MQ=255)

                         ccGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTg                                                               <  1:104842/93‑1 (MQ=255)

                           gTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCt                                                             >  2:894422/1‑93 (MQ=255)

                           gTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCt                                                             >  2:1161429/1‑93 (MQ=255)

                           gTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCt                                                             >  1:2349826/1‑93 (MQ=255)

                           gTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCt                                                             >  1:2016902/1‑93 (MQ=255)

                              gCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGGAAATGCGGTTAAACCACTGCTTGc                                                          >  1:446438/1‑93 (MQ=255)

                               ccTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGcc                                                         >  1:3096676/1‑93 (MQ=255)

                                cttcttCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCt                                                        >  1:1032424/1‑93 (MQ=255)

                                 ttcttcCCACCAGCGATTCTCCACCAAATCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTg                                                       >  1:3323864/1‑93 (MQ=255)

                                 ctcttcCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTg                                                       <  2:3657604/92‑1 (MQ=255)

                                  tcttcCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGc                                                      >  1:3574789/1‑93 (MQ=255)

                                  tcttcCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGc                                                      >  2:3570955/1‑93 (MQ=255)

                                    ttcCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGa                                                    >  2:2454607/1‑93 (MQ=255)

                                    ttcCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGa                                                    >  1:3654574/1‑93 (MQ=255)

                                      cccACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGAcc                                                  <  1:1948973/93‑1 (MQ=255)

                                        caccaGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAg                                                <  1:23049/93‑1 (MQ=255)

                                        caccaGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAg                                                <  1:962962/93‑1 (MQ=255)

                                         cccaGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGt                                               <  2:3151656/92‑1 (MQ=255)

                                          ccaGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTc                                              >  1:2039579/1‑93 (MQ=255)

                                          ccaGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTc                                              >  1:4358069/1‑93 (MQ=255)

                                           caGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTcc                                             <  1:3211957/93‑1 (MQ=255)

                                           caGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAACGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTcc                                             <  1:3984312/93‑1 (MQ=255)

                                             gCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCaa                                           >  2:3820608/1‑93 (MQ=255)

                                              cGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAAc                                          <  1:778539/93‑1 (MQ=255)

                                               gATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAAca                                         >  2:1873748/1‑92 (MQ=255)

                                                aTTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCa                                        >  2:3539564/1‑93 (MQ=255)

                                                aTTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCa                                        <  1:4374382/93‑1 (MQ=255)

                                                  tctcCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAgg                                      <  1:4176094/93‑1 (MQ=255)

                                                  tctcCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAgg                                      <  2:2064905/93‑1 (MQ=255)

                                                    tcCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGtt                                    >  2:4080753/1‑93 (MQ=255)

                                                      caccaAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTaa                                  <  1:2481397/93‑1 (MQ=255)

                                                      caccaAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTaa                                  <  2:2191339/93‑1 (MQ=255)

                                                      caccaAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTaa                                  <  1:2259820/93‑1 (MQ=255)

                                                      caccaAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTaa                                  >  1:2200964/1‑93 (MQ=255)

                                                      caccaAAGCGGGAATGCGCTTGCGCAGATGGATTAACGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTaa                                  <  1:3837028/93‑1 (MQ=255)

                                                       accaAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAAt                                 >  1:474118/1‑93 (MQ=255)

                                                       accaAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAAt                                 <  1:3494320/93‑1 (MQ=255)

                                                        ccaAAGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAAtt                                <  2:2815203/93‑1 (MQ=255)

                                                          aaaGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTgg                              >  2:2229289/1‑93 (MQ=255)

                                                          aaaGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGCTAATTgg                              >  1:2759895/1‑93 (MQ=255)

                                                           aaGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGt                             <  2:1511568/93‑1 (MQ=255)

                                                           aaGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGt                             <  2:1609698/93‑1 (MQ=255)

                                                            aGCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGtt                            >  2:1306461/1‑93 (MQ=255)

                                                             gCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTa                           <  2:3320830/93‑1 (MQ=255)

                                                             gCGGGAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTa                           <  2:4215927/93‑1 (MQ=255)

                                                               gggAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATc                         >  2:1444225/1‑93 (MQ=255)

                                                                ggAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATcc                        <  1:2816374/93‑1 (MQ=255)

                                                                ggAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATcc                        <  1:2889280/93‑1 (MQ=255)

                                                                ggAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATcc                        >  1:4494783/1‑93 (MQ=255)

                                                                ggAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATcc                        <  2:745290/93‑1 (MQ=255)

                                                                 gAATGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGTTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCt                       <  2:1065680/93‑1 (MQ=255)

                                                                   aTGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCTGa                     >  2:900044/1‑93 (MQ=255)

                                                                   aTGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCTGa                     <  2:1522355/93‑1 (MQ=255)

                                                                   aTGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCTGa                     <  1:4551477/93‑1 (MQ=255)

                                                                    tGCGCTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCTGAc                    <  2:3193644/93‑1 (MQ=255)

                                                                      cgcTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCTGACAg                  <  2:115269/93‑1 (MQ=255)

                                                                      cgcTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCTGACAg                  <  2:540931/93‑1 (MQ=255)

                                                                       gcTTGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCTGACAga                 >  1:1280814/1‑92 (MQ=255)

                                                                         ttGCGCAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCTGACAGTAg               <  1:899282/93‑1 (MQ=255)

                                                                            cgcAGATGGATTAACGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCTGACAGTAGGCa            <  1:2532358/93‑1 (MQ=255)

                                                                             gcAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCTGACAGTAGGCAt           <  1:264682/93‑1 (MQ=255)

                                                                              cAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCTGACAGTAGGCAtt          <  2:298937/93‑1 (MQ=255)

                                                                              cAGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCTGACAGTAGGCAtt          <  1:838875/93‑1 (MQ=255)

                                                                               aGATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATACTGACAGTAGGCATTg         >  2:4638314/1‑93 (MQ=255)

                                                                                gATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCTGACAGTAGGCATTGt        >  1:3420744/1‑93 (MQ=255)

                                                                                gATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCTGACAGTAGGCATTGt        <  1:4114857/93‑1 (MQ=255)

                                                                                gATGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCTGACAGTAGGCATTGt        >  2:1702169/1‑93 (MQ=255)

                                                                                 aTGGATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCTGACAGTAGGCATTGtt       >  2:1788715/1‑93 (MQ=255)

                                                                                   ggATTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCTGACAGTAGGCATTGTTAt     <  1:1489090/93‑1 (MQ=255)

                                                                                     aTTAGCGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCTGACAGTAGGCATTGTTATTg   >  1:1353251/1‑93 (MQ=255)

ATATCGATTTAGCCAACGGACATTGCCGTCGCCTTCTTCCCACCAGCGATTCTCCACCAAAGCGGGAATGCGCTTGCGCAGATGGATTAACGCGGCGGTAAATGCGGTTAAACCACTGCTTGCCTGCGACCAGTCCAACCAGGTTAATTGGTTATCCTGACAGTAGGCATTGTTATTGC  >  NC_000913/3567204‑3567382

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